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Abstract 

Irradiation  effects  on  the  wide-bandgap  semiconductor  material  GaN  are  of 
interest  to  the  USAF  due  to  this  material’s  applicability  for  a  wide  range  of  on-orbit  uses. 
Irradiation  is  also  a  valuable  tool  in  analyzing  the  damage  and  defect  formation  dynamics 
of  the  material  which  is  of  great  use  in  determining  and  correcting  deficiencies  in 
material  growth  processes. 

GaN  samples  representing  several  different  growth  methods  and  doping  profiles 
were  irradiated  by  1.0  -  1.5  MeV  electron  beams  to  induce  defects  such  as  vacancies  and 
interstitial  atoms  in  the  material.  Following  irradiation,  the  samples  were  analyzed  using 
Electron  Paramagnetic  Spin  Resonance  (EPR)  spectroscopy  to  determine  effects  of  the 
irradiation,  with  particular  attention  to  the  effects  of  irradiation  on  the  shallow  donors  in 
the  material. 

EPR  allows  the  direct  inspection  of  paramagnetic  impurity  or  defect  sites, 
providing  information  on  site  density  as  well  as,  in  some  cases,  identity.  Samples 
subjected  to  EPR  analysis  prior  to  irradiation  show  a  strong  signal  attributed  to  the 
shallow  donor  in  GaN  that  is  without  resolved  hyperfine  structure.  Following  low- 
temperature  irradiation  with  1.0  MeV  electrons  to  a  total  dose  of  approximately  40  MRad 
(GaN),  the  same  samples  showed  a  marked  decrease  in  the  shallow  donor  signal  and  the 
introduction  of  a  broad  EPR  signal  with  resolved  hyperfine  structure. 

The  decrease  in  the  shallow  donor  signals  is  attributed  to  the  formation  of  a  deep 
band-gap  complex  of  the  nitrogen  interstitial  with  the  shallow  donor  sites,  perhaps  in 


IV 


combination  with  compensation  by  defect-related  centers.  The  damage  constant  for  the 
formation  of  this  process  is  calculated  as  0. 14 cm“' .  This  decrease  in  the  shallow  donor 
concentration  in  undoped  GaN  demonstrates  that  the  native  shallow  donor  in  n-type  GaN 
cannot  be  the  nitrogen  vacancy,  as  has  been  previously  suggested. 

A  broad  EPR  signal  imparted  by  low-temperature  1 .0  MeV  electron  irradiation  is 
identified  as  a  gallium  interstitial  by  modeling  of  the  hyperfine  spectrum.  The  Breit-Rabi 
formula  is  used  with  parameters  characteristic  of  the  two  naturally  occurring  gallium 
isotopes  to  reproduce  the  observed  spectrum.  This  is  the  first  observation  of  the  gallium 
interstitial  by  EPR,  as  well  as  the  first  simultaneous  resolvable  measurement  of  nitrogen 
and  gallium  sublattice  damage  in  a  single  sample  following  irradiation. 
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I.  Introduction 


The  rapidly  developing  field  of  wide-bandgap  semiconductor  research,  and  in 
particular  Ill-Nitride  research,  is  of  current  interest  to  the  United  States  Air  Force  due  to 
the  potential  applications  of  devices  based  upon  these  materials  in  militarily  important 
environments  and  roles.  Gallium  Nitride  (GaN)  specifically  exhibits  material  properties 
that  make  it  well  suited  to  applications  where  high-temperature  operation  [1],  efficient 
high-frequency  operation  [2],  and  high-power  density  properties  are  required  [3].  In 
addition  to  these  unique  properties,  the  wide  bandgap  of  GaN  is  well  suited  for  use  in 
short  wavelength  opto-electronic  devices  [4] . 

The  development  of  high-power  density,  high  frequency  devices  capable  of 
operating  in  high-temperature  environments  is  readily  seen  to  be  of  military  interest  due 
to  the  potential  for  rugged,  all-weather  communications  systems,  on-orbit  processing  and 
communications  capabilities,  and  the  inclusion  of  solid-state  devices  in  aerospace 
applications  where  temperature  and  environmental  conditions  have  previously  limited 
their  use.  With  these  properties,  GaN  and  related  wide-bandgap  semiconductors  would 
seem  to  be  a  panacea  for  shortfalls  in  military  and  commercial  systems  design;  however, 
growth  problems  have  severely  limited  the  production  of  efficient  devices  [5]  and  have 
slowed  the  maturity  of  GaN  based  systems  to  the  point  of  commercial  availability  and 
military  feasibility.  Additionally,  previous  research  has  determined  that  native  GaN 
materials  are  relatively  “radiation  hard”,  or  resistant  to  radiation  effects  [6]  in  comparison 
to  other  materials. 
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Radiation  effects  testing  is  thus  important  for  two  distinct  reasons:  1)  the 
interactions  of  radiation  with  materials  sheds  light  upon  the  formation  and  dynamics  of 
the  types  of  defects  that  have  limited  device  development,  and  2)  the  supposed 
application  of  these  materials  to  military  (specifically  on-orbit)  applications  will  be  aided 
by  the  existence  of  information  on  device  and  material  radiation  hardness.  As  devices 
cross  the  threshold  from  laboratory  experiments  to  commercial  production  and  military 
usefulness,  as  is  happening  at  present,  it  is  appropriate  to  carry  out  radiation  effects 
studies  to  fulfill  these  two  purposes. 

This  dissertation  and  accompanying  research  efforts  were  designed  to  fill  this 
niche,  offering  information  on  the  basic  radiation  interactions  in  GaN  material.  This  will 
hopefully  serve  as  a  solid  basis  for  the  development  of  device  testing  research  efforts 
(which  the  Department  of  Defense  and  Air  Force  Office  of  Scientific  Research  are 
currently  sponsoring),  as  well  as  providing  answers  to  some  of  the  more  contentious 
questions  that  have  developed  in  the  GaN  growth  community.  To  this  end,  we  have 
applied  a  novel  experimental  methodology  to  the  study  of  irradiation  effects  on  GaN 
materials:  the  application  of  Electron  Paramagnetic  Resonance  (EPR)  Spectroscopy  in 
concert  with  Hall  effect  studies  of  irradiated  materials.  The  use  of  these  methods  in  a 
“symbiotic”  fashion  offers  more  information  than  the  application  of  either  could  produce 
independently. 

This  work  provides  conclusive  answers  some  outstanding  questions  concerning 
the  nature  of  growth  problems  in  GaN  layers,  identifies  some  of  the  direct  products  of 
irradiation  in  these  materials,  and  directly  measures  the  impact  of  these  products  on 
material  parameters  which  will  directly  impact  device  performance.  As  such,  it  is  hoped 
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that  this  research  shall  provide  a  bridge  to  span  the  divide  between  material  development 
and  radiation  effects  testing  of  viable,  militarily  useful  devices. 
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II.  Background 


Gallium  Nitride  (GaN)  Applications  and  Testing 

Gallium  Nitride  (GaN),  along  with  the  related  compounds  Aluminum  Nitride 
(AIN)  and  Aluminum  Gallium  Nitride  (AlGaN)  offer  advantages  and  challenges  in  the 
production  of  novel  semiconductor  devices.  GaN  displays  unique  properties  among  the 
III-V  compounds  as  well  as  difficulties  both  in  doping  and  crystal  growth. 

GaN  is  of  interest  to  the  United  States  Air  Force  (USAF)  and  US  Department  of 
Defense  (DoD)  because  of  the  promise  of  high-power,  high-temperature,  and  high- 
frequency  semiconductor  applications  [1],[2],[3]  as  well  as  for  its  unique  optical  uses  in 
short  wavelength  diode  lasers  and  light  emitting  diodes  (LEDs)  [4]  .  These  unique  areas 
of  operational  promise  are  quite  likely  to  be  employed  in  radiation  exposure  situations 
such  as  on-orbit  communications  hardware  which  is  exposed  to  the  space  radiation 
environment,  high-temperature  devices  for  use  in  reactor  monitoring  and  control,  or  high- 
power  switches  which  may  one  day  replace  existing  silicon-based  switches  in  military 
and  civilian  power  grids,  bringing  the  necessity  of  radiation  hardening  against  man-made 
sources  or  solar  disruptions.  GaN  and  its  related  compounds  are  already  coming  into 
widespread  commercial  use  for  many  of  these  optical  applications  [7]  and  are  beginning 
to  be  produced  commercially  for  high-power  density  applications  [8].  Given  the 
radiation  environments  in  which  GaN-based  devices  will  eventually  be  tasked  to  perform, 
it  is  prudent  for  radiation  studies  of  the  underlying  materials  to  be  undertaken  at  the 
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current  stage  of  development  to  aid  in  the  modeling  and  simulation  of  GaN  based  devices 
in  these  environments. 

The  study  of  energetic  electron  interactions  in  GaN  will  form  a  baseline  from 
which  future  radiation  effects  studies  may  proceed.  In  addition  to  the  emphasis  on  device 
performance  in  radiation  environments,  the  radiation  testing  of  GaN  materials  provides 
valuable  insight  into  material  defect  formation  processes  and  dynamics  -  insight  that  can 
be  used  to  improve  growth  processes  and  produce  better  materials.  Irradiation  of 
semiconductor  materials  is  a  valued  means  of  producing  the  same  types  of  point  defects 
in  GaN  [9]  that  material  growth  processes  seek  to  minimize.  Irradiation-induced  defects 
can  be  used  to  study  the  dynamics  and  interactions  of  defects  with  band-gap  states  to 
include  defect  donors,  acceptors,  and  trap  states;  compensation  rates  and  complexation 
behavior  of  defects  with  band-gap  states  can  be  determined  as  well. 

Physical  Properties  of  GaN 

GaN  is  a  wide-bandgap  semiconductor  material  with  a  300K  direct  bandgap  of 
3.43  eV  (361.2  nm)  in  the  thermodynamically  stable  wurtzite  crystalline  phase.  GaN  is 
recognized  as  the  most  widely  studied  Ill-Nitride  compound,  but  even  at  the  current  time, 
much  research  remains  to  be  done  on  the  properties  of  GaN.  Some  selected  physical 
properties  of  wurtzite  GaN  are  given  in  Table  1. 
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Table  1.  Selected  Physical  Properties  of  GaN 


Property  (units) 

Symbol 

Value 

Molecular  weight  (gm/mol)  [10] 

W 

83.728 

Density  (gm/cm^)  [11] 

P 

6.15 

Thermal  Conductivity  (W/cm-K)  [10] 

K 

1.3 

Specific  Heat  (cal/mol  K)  [10] 

Cp 

9.1  +(2.15- 10'^  T) 

Static  Dielectric  Constant  [12] 

ep 

9.0 

High  Frequency  Dielectric  Constant  [13] 

£00 

5.35 

Electron  mobility,  bulk  (cm  /V  sec)  [12] 

Pe 

1000 

Hole  mobility,  bulk  (cm^/V  sec)  [12] 

Ph 

300 

Debye  Fength  @  300K,  Nd=10^^  cm'"* 

(nm) 

Xd 

3.586 

Index  of  Refraction  @  I  pm  [10] 

n 

2.35 

Electron  affinity  (eV)  [14] 

X 

4.1 

The  bandgap  of  the  wurtzite  phase  varies  with  temperature  as  given  in  the 
empirical  Varshni  equation 

aT^ 

£  (T)=E(0)-- - - 

(r+/?) 

where,  for  bulk  GaN,  a  =  -l .1  ■  10  eV jK  and  [5  =  600  K  and  E(0)  is  the  fundamental 
bandgap  energy.  The  most  widely  accepted  value  of  E(0)  is  3.47  eV  [14].  The 
temperature  dependant  bandgap  for  the  wurtzite  phase  of  bulk  GaN  is  shown  in  Figure  1 
as  generated  by  the  analysis  of  Equation  (1).  In  addition  to  its  large  bandgap  energy, 
GaN  is  a  direct  bandgap  semiconductor.  Details  of  the  300K  bandgap  near  the  T 
symmetry  point  are  shown  in  Figure  2,  demonstrating  the  direct-bandgap  nature  of  GaN 
materials. 
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Figure  1 .  Estimated  fundamental  bandgap  energy  dependence  upon  temperature 
in  wurtzite  phase  GaN.  Dotted  vertical  lines  represent  liquid  He,  liquid  N,  and  room 
temperature  points.  Data  shown  here  was  generated  by  equation  1. 


Figure  2.  Wurtzite  GaN  band- structure  as  computed  by  Suzuki  [15],  et  al  via  a  plane- 
wave  pseudopotential  calculation  method.  All  energy  values  are  at  300°K. 
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A  metastable  zinc-blende  phase  also  occurs  with  a  300K  bandgap  of 
approximately  3.2  eV  (387.5  nm)  [16].  Although  the  zinc-blende  phase  can  be  formed  by 
epitaxial  stabilization,  the  wurtzite  crystalline  structure  is  the  dominant  structure  in 
device  construction,  and  is  the  most  thoroughly  studied  of  the  two  forms.  Models  of  the 
wurtzite  structure  in  GaN  are  shown  in  Figure  3  for  two  views  of  interest.  All  further 
references  to  GaN  refer  to  the  wurtzite  phase  unless  otherwise  noted. 


Figure  3.  Calculated  atomic  locations  for  wurtzite  crystal  structure  of  GaN.  The  graph 
on  the  left  illustrates  the  orientation  of  the  (  000 1)  direction,  with  Ga  represented  by  black 
and  N  by  white  balls.  Axis  units  are  in  angstroms.  Plot  on  the  right  is  presented  looking 
down  the  (OOO 1 )  direction  for  a  Ga  face  surface. 
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Wurtzite  GaN  has  lattice  constants  a  =  3. 189  A  and  c  =  5.185A  with  thermal 

coefficients  given  by  —  =  5.59  ■  10“®  K  and  —  =  3.17  ■  10“®  K  ^  [10].  These  values  are 

a  c 

important  in  evaluating  the  effects  of  various  substrate  materials  on  GaN  film  quality  and 
device  performance. 

The  value  of  the  electron  effective  mass  is  commonly  accepted  to  be 
0.22  •  although  the  effective  mass  is  not  presently  determined  to  great  precision. 

Heavy  hole  masses  are  known  to  be  greater  than 0.6  -  mg,  where  mo  is  the  free-electron 

mass,  but  the  value  is  not  well  characterized.  A  sizable  number  of  studies  have 
attempted  to  determine  the  ranges  of  effective  mass  both  for  electrons  and  holes,  and  are 
summarized  for  both  theoretical  and  experimental  studies  in  Table  2  below  [17]. 


Table  2  Electron  and  Hole  Effective  masses  in  GaN  (derived  from  [17]) 


Carrier  type 

Source  (exp/calc) 

Effective  Mass  (xm^) 

* 

Experimental 

0.20  -  0.27 

* 

m. 

Calculated 

0.18-0.22 

mt 

Experimental 

0.20  -  0.23 

Calculated 

0.18-0.23 

ml 

Experimental 

0.20 

ml 

Calculated 

0.17-0.20 

Kh 

Experimental 

1.0 -2.2 

Experimental 

0.75-0.9 

Experimental 

1.3 

m^  (DOS  effective  mass) 

Calculated 

1.5 

9 


GaN  is  typically  strongly  n-type,  with  the  source  of  the  shallow  native  donor 
being  the  subject  of  much  contention.  Early  efforts  seemed  to  point  to  a  nitrogen  vacancy 
(Vn)  [18]as  the  source,  which  was  contested  by  initial  theoretical  results.  The  nature  of 
the  native  donor  in  GaN  is  not  firmly  established  at  present,  with  arguments  having  been 
put  forth  for  both  a  defect  donor  and  for  impurity  doping.  As  recently  as  June  of  2000, 
Van  Nostrand  [19]  laid  out  many  theoretical  and  experimental  studies  which  illustrate  the 
controversy;  after  first  listing  the  sources  supporting  the  “overwhelming  evidence  [that] 
the  residual  donor  concentration  in  native  GaN  is  due  to  N  vacancies”,  the  “substantial 
body  of  works  to  the  contrary”  is  presented.  While  it  is  still  not  conclusively  resolved,  the 
nitrogen  vacancy  (rather  than  impurities  such  as  Si  or  O)  is  the  leading  candidate  for  the 
n-type  auto-doping  (or  defect  doping)  in  GaN,  which  can  reach  carrier  concentrations  of 
10  cm'  at300K[10].  High  quality  GaN  crystals  are  now  available  with  room 
temperature  electron  concentrations  of  5- 10^^  cm'^  [10]. 

Within  the  last  year  experimental  evidence  has  pointed  to  impurity  doping  as  the 
source  of  the  n-type  autodoping  in  GaN  [20] ;  the  impurity  doping  hypothesis  is  therefore 
gaining  acceptance.  There  is  still  some  controversy  associated  with  this  view. 

Figure  4  shows  the  temperature  dependence  of  the  intrinsic  carrier  concentration 
in  GaN  [14].  Fitting  the  temperature  dependence,  as  taken  from  the  data  of  Figure  4, 
yields  the  relationship 

n.  =110^^e  ^  ^  ^  cm  ^  (2) 
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Figure  4.  Temperature  dependence  of  the  intrinsic  carrier  concentration  in  GaN  [14]. 


The  n-type  nature  of  undoped  GaN  has  contributed  to  the  well-noted  difficulty  in 
p-type  doping  of  GaN.  The  large  n-type  conductivity  tends  to  compensate  acceptor 
dopants  and  results  in  highly  resistive  materials.  Historically,  the  only  acceptor  that  has 
been  successfully  used  is  Mg  [21],  when  used  in  conjunction  with  electron  beam 
irradiation  or  thermal  annealing  to  convert  the  resistive  compensated  material  to 
conductive  p-type.  The  thermal  annealing  process  has  been  observed  to  remove 
hydrogen,  the  primary  compensating  agent.  These  processes  were  observed  to  be  of  little 
utility  in  reducing  the  compensation  of  other  group  II  elements,  and  to  date  no  reliable 
process  has  been  found  to  effectively  acceptor  dope  GaN  with  any  element  besides  Mg 
[10]. 
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Even  though  Mg  is  the  most  prominent  p  type  dopant  in  GaN,  the  acceptor  level 
of  Mg  in  GaN  has  been  shown  to  be  hundreds  of  meV  above  the  valence  band  level  - 
approximately  220  meV  for  optical  determinations  of  ionization  energy  and  varying  with 
acceptor  concentration  from  125  meV  to  180  meV  for  thermal  ionization  energies.  These 
values,  when  extrapolated  to  a  zero  acceptor  density,  are  consistent  with  the  optical  value 
of  220  meV  to  within  the  experimental  uncertainty  of  ±20  meV  [22] . 

The  large  band-gap  (3.49  eV)  of  GaN,  which  leads  to  a  thermally  insensitive 
semiconductor  material,  also  leads  to  a  large  breakdown  field.  This  large  breakdown 
field  strength,  calculated  to  be  greater  than  3  MV/cm  [23]  and  indicated  by  some  sources 
to  be  greater  than  5  MV/cm,  makes  GaN  a  good  candidate  for  high  power  operation. 
These  values  are  between  five  and  eight  times  the  values  of  Si  and  GaAs  [12].  Measured 
values  of  the  breakdown  field  and  voltage  for  GaN  are  shown  in  Figure  5  as  a  function  of 
sample  temperature. 

GaN  also  has  large  carrier  velocities,  which  enable  high  speed  and  high  frequency 
operation  for  switching  and  RF  devices.  GaN  has  exhibited  peak  electron  velocities  of 
3.010’  cm! sec  and  saturation  velocities  of  1.510’  cm! sec  [12]. 
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Figure  5.  Normalized  breakdown  field  (Fi)  and  voltage  (Vi)  in  GaN  samples. 
Normalization  is  to  1-2- 10^  V/cm  and  42  V  respectively.  Samples  are  p'^-p-n  diodes.  [14] 


GaN,  like  all  of  the  Ill-Nitrides,  is  highly  piezoelectric  and  is  highly  polarized,  by 
semiconductor  standards.  The  piezoelectric  and  spontaneous  (zero  field)  polarization  of 
GaN  is  ten  times  that  of  conventional  compound  semiconductors.  The  field  due  to 
spontaneous  polarization  of  Ill-Nitrides  can  reach  3.0  MV/cm  and  for  AlGaN/GaN 
heterostructures  the  piezoelectric  field  may  reach  2.0  MV/cm  [24]. 

One  of  the  biggest  challenges  in  the  use  of  GaN  for  electronic  devices  has  been 
the  development  of  suitable  substrates  for  GaN  device  growth.  GaN  has  proven  very 
difficult  to  produce  in  a  large  crystal  formats,  primarily  due  to  a  low  solubility  of  N2  in 
Ga,  and  is  still  not  commercially  available.  SiC  and  AI2O3  (sapphire)  are  the  substrate 
materials  of  choice  today  for  GaN  devices,  although  each  has  associated  difficulties.  AIN 
substrates  may  soon  solve  some  of  these  difficulties  with  lattice  mismatch  and  thermal 
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conductivity  issues  and  become  commercially  available.  Fifteen  millimeter  (mm) 
substrates  are  commercially  available,  and  substrates  up  to  50  mm  are  projected  to  be 
available  soon  [25].  Thick  (over  200  |im)  free-standing  GaN  films  with  low  defect 
densities  have  been  produced  by  HVPE  growth  and  subsequent  removal  of  sapphire 
substrate  material.  These  films  have  low  defect  densities  due  to  the  removal  of  the 
substrate  interface  region  and  hold  promise  for  production  of  large  GaN  wafers  [26] 

Sapphire  substrate  materials  offer  the  advantages  of  widespread  availability,  low 
cost,  hexagonal  crystal  symmetry,  and  ease  of  handling.  Due  to  the  low  solubility  of  N2 
in  Ga  and  high  vapor  pressure  of  N2  over  Ga,  extremely  high  temperatures  and  pressures 
(in  excess  of  1200°  and  12  kbar)  are  required  for  HVPE  and  MOVPE  GaN  growth  - 
making  the  high  temperature  stability  of  sapphire  an  important  factor  [10]. 

Sapphire  substrate  issues  are  primarily  centered  around  the  lattice  mismatch  (~ 
16%  on  the  (0001)  plane)  [10],  which  induces  large  dislocation  defect  densities  near  the 
interface,  and  the  extremely  low  thermal  conductivity  of  sapphire  (0.42  W/cm-°K),  which 
is  untenable  for  the  high-power  devices  that  GaN  is  suited  for.  It  is  possible  to 
compensate  for  the  interface  defect  density  by  growing  a  thin  GaN  buffer  layer  before  the 
device  layers  are  grown.  This  process  tends  to  isolate  the  defects  from  the  active  device 
region,  but  adds  an  additional  processing  step  to  the  growth  process,  increasing  time  and 
cost  to  produce  devices.  The  difficulty  with  the  low  thermal  conductivity  of  sapphire  is 
not  easily  solved;  as  a  result,  most  high  power  devices  where  thermal  energy  dissipation 
is  an  issue  are  grown  on  SiC  substrates  to  take  advantage  of  its  high  thermal  conductivity. 

Eactors  that  SiC  has  in  its  favor  as  a  GaN  substrate  include  a  smaller  lattice 
mismatch  than  sapphire  (>  3.5%)  and  good  binding  qualities  to  GaN  epilayers  [27],  much 
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better  thermal  conductivity,  and  easier  formation  of  facets  than  in  sapphire.  However; 

SiC  is  much  more  expensive  to  produce  than  sapphire  substrates  and  there  are  issues 
associated  with  differing  crystal  symmetries  (which  produce  structural  defects  in  the  GaN 
epilayer),  as  well  as  surface  preparation  difficulties  that  must  be  considered  [10].  SiC  is 
the  substrate  of  choice  for  high-power  applications  due  primarily  to  its  ability  to  conduct 
thermal  energy  away  from  the  device. 

GaN  Growth  Methods 

A  variety  of  growth  methods  and  techniques  are  used  to  produce  epitaxial  GaN 
layers  on  the  substrates  mentioned  previously.  Short  descriptions  of  the  methods  used  to 
produce  the  materials  in  this  dissertation  are  provided  below. 

Molecular  Beam  Epitaxy  (MBE) 

Much  of  the  material  studied  in  this  dissertation  was  produced  via  MBE,  a  growth 
method  used  to  produce  thin  layers  of  GaN  on  a  prepared  substrate  material  via  atomic 
deposition  in  ultra-high  vacuum  (UHV)  conditions.  The  primary  limiting  factor  inherent 
in  this  technology  is  the  relatively  low  growth  rate.  Growth  rates  are  limited  indirectly  by 
the  requirement  to  produce  atomic  N;  cracking  of  the  N2  precursor  into  atomic  nitrogen  is 
a  high  energy  process  requiring  formation  of  a  plasma  to  extract  atomic  N.  The 
requirement  to  dissociate  the  N2  precursor  limits  the  use  of  conventional  effusion  cell 
MBE  systems  and  drives  the  GaN  production  process  to  an  RE  or  Electron  Cyclotron 
Resonance  (ECR)  microwave  plasma  source  for  the  N  cell  [28].  The  resulting  trade-off 
is  between  the  growth  rate,  which  demands  higher  plasma  excitation  powers,  and  ion- 
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beam  damage,  which  is  limited  by  lower  excitation  power.  If  higher  temperatures  are 
used  (700°  C  and  above)  an  ammonia  source  can  be  used  as  the  N2  precursor,  limiting  ion 
beam  damage  [29] ;  however,  this  introduces  defects  due  to  thermal  cycling. 

Primary  advantages  of  MBE  epitaxial  growth  are  the  relatively  low  growth 
temperature  and  the  low  impurity  content  (particularly  with  respect  to  H  impurities)  due 
to  the  UHV  deposition  [30] .  Low  growth  temperatures  lead  to  lower  thermal  stresses  and 
thus  to  lower  levels  of  gross  defects  than  occur  with  high  temperature  growth  processes. 
The  UHV  deposition  process  limits  the  amount  of  impurities  in  the  material,  particularly 
with  respect  to  oxygen  and  hydrogen  at  the  cost  of  introducing  lattice  defects  due  to  ion 
beam  damage. 

Hydride  Vapor  Phase  Epitaxy  (HVPE) 

HVPE  is  was  one  of  the  first  growth  methods  used  to  prepare  GaN  epitaxial  layers 
and  is  today  a  leading  contender  for  the  growth  of  large  diameter,  thick,  high  quality  GaN 
layers  [30].  HVPE  is  conducted  in  a  hydrogen  environment  at  temperatures  near  lOOOC, 
where  the  precursor  compounds,  typically  NH3  and  GaCl  (produced  by  flowing  HCl  over 
metallic  gallium),  react  via  the  following  equation 

GaCl  +  NH^  GaN  +  HCl  +  (3) 

and  are  deposited  on  the  substrate  material.  The  hydrogen-rich  environment  leads  to 
incorporation  of  H  impurities  into  the  material. 

HVPE  material  is  characterized  by  low  defect  densities  and  correspondingly  good 
material  properties.  GaN  material  grown  by  HVPE  and  dissociated  from  the  underlying 
substrate  is  currently  the  world’s  highest  mobility  GaN  [20].  Dissociation  from  the 
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underlying  substrate  layer  insures  that  conductive  interface  layers  are  not  formed,  and  the 
high-quality,  relatively  defect  free  bulk  material  is  retained.  Free-standing  GaN  layers 
over  200  |im  in  thickness  can  be  produced  by  this  means,  with  carrier  concentrations  in 
the  mid  to  high  10^^  range,  and  have  been  used  in  this  study. 

Metal-Organic  Vapor  Phase  Epitaxy  (MOVPE) 

MOVPE  is  a  method  of  epitaxial  growth  using  metal  alkyls  (such  as 
trimethylgallium,  Ga(CH3)3 )  as  group  III  precursors  and  ammonia  as  a  nitrogen  precursor 
for  a  high-temperature,  vapor  phase  deposition  of  GaN.  MOVPE  was  used  extensively 
early  in  the  production  of  the  Ill-nitrides  and  is  still  the  prevalent  method  of  commercial 
production  [31]. 

MOVPE  GaN  deposition  using  trimethylgallium  and  ammonia  as  precursors 
proceeds  by  a  complicated  set  of  reactions,  based  on  the  intermediary  formation  of  acid- 
base  pairs  [31].  The  complex  chemical  nature  of  this  growth  process,  which  is  still  not 
completely  understood,  is  a  source  of  defect  formation  in  MOVPE  GaN  layers  since  the 
adduct  pairs  created  may  be  deposited  in  the  material  intact.  Other  defect  sources 
inherent  in  this  growth  method  include  the  formation  of  nitrogen  vacancies  due  to  the 
high  temperatures  required  to  dissociate  the  ammonia  precursor,  and  the  inclusion  of 
carbon  and  methyl  molecules  into  the  lattice  [32]. 

MOVPE  grown  materials,  while  not  used  in  this  dissertation,  were  analyzed  by 
Carlos  using  EPR  [33].  Measurements  in  MOVPE  GaN  were  shown  to  be  very  similar  to 
those  made  in  HVPE  GaN. 
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Radiation  interactions 


The  expected  effects  of  ionizing  radiation  on  solid  state  materials  such  as  GaN  are 
two-fold:  dose  rate  dependent  ionization  effects  and  total  dose  dependent  persistent 
effects.  In  the  case  of  energetic  electron  bombardment  both  effects  are  produced; 
however,  we  are  concerned  primarily  in  this  study  with  the  persistent  effects  of  the 
radiation  interaction  and  will  not  consider  dose-rate  dependent  transient  effects  (mainly 
photo-ionization) . 

The  primary  mechanism  for  introduction  of  persistent  radiation  effects  is  the 
displacement  of  lattice  atoms  (Ga  or  N)  through  elastic  collisions  with  energetic  particles. 
Binding  energies  of  the  lattice  constituents  vary  with  the  particular  sublattice  considered 
and  are  anisotropic  in  nature  [34] .  The  rates  of  displacement  damage  for  Ga  and  N 
sublattices  are  a  function  of  both  atomic  binding  energy  and  collisional  energy  transfer. 

In  general,  Ga  atoms  are  bound  less  tightly  than  the  N  atoms  in  GaN;  however,  the 
energy  deposited  through  an  elastic  collisional  process  is  heavily  dependent  on  the  lattice 
atom’s  mass  and  so  much  more  energy  can  be  imparted  to  the  less  massive  nitrogen 
atoms.  Over  the  ranges  of  energy  considered  in  this  study  (0.5- 1.5  MeV)  both  types  of 
damage  are  expected  to  occur  simultaneously. 

Displacement  energies  for  the  lattice  constituents  in  GaN  have  been  calculated  by 
Nord  [35].  The  minimum  displacement  energies  are  22±leV  for  gallium  and  25  ±  leV 
for  nitrogen,  where  the  minimum  is  taken  over  all  angles.  For  an  average  over  all 
possible  angles,  the  displacement  energies  are  45  ±  leV  for  gallium  and  109  ±  2eV  for 
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nitrogen.  The  maximum  energy  that  may  be  imparted  to  a  lattice  atom  by  an  energetic 
electron  via  coulomb  scattering  is  given  by 


{e  +  2m  _c^) 

^max  _  2,-—-^ _ ^—-E 


(4) 


so  that  the  maximum  energy  available  to  a  lattice  atom  via  an  interaction  with  a  1  MeV 
electron  is  62  eV  for  gallium  atoms  and  309  eV  for  nitrogen  atoms.  The  maximum 
energy  transfer  as  a  function  of  the  incident  particle  energy  is  shown  in  Figure  6  for  both 
gallium  and  nitrogen.  Calculations  of  equation  4  shown  in  Figure  6  indicate  that  not 
only  can  one  expect  point  defects  on  the  nitrogen  sublattice  at  incident  electron  energies 
above  0.3  MeV,  but  that  gallium  sublattice  point  defects  should  appear  at  incident 
electron  energies  around  0.5  MeV.  Thus,  point  defects  on  both  sublattices  are  expected 
for  electron  irradiation  at  1.0  MeV.  The  possibility  of  knock-on  damage  exists,  although 
large  damage  cascades  are  not  expected  from  the  knock-on  particles,  which  are  limited  to 
about  290  eV  of  kinetic  energy  for  nitrogen  atoms  and  41  eV  for  gallium  atoms.  This 
follows  from  Nord’s  calculations  of  secondary  defect  production  which  show  that  for 
recoil  energies  of  200  eV,  less  than  one  additional  point  defect  in  either  sublattice  (per 
recoil  atom)  is  expected  from  either  species  of  recoil  atom  [35]. 
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Figure  6.  Maximum  energy  transfer  for  gallium  and  nitrogen  as  a  function  of  the  incident 
electron  energy.  The  average  and  minimum  displacement  energies  for  both  gallium  and 
nitrogen  are  shown  as  horizontal  lines.  Values  are  calculated  via  Equation  4. 


This  evaluation  is  reinforced  by  the  reported  observation  of  dislocation  damage 
in  GaN  following  electron  irradiation  at  approximately  300-350  keV  irradiation  energies 
[36]. 

The  types  of  radiation-induced  point  defects  produced  will  primarily  be  nitrogen 
vacancies  (V^),  nitrogen  intersitials  (/^),  and  the  corresponding  gallium  sublattice  point 

defects  Vca  and  Ica-  These  point  defects  are  expected  to  be  mobile  in  the  lattice  for  all 
but  the  lowest  temperatures  [37],  leading  to  the  potential  formation  of  a  variety  of 
complexes  with  impurities  in  the  lattice  or  with  other  defects. 
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Introduction  rates  of  electron  irradiation  induced  defects  in  GaN  are  not  well 


characterized  due  to  the  inherent  difficulty  in  separating  the  effects  of  different  types  of 
sublattice  damage  and  the  compensating  effects  of  the  various  defects.  While  the  exact 
classification  of  defect  types  in  GaN  is  uncertain,  many  defects  (both  shallow  and  deep 
states)  have  been  experimentally  observed  by  Hall  effect  measurements,  DLTS,  PL,  and 
various  magnetic  resonance  techniques,  including  radiation  induced  states.  The 
activation  energies  and  identity  of  these  states  are  catalogued  graphically  in  Figure  7. 
While  many  of  the  states  depicted  are  identified  with  a  particular  impurity  or  defect  type, 
one  should  bear  in  mind  the  tentative  nature  of  many  such  identifications  and  the 
controversy  which  surrounds  some  of  these  assignments. 
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Electron  Paramagnetic  Resonance  (EPR)  Spectroscopy 

Basic  elements  of  EPR  spectroscopy  will  be  discussed  in  this  section,  while 
elements  of  EPR  theory  appropriate  to  the  analysis  of  measurements  in  this  dissertation 
will  be  discussed  in  the  EPR  Theoretical  Considerations  section.  Elements  of  EPR  theory 
introduced  in  this  section  will  present  a  meaningful  introduction  to  the  physics  of  the 
technique  and  provide  a  background  for  the  theoretical  methods  presented  later. 

EPR  spectroscopy  is  a  powerful  experimental  technique  allowing  direct 
observation  of  unpaired  electron  spins,  including  unpaired  spins  associated  with  a  defect 
or  impurity  site  in  a  solid  state  material  via  microwave  absorption.  In  some  instances, 
identification  of  the  host  site  or  atom  is  available  via  interpretation  spectral  features 
created  by  the  spin  coupling  with  the  host  nucleus;  however,  this  is  not  guaranteed. 
Electron  Spin  Resonance  (EPR)  and  Electron  Paramagnetic  Resonance  (EPR)  are 
generally  accepted  as  synonymous  terms;  EPR  will  be  used  in  the  balance  of  this  work 
since  it  succinctly  captures  the  essence  of  the  technique. 

Since  the  recording  of  the  first  EPR  spectrum,  by  Zavoisky  in  1945,  studies  of 
impurities  in  solid-state  crystal  structures  have  been  an  important  part  of  this 
experimental  field.  The  early  development  of  EPR  techniques  was  performed  primarily 
at  Clarendon  Eaboratory  in  Oxford,  where  the  primary  application  was  the  study  of  metal 
ions  in  crystalline  lattices  [39].  Abragam  was  similarly  engaged  in  crystalline  studies 
when  he  pioneered  the  analysis  of  hyperfine  spectra,  particularly  with  regard  to  analysis 
of  large  hyperfine  coupling  constant  materials  [40].  Indeed,  by  1967  Bleaney  would 
write  that. 
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“Electron  Paramagnetic  Resonance  has  developed 
primarily  into  a  tool  for  investigation  of  effects  in  the 
solid  state,  in  which  additional  information  is  provided 
by  hyperfine  structure;  for  example  the  latter  can 
identify  the  nucleus  (or  nuclei)  with  which  magnetic 
electrons  interact,  and  give  the  strength  of  the 
interaction”  [41] 


Since  Bleaney’s  writing,  EPR  has  become  an  important  tool  in  molecular 
chemistry  where  it  is  used  as  a  sensitive  probe  of  the  bond  structures  and  bonding 
electron  states;  it  remains  a  powerful  tool  for  investigation  of  the  solid  state. 

Today,  EPR  spectroscopy  remains  a  valued  tool  for  defect  studies  in  crystalline 
materials,  along  with  methods  that  have  been  derived  from  traditional  EPR.  Derivative 
methods  of  EPR  include:  Optically  Detected  Magnetic  Resonance  (ODMR),  a  method 
that  obtains  increased  sensitivity  but  must  depend  upon  photonic  transistions  between 
states  [42];  Electronic  Nuclear  Double  Resonance  (ENDOR),  a  method  that  stimulates 
both  electronic  and  nuclear  resonances  [43];  and  Electrically  Detected  Magnetic 
Resonance  (EDMR),  a  method  using  electrical  measurements  to  enhance  detection 
sensitivity  [44].  Many  of  these  methods  trade  increased  signal  sensitivity  for  the 
inclusion  of  secondary  processes,  such  as  photonic  transitions  or  nuclear  magnetic 
resonances,  that  can  convolute  the  analysis  of  the  resulting  data.  Many  good  references 
are  available  on  EPR  spectroscopy  and  related  methods,  including  Atherton  [39],  Weil 
[45]  and  Poole  [46]. 

At  its  most  basic  level,  EPR  spectrometry  is  based  on  measuring  the  transition 
energies  and  probabilities  of  photon  induced  transitions  between  Zeeman  levels  of  a 
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charged  particle  in  a  magnetic  field  (B).  Zeeman  splitting  for  an  electron  in  a  magnetic 
field  of  magnitude  |5|  is  described  by  the  expression 

AE  =  ±ju^  |b|  (5) 

where  the  Bohr  magneton  is  given  by 

Bb=- —  (6) 

2m^ 

and  =  2.00233  is  the  measured  gyromagnetic  ratio  of  the  free  electron.  The  g-value 
is  a  scaling  factor  that  describes  the  magnitude  of  the  energy  splitting  of  a  charged 
particle  in  a  magnetic  field.  A  naive  classical  electromagnetic  model  of  a  free  electron  as 
a  spinning  charged  sphere  interacting  with  an  external  field  predicts  energy  levels  of 

AE  =  ±ju,m  (7) 

Therefore,  when  experiments  determined  the  scaling  factor  between  AE  and  |B|  to  be 

2.0023  ■  ju^ ,  the  value  2jUg  was  referred  to  as  the  “anomalous  electron  magnetic 
moment”  of  the  electron  [39],  a  term  which  has  fallen  out  of  use.  The  g-value  can  thus  be 
thought  of  as  the  scale  of  the  departure  from  the  classically  expected  magnetic  moment 
values  of  a  magnetic  particle. 

The  term  “electron  paramagnetic  resonance”  arises  from  the  requirement  for  the 
electron  to  be  in  an  unpaired,  or  paramagnetic,  spin  state  in  order  to  carry  out 
spectroscopy.  For  electrons  in  diamagnetic  (paired)  states,  all  quantum  numbers  are 
identical  with  the  exception  of  the  spin  quantum  number  ,  which  takes  on  values  of 
-1-1/2  and  -1/2  for  the  paired  electrons.  In  the  diamagnetic  state,  photon  absorption  is 
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prohibited  by  the  exclusion  principle;  thought  of  in  another  way,  absorption  and  emission 
at  the  excitation  frequency  would  be  exactly  balanced,  making  spectroscopy  impossible. 

In  a  real  system  such  as  transition  metals  in  crystalline  structures,  organic 
radicles,  or  donor  bound  electrons,  the  g-value  can  and  does  vary  from  the  free-electron 
g-value.  These  variations  in  the  g-value  are  attributed  to  the  admixture  of  the  orbital 
angular  momentum  with  the  spin  angular  momentum  [47]  (even  though  the  application  of 
the  external  field  quenches  the  orbital  angular  momentum  of  the  ground  state  due  to 
arguments  stemming  from  the  non-degenerate  nature  of  the  states  [48]).  Non-spherical 
symmetry  is  often  associated  with  this  admixture,  producing  anisotropy  in  the  g-value, 
which  is  therefore  properly  termed  the  g-tensor,  with  respect  to  the  orientation  of  the 
applied  magnetic  field.  In  crystalline  systems,  g  value  anisotropy  is  typically  present, 
with  the  relationship  between  angle  6  and  g-value  given  by  [33] 

g(0)  =  (g^sin^(0)+g^cos^(0))  (8) 

for  the  case  of  two  dimensional  anisotropy  (for  a  more  general  treatment  in  three 

dimensions,  see  [46]).  Here  the  parallel  and  perpendicular  subscripts  indicate  the  angle 
of  the  semiconductor  c  axis  with  the  applied  magnetic  field. 

Calculation  of  the  g-value  from  the  experimental  data  is  performed  by 
comparing  the  measured  Zeeman  transition  energy  (derived  from  the  microwave 

frequency  via  AE  =  hv and  the  measured  magnetic  field  strength  |.b|  at  the  microwave 
absorption  peak.  Given  the  experimentally  determined  values  of  and  Bres  (the 
resonant  magnetic  field  strength),  the  gyromagnetic  ratio  is  calculated  by 
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g  = 


hv 


JilW 


jd  B 


(9) 


This  g-value  is  unique  to  the  species  and  site  depending  upon  binding  status  and 
surroundings.  According  to  equation  9,  the  g-value  of  a  site  may  also  be  expressed  in 
terms  of  the  resonant  magnetic  field  value  and  this  is  used  extensively  in  laboratory 
practice.  In  certain  instances,  coupling  of  the  electron  and  nuclear  spins  associated  with 
the  site  lead  to  a  splitting  of  the  Zeeman  energy  levels  into  “hyperfine”  levels.  The 
magnitude  of  this  splitting  is  dependent  upon  the  magnetic  properties  of  the  nucleus, 
allowing  the  identification  of  the  nuclear  species  involved  in  the  coupling.  When 
instrument  sensitivity  and  resolution  are  sufficient,  additional  splittings  due  to  coupling 
with  surrounding  (nearest  neighbors)  nuclei  may  be  resolved  as  well. 

As  described  above,  coupling  of  the  electron  spin  (s  )  with  a  nuclear  spin  or  spins 
(/  )  leads  to  a  splitting  of  the  Zeeman  levels  into  (2/  -i- 1)  levels.  To  first  order  the 
splitting  is  described  by  the  hyperfine  coupling  constant  (A)  and  the  nuclear  spin  of  the 
interacting  nuclei,  with  energy  levels  given  by 

AE  = //gg|5| -I- (10) 
where  the  nij  =  (- 1 -I  + 1,.../),  producing  the  (2/  -I- 1)  equally  spaced  equal  magnitude 
lines  in  the  resulting  spectrum.  For  traditional  EPR  measurements,  where  nuclear  spin 
modes  are  not  excited,  the  allowed  transitions  measured  are  determined  by  the  selection 
rules  Am^j  =  ±1,  Am^  =  0 .  This  simple  view  of  the  effects  of  hyperfine  coupling  is 
insufficient  to  explain  all  spectral  features  and  must  be  modified;  however,  it  serves  to 
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communicate  the  basic  nature  of  the  hyperfine  interaction.  Further  developments  will  be 
left  until  the  section  on  EPR  Theoretical  Considerations  in  the  theory  section,  in  which 
derivations  appropriate  to  the  cases  in  this  dissertation  will  be  provided. 

Since  the  sites  under  study  are  not  truly  isolated  spin  states  or  systems  of  states, 
interactions  with  the  surroundings  (termed  the  lattice)  are  important  in  EPR  spectroscopy, 
and  enter  into  our  consideration  primarily  through  the  process  of  spin-lattice  relaxation. 
The  process  of  microwave  absorption  disturbs  the  thermodynamic  equilibrium 
concentrations  of  the  electrons  under  study  and  relaxation  to  the  equilibrium  levels 
proceeds  through  the  exchange  of  energy  with  the  surrounding  lattice.  This  relaxation  is 
characterized  by  a  lattice  relaxation  time  constant  (t ) .  This  relaxation  time  is  related  to 
the  linewidth  (or  the  uncertainty  in  the  energy  of  the  transition)  by  the  uncertainty 
relationship 

AE At>h  orTAffl>l  (11) 

which  asserts  that  as  the  lifetime  (r)  becomes  short,  to  the  point  of  equaling  the  inverse 
of  the  microwave  frequency,  linewidths  can  become  so  large  as  to  make  the  signal 
unmeasurable  [49] .  Eor  this  reason,  EPR  measurements  on  solid  state  materials  are 
typically  made  at  cryogenic  temperatures  (<  20 K)  to  limit  the  phonon-mediated 
interactions  of  the  spin  states  with  the  crystalline  lattice  [50].  A  benefit  of  conducting 
measurements  at  these  temperatures  is  the  enhanced  population  difference  between  the 
upper  and  lower  states.  The  EPR  signal  is  proportional  to  the  population  imbalance,  and 
is  given  by 
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where  Yo  is  the  signal  as  the  temperature  approaches  0  °K.  For  X  band  (10  GHz)  cavities, 
and  g  values  near  2.0,  the  population  of  the  lower  levels  is  greatly  enhanced  at  cryogenic 
temperatures.  In  this  range  the  sensitivity  is  not  a  strong  function  of  temperature,  as 
h(o»kT  and  the  signal  strength  approaches  the  constant  value  Yq. 

EPR  spectroscopy  systems  are  tasked  with  measuring  a  relatively  low  intensity 
signal  against  a  noisy  background,  prompting  the  use  of  magnetic  field  modulation  and  a 
lock-in  amplifier  system  to  increase  signal  to  noise  ratios.  The  magnetic  field  strength  is 
varied  in  X-band  EPR  systems  and  the  microwave  frequency  held  constant  due  to 
difficulty  in  designing  variable  frequency  microwave  sources  and  cavities.  The  signal 
resulting  from  this  magnetic  field  modulation  is  the  derivative  of  the  absorption  profile 
and  must  be  integrated  to  obtain  the  absorption  behavior.  Given  this  magnetic  field 
sweeping  and  the  magnetic  field  modulation,  the  evolution  of  the  first  derivative 
absorption  signal  becomes  clear.  This  is  illustrated  in  Eigure  8,  borrowed  from  [46]. 

Note  that  in  Eigure  8,  the  independent  variable  is  the  magnetic  field  intensity  H  rather 
than  the  magnetic  flux  density  B  that  has  been  used  to  this  point;  the  two  values  are 
simply  related  by  the  magnetic  permeability  via  B  =  juH  . 


A 


Figure  8.  Relationship  between  magnetic  field  modulation  (vertical  waveform), 
absorption  spectrum  (Y(x)),  and  resultant  signal  (horizontal  waveform).  [46] 


Applicability  of  EPR  Spectroscopy  to  Radiation  Effects  Studies 

EPR  spectroscopy  is  distinguished  from  many  of  the  derivative  methods  (such  as 
ODMR,  EDMR,  etc.)  by  the  direct  observation  of  paramagnetic  states  without  the 
requirement  for  corresponding  transitions  to/from  other  discrete  energy  levels  in  the 
material.  This  fact  simplifies  the  analysis  of  EPR  data  with  respect  to  relative  changes  in 
the  population  of  the  sites  under  investigation;  the  signal  observed  is  not  influenced  by 
changes  in  the  population  states  of  other  energy  levels,  an  important  distinction  to 
remember  when  comparing  EPR  and  Hall  effect  data.  The  direct  measurement  of  a 
particular  site’s  (such  as  the  unionized  shallow  donors)  population  makes  the  application 
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of  the  EPR  method  to  radiation  studies  attractive  because  of  the  inherent  ability  to 
determine  relative  changes  in  the  pre-irradiation  and  post- irradiation  site  populations. 

Although  ODMR  or  PL-ODMR  both  can  offer  enhanced  signal  sensitivity,  they 
suffer  from  the  inter-dependence  of  different  states  to  produce  the  observed  resonance. 
Resonance  signals  in  these  methods  must  meet  the  resonance  requirements  in  the  applied 
magnetic  field  and  must  meet  resonance  requirements  for  the  optical  transitions  between 
inter-bandgap  energy  transitions.  While  the  ability  to  selectively  excite  these  optical 
transitions  allows  greater  sensitivity,  at  least  two  states  (and  often  more)  are  involved  in 
the  transition  process,  with  all  of  the  participating  states  conceivably  exhibiting  spin 
dependent  excitation  properties.  Unfortunately,  the  population  statistics  of  any  of  these 
states  may  prove  a  transition  rate-limiting  factor,  making  determination  of  the  actual 
population  of  the  state  of  interest  impossible. 

EPR  Measurements  in  GaN 

GaN  has  been  studied  via  EPR  and  ODMR  spectroscopy  since  the  early  1990’s, 
due  primarily  to  interest  in  the  nature  of  the  native  defects  which  plagued  early  growth 
efforts.  Eor  effective  mass  (hydrogenic)  donors  in  GaN  the  g-value  is  typically  near  1.95 
(1.93  <  g  <  1.96)  and  is  anisotropic.  Gyromagnetic  ratios  of  effective  mass  donors  are 
slightly  anisotropic  with  respect  to  magnetic  field  direction,  with  values  of 
g_i=  1.9485±0.0002  and  gy=  1.95 15 ±0.0002  as  determined  by  several  studies  [51]. 
Several  other  signals,  characteristic  of  various  defects  or  impurities,  have  also  been 
detected  via  EPR  and  ODMR,  and  are  shown  in  Table  3.  GaN  grown  on  sapphire 
substrates  also  displays  strong  signals  from  substrate  transition  metals,  and  these  are 
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differentiated  by  their  strong  angular  dependencies  and  saturation  at  low  microwave 


powers  [33].  It  is  important  to  note  that  effective  g-value  of  hydrogenic  donors  in  GaN  is 
dominated  by  the  conduction  band  structure;  therefore,  any  effective  mass  donor  should 
exhibit  behavior  similar  to  the  anisotropic  1.95  g-value  of  the  shallow  native  donor  in 
GaN  reported  by  Carlos  [33]. 


Table  3.  Observed  EPR  Signals  in  GaN 


^11 

g± 

Linewidth 
(AB  Gauss) 

Ascribed  Nature  /  Source 

1.9515  ±0.0002 

1.9485  ±0.0002 

22-170 

EM  Donor  [51] 

1.989  ±0.001 

1.992±  0.001 

130 

Deep  Donor  [51] 

2.004  ±0.001 

2.008  ±  0.001 

* 

Ga,-  [52] 

1.960  ±0.002 

-2.03 

* 

Unknown  [52] 

1.989 

1.992 

* 

Unknown  [52] 

2.08  ±0.01 

2.00  ±  0.01 

260 

Mg  [51] 

*  Reported  after  2.5  MeV  electron  irradiation,  via  ODMR. 


Hyperfine  splitting  in  GaN  has  been  observed  via  ODMR  in  GaN  [52].  Hyperfine 
splitting  in  GaN  is  expected  to  occur  predominantly  with  the  primary  lattice  constituents 
exhibiting  non-zero  nuclear  spin.  In  the  case  of  GaN,  both  of  the  primary  lattice 
constituents  have  non-zero  nuclear  spins,  and  both  have  two  naturally  occurring  isotopes. 
The  isotropic  hyperfine  coupling  constants,  nuclear  spin  values,  nuclear  g-values,  and 
isotopic  abundances  are  shown  in  Table  4  for  each  of  the  isotopic  species  occurring  in  the 
GaN  lattice. 
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Table  4.  EPR  Hyperfine  Parameters  for  Ga  and  N  [53] 


Nuclide 

Nuclear 

Spin 

Natural 

Abundance 

Isotropic  Hyperfine 
Splitting  Constant  (mT) 

Nuclear 

g-value 

1 

99.63% 

64.62 

0.4038 

1/2 

0.366% 

-90.65 

-.5664 

3/2 

60.1% 

435.68 

1.3444 

^'Ga 

3/2 

39.9% 

553.58 

Simple,  first  order  approximations  to  the  hyperfine  splitting  )  are 
insufficient  to  predict  the  resonance  positions  )  when  the  hyperfine  coupling 
constant  becomes  large  with  respect  to  the  Zeeman  splitting.  Wiel  gives  a  rule  of  thumb 
for  comparing  these  magnitudes  at  microwave  frequencies  near  10  GHz  as  [54] 

Ajg^^^>l0mT.  (13) 

where  /  g^ju^  is  the  hyperfine  coupling  constant  expressed  in  magnetic  field  units. 

Since  the  measured  values  of  the  isotropic  hyperfine  splitting  constants  for  both 
gallium  isotopes  are  in  the  450-500  mT  range  and  is  approximately  64  mT  for  the 
predominant  isotope  of  nitrogen,  it  is  immediately  obvious  that  a  more  sophisticated 
approach  to  computing  hyperfine  splittings  will  be  necessary.  This  approach  will  be 
outlined  in  the  Theoretical  Considerations  section. 

Hall  Effect  Measurements 

In  addition  to  EPR  measurements,  temperature  dependent  Hall  effect 
measurements  of  samples  were  performed  to  determine  carrier  density  and  mobility  of 
samples  before  and  after  irradiation.  Hall  measurements  are  used  to  determine  the  energy 
levels  and  concentrations  of  donor  and  acceptor  populations  before  and  after  irradiation. 
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The  most  powerful  use  of  Hall  effect  measurements  in  this  study  occurs  when  the  data  is 
linked  with  the  EPR  measurements  to  determine  the  relative  contributions  of  the  EPR 
active  centers  to  carrier  concentration.  The  history  and  theory  of  Hall  effect 
measurements  are  widely  available  [55]  and  will  not  be  repeated  here  except  as  is 
necessary  for  the  current  analysis. 

Carrier  concentration  measurements  in  this  study  are  made  more  difficult  by  the 
existence  of  a  degenerate,  conductive  channel  at  the  sapphire/GaN  interface,  a  result  of 
the  high  defect  density  found  at  this  interface  [56].  After  analyzing  the  impact  of  this 
conductive  channel  and  isolating  the  behavior  of  the  GaN  epilayers  as  outlined  in  the 
theory  section,  estimates  of  the  GaN  carrier  concentration  are  possible.  Donor  and 
acceptor  concentrations  in  a  sample  can  be  determined  by  fitting  a  donor  activation 
model  to  the  carrier  concentration  data.  The  resulting  donor  concentrations,  when 
compared  with  the  relative  concentration  changes  in  EPR  active  sites,  may  substantiate 
the  EPR  data  and  allow  an  estimate  of  the  relative  contribution  of  the  EPR  site  to  the 
carrier  concentrations. 

The  most  important  distinction  between  the  concentrations  measured  via  the  Hall 
effect  and  the  EPR  measured  concentrations  is  in  the  source  of  the  concentration  data: 
the  Hall  effect  measures  the  carrier  density  and  mobility  whereas  the  EPR  process 
measures  the  concentration  of  the  donor  or  acceptor  sites  themselves,  as  identified  by  a 
unique  g-value.  Donor  and  acceptor  concentrations  must  be  extracted  from  Hall  effect 
data  by  modeling,  while  they  may  be  directly  measured  by  EPR  if  paramagnetic.  Thus, 
the  combination  of  these  two  experimental  methods  allows  a  complete  analysis  of  both 
donor  and  acceptor  populations  and  the  conduction  band  electron  concentration  derived 
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from  these  sites.  When  these  data  are  combined,  there  is  much  less  ambiguity  regarding 
the  nature  of  the  changes  in  carrier  concentrations  and  the  changing  donor  site 
populations  may  be  correlated  with  the  changes  observed  in  the  conduction  band 
populations. 
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III.  Theoretical  Considerations 


Theoretical  discussion  in  this  section  are  limited  to  development  of  theoretical 
models  found  to  be  useful  for  analysis  of  the  experimental  data  collected  during  the 
research.  The  development  of  a  mathematical  model  to  explain  the  observed  changes  in 
shallow  donor  concentrations  following  irradiation  is  presented  first.  Development  of  the 
spin  Hamiltonian  and  the  hyperfine  coupling  analysis  of  a  spin  3/2  nucleus  and  a  single 
electron  is  then  presented.  Both  of  these  topics  prove  to  be  useful  in  later  analysis. 

EPR  Theoretical  Considerations 

Donor  Passivation  Modeis 

EPR  measurements  measure  absorption  due  to  spin  state  transitions  of  unpaired 
electrons.  Since  EPR  measurements  in  solid-state  materials  are  typically  made  at 
temperatures  below  lOK,  almost  all  EPR  active  impurity  donors,  such  as  silicon  or 
oxygen  in  GaN,  will  be  in  their  un-ionized  state,  and  the  EPR  signal  can  safely  be 
assumed  to  arise  from  the  excitation  of  the  loosely  bound  donor  electrons.  Ionization  of 
these  shallow  donors  thus  promotes  the  donor  bound,  unpaired  electrons  to  the 
conduction  band,  rendering  the  resulting  positively  charged  donor  site  diamagnetic  and 
un-measurable  by  EPR.  Any  process  which  isolates  these  donors  in  the  positively 
charged  state  at  low  temperatures  or  which  involves  the  unpaired  electron  in  a  bonding 
arrangement  will  render  the  site  insensitive  to  EPR  absorption. 

Particle  balance  considerations  dictate  that  the  number  of  ionized  donors  is  given 
by 
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(14) 


where  n‘^'  is  the  total  number  of  shallow  donors,  is  the  number  of  neutral  donor 
atoms  at  low  temperature,  and  is  the  number  of  ionized  donors.  The  doubly  ionized 
donor  state  is  not  considered  as  a  possibility  at  the  temperatures  of  interest.  Donor  sites 
in  the  material  can  be  compensated  by  acceptors  or  electron  traps,  occurring  during  the 
growth  process  or  as  a  result  of  irradiation  respectively;  while  typically  goes  to  zero 
at  low  temperatures,  the  inclusion  of  acceptor  sites  or  electron  traps  may  increase  their 
concentration  by  reducing  conduction  band  electron  populations  to  the  point  that 
electrons  are  not  available  to  fill  the  ionized  donor  sites.  This  effectively  results  in  the 
pinning  of  the  Fermi  level  near  the  donor  level  to  account  for  the  increased  population  of 
ionized  donors  at  a  particular  temperature  [57]. 

From  charge  balance  considerations  the  carrier,  donor  and  acceptor  concentrations 
may  be  expressed  as 

«  +  +  (15) 

k  j  j 

where  k  and  j  are  summed  over  the  acceptor  sites  with  f  <  f  ^  and  donor  sites  with 

£j^j  >£^  respectively,  n  represents  the  total  conduction  band  electron  concentration  and 

p  is  the  hole  concentration.  Again  the  and  represent  the  ionized  acceptor  and 

donor  concentrations  respectively.  As  the  temperature  drops  towards  OK,  the  expression 
above  reduces  to 
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(16) 


no-Nl 


0 


where  is  the  concentration  of  acceptor  sites  (summation  over  k  levels  in  equation  15) 
occupied  by  an  electron  and  assuming  a  single  donor  level  is  at  or  above  the  Fermi  level. 
At  low  temperatures,  we  may  assume  that  all  of  the  acceptors  are  in  the  “occupied”  state, 
at  least  to  the  point  where  there  is  a  higher  acceptor  concentration  than  there  are  donor 
electrons.  When  equations  (14)  and  (16)  are  combined  and  solved  for  the  un-ionized 
donor  concentration  the  result  is 


0  + 

n  =  n  —  ti  =  n 
Hd  '‘'D  '‘'D 


■N] 


(17) 


so  that  the  EPR  measured  quantity  is  the  uncompensated  donor  concentration,  as  is 
expected.  This  result  suggests  that  two  types  of  processes  can  change  the  measured  EPR 
signal:  processes  which  remove  or  change  the  substitutional  donor  sites  by  involving  the 
loosely  bound  donor  electron  in  a  donor-defect  complex  bond,  or  processes  which  change 
the  acceptor  concentration  and  thereby  change  the  Eermi  level  pinning  energy  with 
respect  to  the  shallow  donor  energy  level. 

In  the  case  of  complex  formation,  mobile  point  defects  which  exhibit  acceptor¬ 
like  properties  and  are  negatively  ionized  will  bond  with  the  positively  charged,  ionized 
donor  sites  due  to  the  coulombic  attraction  between  the  sites.  An  example  of  this  process 
would  be  a  bound  state  of  silicon  substitutional  donors  with  mobile  nitrogen  interstitials 
in  GaN.  The  resulting  donor-defect  complex  may  not  exhibit  paramagnetic  properties 
and  certainly  will  be  distinct  in  those  properties  from  the  donor  site  which  it  replaces. 

This  process  will  clearly  require  a  level  of  mobility  by  the  defect  sites,  and  would  not  be 
expected  to  occur  for  very  low  temperature  irradiations;  this  has  been  experimentally 
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verified  by  ODMR  measurements  on  electron  irradiated  GaN  [37]  by  observations  of 
reduced  complex  formation  with  decreasing  temperature.  The  marked  temperature 
dependence  observed  in  this  study  is  indicative  of  a  mobility  limited  process  such  as 
donor-defect  complex  formation. 

Donor-defect  complex  formation  will  reduce  donor  concentration  in  a  manner  that 
is  dependent  upon  both  radiation  fluence  (by  way  of  the  defect  formation  rates)  and  donor 
concentration,  as  given  by  Boudinov,  et  al  [58]  as: 

n  ,  ^ 

i -  (is; 

where  ct(T)  is  a  temperature  dependent  “rate  constant”  for  the  donor  complex  formation 
process,  no  is  the  number  of  ionized  donors,  n^e/is  the  number  of  point  defects  available 
to  react  with  the  donors,  and  ricompiex  is  the  number  of  complexes  present  in  the  material. 
While  Boudinov  did  not  cast  the  rate  constant  as  a  temperature  dependent  function,  the 
dependence  of  the  process  on  the  defect  mobility  clearly  requires  that  this  rate  constant  be 
temperature  dependent.  It  is  also  clear  from  the  above  equation  that  the  rate  of  donor- 
defect  complex  formation  is  dependent  upon  the  number  of  unbound  impurity  donor  sites 
)  and  that  the  process  must  slow  as  more  of  these  donors  are  converted  to 
the  complexed  state. 

Donor  compensation  by  acceptor  sites  is  observed  in  as-grown  materials; 
irradiation  is  also  expected  to  produce  additional  acceptor  states  in  the  material  since 
some  radiation-induced  point  defects,  such  as  the  N  interstitial,  are  expected  to  be 
acceptor-type  sites  [6].  To  first  order  there  is  no  reason  to  expect  that  the  process  of 
defect  formation  is  dependent  upon  the  doping  density  in  the  virgin  material,  thus  the 
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process  of  acceptor  or  trap  compensation  of  shallow  donors  is  a  function  only  of  the 
radiation  dose  (or  fluence).  The  distinction  between  the  acceptor  state  and  the  acceptor 
trap  state  is  subtle:  an  acceptor  level  lies  below  the  Fermi  level,  while  an  acceptor  trap  is 
a  state  above  the  Fermi  level  which  exhibits  a  tendency  to  “trap”  an  electron  in  an 
unfilled  electron  orbital  or  bond.  The  trap  state,  being  above  the  Fermi  level,  will  trap 
electrons  which  are  then  thermally  excited  to  the  conduction  band,  a  process  which  is 
characterized  by  the  state  lifetime. 

Formation  of  compensating  acceptor  trap  sites  is  thus  independent  of  the  donor 
concentration  and  is  of  a  completely  different  nature  than  the  complexation  process 
described  above.  This  process  may  be  described  by  the  following  relationship 

NT=ccJ  (19) 

where  is  the  concentration  of  radiation  induced  acceptor-type  defects,  is  the 
damage  constant  associated  with  this  defect  type,  and  (j)  is  the  radiation  fluence.  As 
shown  previously,  donor  sites  are  compensated  by  impurity  acceptors  in  the  virgin 
material  and,  to  first  order,irradiation  should  not  affect  this  native  acceptor  concentration. 
Therefore,  in  any  relative  comparison  the  native  impurity  acceptor  concentrations  may  be 
ignored  and  the  radiation  induced  acceptor  states  treated  as  the  predominant  contributor 
to  the  change  in  donor  charge  state.  Unlike  the  complexation  process  described  above, 
the  compensation  process  does  not  change  the  fundamental  identity  of  the  donor  site  nor 
its  bandgap  level,  the  acceptors  instead  acting  to  affect  the  Fermi  level  and  thus  the 
charge  states  of  the  donors. 
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The  rate  equation  governing  the  implantation  and  population  of  radiation-induced 
defects  and  donor-defect  complexes  is  given  by  Titov  and  Kucheyev  as  [59] 

n  \ 

1 - ^20) 

where  their  formulation  has  been  modified  slightly  to  allow  the  inclusion  of  a  “stable” 
radiation-induced  acceptor  concentration  n  .  In  equation  (20)  p  is  the  defect 
introduction  rate, /is  the  radiation  flux  and  other  variables  are  as  they  have  been  defined 
previously.  In  this  formulation,  the  strict  rate  balance  between  defect  production, 
annealing,  and  donor-defect  complex  formation  is  not  preserved  since  the  build-up  of  an 
uncomplexed  defect  concentration  is  allowed  due  to  the  introduction  of  the  n''^  term. 
This  expression  is  amenable  to  a  numerical  solution,  which  is  presented  later  in  this 
discussion. 

If  the  formulation  of  equation  (20)  for  the  radiation-induced  trap  density  and 
donor-defect  complex  formation  are  considered  in  the  context  of  charge  balance,  the 
following  equation  is  the  result 

n  )-  (iv;  +  (iV  ).  (21) 

and  the  difference  between  the  pre-irradiation  and  post-irradiation  EPR  measured  donor 
populations  must  be  given  by 

(22) 

so  that  the  change  in  EPR  signal  is  a  function  of  the  fluence,  which  is  implicit  in  both  of 
these  terms,  and  the  initial  donor  concentration,  which  is  implicit  in  the  term. 


pfAt-n^ 


)^D  ^def 
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When  the  effects  of  the  defect  interactions  listed  above  are  evaluated,  it  is  clear 


that  there  are  three  cases  of  interest: 

a.  Donor-defect  complex  formation  is  the  only  process  affecting  the  donor 
population  (“complex  only”  model); 

b.  Radiation  induced  acceptors  limit  the  donor  population  that  can  be  maintained 
in  the  neutral,  “EPR  visible”  charge  state  (“compensation  only”  model);  or 

c.  Both  complex  formation  and  trap  introduction  combine  to  reduce  the  “EPR 
visible”  donor  population  (“complex  -i-  compensation”  model). 

The  effect  of  these  competing  processes  can  be  determined  to  some  extent  by  the 
evaluation  of  equation  (20)  by  numerical  means  and  determining  the  effect  on  EPR  or 
Hall  measurements  by  use  of  equation  (21). 

Results  of  this  modeling  indicate  that  the  effects  of  the  processes  listed  above  on 
the  donor  population  can  be  identified  to  some  extent.  Eigure  9  supports  this  contention 
by  illustrating  the  difference  between  the  “donor  complex  only”  and  “complex  -i- 
compensation”  models  presented  above.  Erom  the  simulation  results  shown  in  Eigure  9  it 
is  clear  that  for  the  correct  radiation  doses  a  distinction  can  be  drawn  between  the  models, 
with  the  “complex  only”  model  displaying  a  dependence  upon  the  initial  donor 
concentration  that  the  “complex  -i-  compensation”  results  do  not  show.  The  model  for 
“compensation  only”  donor  passivation  predicts  a  constant  decrease  in  the  donor 
population  for  a  given  fluence,  up  to  the  point  of  saturation.  The  same  is  true  for  the 
“complex  -I-  compensation”  model,  as  all  of  the  defects  not  involved  in  complex 
formation  are  available  to  compensate  any  remaining  donors.  In  both  of  these  situations, 
the  effect  of  the  irradiation  is  to  simply  reduce  the  donor  concentration  by  a  constant 
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amount  equal  to  the  number  of  radiation-induced  donors.  The  linear  relationship  of  the 
pre-irradiation  data  is  thus  reproduced  in  the  post-irradiation  case,  with  the  slope  of  the 
linear  relationship  preserved  and  the  intercept  moved  from  zero  to  the  donor  saturation 
value  (ie,  the  initial  donor  concentration). 
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Figure  9.  Results  of  the  numerical  solution  of  equation  20  for  different  initial  donor 
concentrations.  For  each  of  the  pre-irradiation  donor  concentrations  the  virgin  EPR 
system  response  is  shown  along  with  the  post-irradiation  results  for  the  “complex  only” 
model  and  the  “complex  -i-  compensation”  model. 


In  contrast,  the  “complex  only”  model  predicts  that  the  post-irradiation  donor 
concentrations  vary  from  the  linear  pre-irradiation  values.  This  relationship  is  not  linear, 
and  the  limiting  value  in  the  high  donor  concentration  regime  is  the  total  defect 
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concentration,  as  expected  since  a  high  initial  donor  concentration  raises  the  probability 
of  all  radiation  induced  defects  being  bound  into  donor-defect  complexes,  especially 
when  the  total  number  of  defects  is  much  smaller  than  the  initial  donor  concentration. 

The  impact  of  the  a  parameter  appearing  in  equation  (20)  is  important  in  this 
analysis,  particularly  since  it  is  a  free  parameter  in  this  model,  not  having  been 
experimentally  determined.  As  a  approaches  zero,  the  complexation  process  is  “turned 
off’  and  does  not  occur  in  appreciable  concentrations  regardless  of  the  defect  and  donor 
concentrations.  Conversely,  as  a  goes  to  very  large  values  the  complexation  process 
looses  the  non-linear  dependence  upon  the  donor  concentration  since  almost  every  defect 
site  available  will  be  complexed  despite  dwindling  concentrations  of  uncomplexed  donor 
sites.  It  is  only  for  values  of  a  on  the  order  of  10’  to  10’  cm  /sec  that  the  nonlinear 
donor  density  dependence  is  obvious.  Therefore,  differences  in  complex  formation  and 
defect  compensation  effects  may  be  obvious  only  with  a  priori  knowledge  about  the 
magnitude  of  this  interaction  rate  constant;  however,  an  observed  deviation  from  linearity 
in  the  pre-  and  post- irradiation  donor  concentrations  would  positively  indicate  that 
complex  formation  is  the  predominant  process  arising  from  irradiation. 

These  results  imply  that  it  will  be  very  difficult  to  differentiate  between  these 
models  in  an  experimental  setting,  since  the  nitrogen  interstitial  is  known  to  be  an 
acceptor  site  and  the  presence  of  any  of  these  defects  outside  of  a  donor-defect  complex 
will  mask  the  effects  of  the  complex  formation  by  compensating  the  uncomplexed  donors 
to  the  level  of  the  total  defect  concentration. 
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Development  and  Application  of  the  Hyperfine  Coupling  Constants 

The  A  tensor,  as  described  previously  in  the  introductory  section  dealing  with 
EPR  spectroscopy,  determines  the  magnitude  of  the  EPR  signal  splittings  due  to  coupling 
of  electronic  and  nuclear  spins.  As  was  previously  demonstrated,  in  the  case  of  GaN  the 
hyperfine  coupling  constants  of  the  lattice  constituents  are  large  enough  to  require  a 
complete  analysis  of  the  hyperfine  splitting.  We  proceed  with  the  analysis  of  the 
hyperfine  spectra  in  GaN  by  first  examining  the  development  of  the  hyperfine  coupling 
constant  and  then  applying  this  development  to  the  particular  case  of  Ga  or  N  nuclear 
spin  coupling.  Einally,  the  determination  of  the  spin  Hamiltonian,  analysis  of  the 
resulting  eigenvalues,  and  identification  of  resonant  magnetic  field  values  is  presented. 

The  A  tensor  is  comprised  of  two  primary  components,  an  isotropic  term  (also 
called  the  {Eermi}  contact  term)  and  a  dipole  term  which  is  anisotropic  due  to  the 
asymmetry  of  the  dipole  interaction.  The  derivation  of  these  terms  and  their  development 
is  available  in  several  texts  [60]  [54]  and  will  not  be  repeated  here;  however,  the 
dependence  of  the  interaction  on  electron  spin  states  and  orbital  configurations  is 
pertinent  to  the  discussion  and  will  be  addressed. 

The  dipole  term  is  derived  from  a  quantum  mechanical  computation  of  the 
magnetic  dipole  interaction  and  is  given  by  a  matrix  with  the  elements 
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(23) 


A 


a,a 


An 


ScTn 


A  =^srl^ 


where  the  brackets  indicate  integration  over  the  spatial  variables,  and  a  and  take  on 


values  of  these  spatial  values  (such  as  xx,  xy,  yy,  etc).  The  matrix  resulting  from  the 
evaluation  of  these  terms  is  traceless  (i.e.  zero  trace)  and  symmetric,  representing  the 


tensor  coupling  of  the  angular  momentum  operators  I  and  S  ;  therefore,  we  borrow 
Atherton’s  notation  and  call  this  matrix  the  dipolar  hyperfine  coupling  tensor  (a“).  This 
matrix  can  be  diagonalized  to  recover  the  spatial  dependence  of  the  anisotropy.  Atherton 
points  out  that  the  integration  over  the  electron  wavefunction  leading  to  the  development 
of  the  matrix  values  above  exhibits  a  singularity  at  the  origin.  This  singularity  is  of  no 
consequence  when  the  electron  orbital  momentum  is  greater  than  zero  (all  but  the  s 
orbitals),  as  the  electron  wave  functions  go  to  zero  at  the  origin  much  more  quickly  than 
the  singularity.  For  s  electrons,  the  symmetry  of  the  electron  wave  function  insures  that 
the  dipole-dipole  interaction  upon  which  the  expression  is  based  is  zero.  Thus,  the 
anisotropic  portion  of  the  hyperfine  coupling  constant  arises  from  the  interactions  of 
electrons  with  nonzero  orbital  angular  momentum  (l  >  O) . 

It  was  initially  recognized  from  experimental  data  that  there  must  be  a  non-zero 
isotropic  hyperfine  coupling  term  due  to  observations  of  hyperfine  splitting  in  EPR 
spectra  of  samples  in  solution.  The  isotropic  component  of  the  hyperfine  coupling  arises 
from  the  interaction  of  the  spherically  symmetric  s  orbital  electrons  with  the  nucleus. 

The  contributions  of  these  electrons  are  deduced  following  a  quasi-classical  analysis  after 
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Atherton  [60]  which  was  first  proposed  by  Fermi.  Treating  the  nucleus  as  a  spinning 
charged  shell  of  radius  r,  angular  velocity  05,  and  charge  density  a,  the  magnetic  field 


imposed  for  r  <a  is  given  by 


2jU^o  a  _ 

B  =  -^ - Q)  = 
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Bo 

Ina^ 


Bn 


(25) 


where  the  nuclear  magnetic  moment  is  determined  by  calculating  the  field  at  r  >  a  and 
determining  the  magnitude  of  an  equivalent  point  dipole  required  to  produce  this  field. 
The  energy  of  the  interaction  between  the  electron  and  the  nucleus  at  any  point  within  the 
shell  is  approximated  by  the  product  of  the  respective  magnetic  moments  multiplied  by 
the  probability  of  the  electron  being  found  at  the  point  under  consideration.  The  electron 
wave  function  ('F(r  <  a))  is  taken  to  be  approximately  the  value  at  the  origin  (since  a  is 
very  small)  and  the  energy  is  integrated  over  the  volume  of  the  spherical  shell,  yielding 

E  =  -[^lx|'P(0f  (26) 

I  3  J  [  \27,a  j  J 


where  'F(o)  is  the  value  of  the  electron  wave  function  at  the  origin,  and  is  the  Bohr 
magneton.  When  the  magnetic  moments  are  recast  in  the  form  of  the  equivalent  spin 
operators, 

(27) 

The  terms  multiplying  the  spin  operator  dot  product  are  collected  and  termed  the 
isotropic  hyperfine  coupling  constant: 
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a=\^y«s.rA'v{of 


(28) 


Since  the  isotropic  hyperfine  constant  is  obviously  dependent  upon  the  electron 
wave  function  probability  density  at  the  origin,  it  has  also  been  termed  the  Fermi  contact 
interaction.  This  term  goes  to  zero  for  any  orbital  possessing  orbital  angular  momentum, 
since  the  electron  wave  functions  must  go  to  zero  at  the  origin;  thus  the  isotropic 
hyperfine  coupling  is  due  solely  to  the  interactions  of  the  l  —  Os  electrons.  The  total 
hyperfine  coupling  term  may  thus  be  represented  by  summing  the  isotropic  and 
anisotropic  portions: 

A  =  A“  +  a/  (29) 

where  1  is  the  identity  matrix  rather  than  the  nuclear  spin  operator. 

The  effective  hyperfine  coupling  constant  A  is  the  sum  of  isotropic  and 
anisotropic  components.  Since  for  gallium  and  nitrogen  the  unpaired  electrons  are 
expected  to  occur  in  the  p  shells,  the  first  expectation  is  that  the  hyperfine  coupling 
constant  is  comprised  of  the  anisotropic  dipole  term  alone;  however,  it  has  been 
determined  that  in  most  cases  a  hyperfine  component  exists  even  for  filled  s  shell  atoms 
due  to  screening  effects  of  the  partially  filled  outer  shells  which  shield  the  spin  up  and 
spin  down  5  states  differentially,  producing  a  non-zero  5  type  spin  density  at  the  nucleus. 
This  result  implies  that  there  will  always  be  some  admixture  of  isotropic  and  anisotropic 
components  in  the  measured  hyperfine  components  [39] 
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Spin  Hamiltonian  and  Resonance  Analysis  for  the  GaN  system 

The  spin  Hamiltonian  describing  nuclear  coupling  of  the  EPR  active  5=1/2  site  is 
composed  of  several  terms  which,  depending  upon  the  relative  magnitudes  of  the 
magnetic  field  strength  and  the  hyperfine  coupling  constant,  may  or  may  not  be  important 
to  the  analysis.  These  terms  are  included  in  the  Hamiltonian  below  and  will  be  discussed 
in  turn: 

H,^,^=l^,gBS  +  r,BI+IAS  (30) 

The  first  term,  /j,^gB  S  the  Zeeman  splitting  term  and  is  typically  the  largest  term  in 
the  Hamiltonian.  For  small  values  of  the  hyperfine  coupling  constant  (A)  or  in  the  high- 
magnetic  field  limit,  the  Zeeman  splitting  term  dominates  and  the  other  terms  may  be 
treated  as  perturbations  or  neglected. 

The  second  term,  ■  7  ,  is  the  nuclear  Zeeman  term  (the  measured  term  in 
NMR)  and  is  small  for  most  nuclei  since  7^  «  BbS  can  generally  be  disregarded  or 
treated  as  a  perturbation.  However,  when  diagonalization  of  the  full  Hamiltonian  is  called 
for,  the  inclusion  of  this  term  admits  no  additional  complexities  in  the  analysis  and  is 
typically  preserved. 

For  the  case  of  a  large  hyperfine  coupling  constant  (or  conversely  of  a  small  value 
of  the  magnetic  field)  the  I  ■  ■  S  term  must  be  included  in  the  computation.  The 

inclusion  of  the  hyperfine  coupling  term  introduces  off-diagonal  terms  into  the  matrix 
representation  of  the  Hamiltonian  because  the  operator  representation  used  to  calculate 
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the  contribution  from  the  hyperfine  coupling  constant  is  comprised  of  the  raising  and 
lowering  operators  for  both  the  nuclear  and  electronic  spins 

T  S  ^I^S^+{I_S^  +  IJ_)  (31) 

where  we  have  not  included  the  hyperfine  coupling  tensor,  implicitly  assuming  that  for  a 
single  orientation  it  may  be  represented  by  a  scalar  constant,  a,  the  hyperfine  coupling 
constant.  The  operators  7+  and  S+  represent  the  application  of  the  generic  raising  and 
lowering  operators  7+  to  the  nuclear  and  electronic  spin  states,  respectively.  The  action 
of  the  raising  and  lowering  operators  on  a  given  spin  state  is  expressed  as 


j+  j, m. )  =  + 1)  -  m .  (m .  + 1)  j, m .  + 1) 

(32) 

7_  J, m. ^  =  ^7( J  + 1)  -  tn. (nij  - 1)  j, m.  - 1) 

(33) 

where  |  j,mj^  is  an  arbitrary  spin  eigenvector  [61]. 

In  this  operator  representation  the  Hamiltonian  for  a  single  orientation  with  respect  to  the 
external  magnetic  field  becomes 

=jUB8B  S  +  r,B-  I+aI,S,  +  a{l_S,  +  I,S_  )  (34) 


This  Hamiltonian  can  now  be  represented  in  matrix  form  in  the 


basis,  which 


will  be  referred  to  as  the  |  ^  basis  for  notational  simplicity,  with  the  S  and  I 

understood.  For  the  case  of  a  spin  %  nuclei  and  a  single,  spin  V2  electron  the 
Hamiltonian  may  be  expressed  in  matrix  form  as: 
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(35) 

There  are  three  immediately  apparent  results  of  this  matrix  formulation:  1)  the 
matrix  is  still  block  diagonal,  making  an  analytical  solution  possible  without 
unreasonable  complexity,  2)  the  eigenvalues  will  be  explicit,  nonlinear  functions  of  the 
magnetic  field  strength,  and  3)  the  eigenvalues  associated  with  the  block  diagonal 
elements  will  be  associated  with  mixed  quantum  states  of  the  |  )  wavefunctions. 

The  eigenvalues  associated  with  the  matrix  above  were  originally  given  in  analytical 
form  by  Breit  and  Rabi  [62] .  The  eigenvalue  formula,  or  Breit-Rabi  Formula,  is 
expressed  in  a  convenient  form  by  [63]  as 
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Where 
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X  = 


(37) 
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where  the  parameterization  is  by  the  quantum  numbers  F  {F  =  I  +  S\  I,  and  nip  ;  these 
are  good  quantum  numbers  only  in  the  spin  coupled  low-field  region  but  are  commonly 
used  to  label  the  energy  levels  throughout  the  magnetic  field  regime.  The  A  parameter 
appearing  in  equation  (38)  is  an  effective  hyperfine  coupling  constant  representing  some 
combination  of  the  isotropic  and  anisotropic  coupling  constants.  The  nonlinear 

dependence  of  the  energy  eigenvalues  on  |.b|  appears  in  the  Breit-Rabi  formula  in  the 


treatment  of  the  x  parameter.  The  eigenvalues  F\B\]  are  shown  in  Figure  10  (page  54) 


IbI); 


for  a  spin  %  nuclei  coupled  with  a  single  spin  Vi  electron,  calculated  using  the  Breit-Rabi 
formula  above  and  parameters  appropriate  to  a  gallium  nucleus. 

The  wavefunctions  that  result  from  the  eigen-solution  of  the  Hamiltonian  matrix 
of  equation  (35)  are  properly  labeled  by  the  total  angular  momentum,  F,  only  at  zero 
applied  field  since  the  electron  and  nuclear  spins  are  completely  coupled.  In  the  high- 
field  regime,  the  electronic  Zeeman  term  dominates  and  the  resultant  wave  functions  are 
approximately  given  by  the  simple  product  wavefunctions  of  the  electron  and  nuclear 
spins  I  ) .  However,  in  the  intermediate  field  range  where  the  hyperfine  splitting 

constant  is  of  the  order  of  the  applied  magnetic  field  mixing  of  the  spin  wavefunctions  is 
appreciable.  In  this  case,  the  mixed  eigenstates  that  correspond  to  the  eigenvalues 
obtained  from  the  diagonalization  of  the  Hamiltonian  may  be  determined  by  application 
of  the  Clebsch-Gordan  coefficients  [64]  for  the  spin  coupling  of  interest,  but  common 
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practice  in  the  EPR  community  is  to  label  the  manifold  of  eigenvalues  by  their  zero-field 
spin  labels  F  and  nip. 

Mixing  of  the  eigenstates  of  the  nuclear  and  electron  spin  in  the  energy 
eigenvalue  solution  leads  to  ambiguity  in  the  determination  of  which  transitions  are 
allowed  EPR  transitions.  In  the  low  field  regime,  where  the  Hamiltonian  may  be 
approximated  by  a  diagonal  matrix  and  corrections  treated  as  perturbations,  the 
transitions  may  be  classified  as  EPR  transitions  or  NMR  transitions  based  upon  the 
following  rules: 

NMR  :  Am,  =  =  0  (39) 

EPR  :  Anij  =  0,  Am^  =  ±1  (40) 

However,  the  admixture  of  nuclear  and  electron  spin  states  in  the  eigenvalues  in  the 
intermediate  field  region  dictates  that  transitions  cannot  be  so  simply  labeled,  since  the 
labels  nil  and  nis  are  no  longer  unambiguous.  The  final  arbiter  of  which  states  are 
allowed  is  the  experimental  data;  many  instances  of  “EPR  forbidden”  transitions  occur  in 
EPR  measurements  and  the  appearance  of  “NMR  transitions”  in  EPR  spectra  has 
precendent  [47].  Eor  the  purpose  of  this  analysis,  no  possible  transitions  are  discarded 
except  those  that  clearly  fall  outside  of  the  magnetic  field  regime  of  the  experiment. 

The  magnetic  field  values  at  which  transitions  between  the  energy  eigenstates,  or 
“Breit-Rabi  levels”,  occur  is  determined  by  the  resonance  condition  with  respect  to  the 
microwave  source  energy. 
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SO  that  once  the  eigenvalues  are  determined  as  a  function  of  the  magnetic  field,  the  values 
of  the  magnetic  field  resonance  locations  )  may  be  calculated  for  each  transition  as 
shown  in  Figure  10. 


Figure  10.  Breit-Rabi  energy  levels  for  the  coupling  of  the  a  single  electron  with  the  spin 
3/2  nucleus  ^^Ga.  Resonance  locations  for  9.5GHz  microwave  photons  are  shown  as 
dashed  vertical  lines.  The  spacing  between  resonance  locations  is  non-symmetric,  as 

shown  in  the  indicated  intervals. 


The  model  outlined  above  for  determining  the  resonance  values  of  the  magnetic  field  is 
parameterized  by  three  values:  the  electronic  g-value,  the  nuclear  g-value,  and  the 
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hyperfine  coupling  constant.  The  values  of  these  three  parameters  which  produce  the 
best  fit  to  the  experimental  data  are  determined  and  from  these  values  a  comparison  with 
known  nuclear  parameters  is  made  to  identify  the  nucleus  involved  in  the  coupling. 
Other  factors  may  affect  the  fitting  of  the  experimental  data,  as  will  be  discussed  in  the 
following  paragraphs,  but  these  three  parameters  remain  the  basis  of  the  nuclear 
identification. 


The  presence  of  more  than  one  naturally  occurring  isotope  of  a  spin-coupled 
nucleus  adds  some  complexity  to  the  analysis  of  the  hyperfine  coupling  resonances.  In 
the  case  where  isotopic  species  have  differing  nuclear  spin  values,  a  different 
Hamiltonian  matrix  must  be  determined  for  each  of  the  isotopes  and  solutions  unique  to 
the  isotope  determined.  Fortunately,  in  many  cases  the  nuclear  spins  are  the  same  and  the 
spectral  contributions  of  different  isotopes  may  be  calculated  by  simply  changing  the 
parameters  associated  with  the  isotope,  primarily  the  hyperfine  coupling  constants.  In  the 
case  of  a  nucleus  having  two  naturally  occurring  isotopes  (labeled  a  and  b)  with  identical 
spin,  the  hyperfine  coupling  constants  must  be  related  by  the  relationship 


Sn 

h 

S  N 


(42) 


so  that  the  relative  magnitudes  of  the  hyperfine  coupling  constants  are  constrained  and  do 
not  introduce  any  additional  degrees  of  freedom  to  the  modeling  process.  The 
relationship  between  the  ratios  of  the  nuclear  g-values  and  hyperfine  coupling  constants 
is  a  consequence  of  the  definition  of  the  hyperfine  coupling  constants  discussed 
previously,  where  the  hyperfine  coupling  constants  were  shown  to  be  linearly  dependent 
upon  the  nuclear  g- value 
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After  the  resonance  locations  (in  magnetic  field  units)  are  determined,  the 
hyperfine  spectrum  is  modeled  by  summing  the  contributions  of  the  various  naturally 
occurring  isotopes,  weighted  by  their  isotopic  abundances.  The  resonance  locations 
calculated  by  the  analysis  method  above  represent  the  peak  locations  of  the  microwave 
absorption  spectrum,  or  the  null  locations  (zero  crossings)  of  a  first-derivative  spectrum. 
In  order  to  facilitate  simple  visual  comparison  a  Gaussian  or  Lorentzian  derivative 
function  is  placed  at  each  of  the  constituent  resonant  B  values  and  it  is  these  functions 
that  are  summed.  The  placement  and  summing  of  the  derivative  functions  admits  an 
additional  fitting  parameter  to  the  model,  the  width  of  the  underlying  distribution 
function.  It  is  not  uncommon  for  the  components  of  an  isotropic  hyperfine  spectrum  to 
display  different  linewidths  [65],  and  although  this  is  typically  observed  in  liquid 
solutions  where  tumbling  of  the  molecules  may  broaden  lines  anisotropically,  chemical  or 
physical  processes  may  exist  that  cause  anisotropic  line  broadening  in  solid  state  centers. 
The  most  conservative  approach  to  choosing  linewidths  is  to  assume  isotropic  broadening 
of  the  hyperfine  spectral  components  unless  theoretical  considerations  dictate  otherwise, 
limiting  the  added  fitting  parameters  to  a  single  linewidth  parameter.  Analyses  presented 
in  this  work  are  based  upon  a  single  linewidth  fit  unless  otherwise  noted. 

Hall  Effect  Analyses 

Hall  measurements  performed  on  GaN  epilayers  grown  upon  sapphire  substrates 
often  show  anomalous  low  temperature  behavior  characterized  by  large  measured  values 
of  the  carrier  concentration  [66] .  The  source  of  this  anomalous  concentration  is  a  highly 
conductive,  degenerate,  high  defect  density  region  at  the  GaN/sapphire  interface.  This 
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layer  is  observed  in  MBE  and  HVPE  grown  layers  as  a  consequence  of  the  lattice  spacing 
mismatch  between  the  epilayer  and  substrate  material. 

The  major  impact  of  this  degenerate  layer  is  to  mask  the  carrier  concentration 
decrease  at  low  temperatures  and  thus  skew  the  interpretation  of  shallow  donor 
concentrations,  although  Eook  [67]  has  shown  that  this  interfacial  conduction  layer 
influences  even  high-temperature  data.  Analysis  of  samples  exhibiting  this  degenerate 
layer  is  accomplished  by  either  physically  removing  the  layer  to  conduct  Hall 
measurements  on  the  remaining  material  or  by  compensating  for  the  effects  of  this  layer 
by  a  multi-layer  analysis  which  isolates  the  behavior  of  the  bulk  region.  This  model  has 
been  developed  by  Eook  [67]  and  salient  details  are  provided  here  as  they  related  to  the 
analysis  of  Hall  measurements  performed  in  this  research. 

Measurement  of  sample  conductivity  and  Hall  coefficients  as  a  function  of 
temperature  are  typically  used  to  determine  the  carrier  concentration  and  mobility  of  the 

samples  via  the  relations  =  Rh^  =  {^Rh  )  '  •  Application  of  a  multiple 

conducting  layer  model  (after  Eook)  will  allow  the  simultaneous  correction  of  the 
measured  mobility  and  carrier  concentration  data.  The  basic  relationships  for  a  multi¬ 
layer  model  are  simple  summations 

=11^'  =  'LRhj  Rhj  ='Z^Hj  j )  (43) 

j  j  j 

and 

R'(a^J  =J^(MhjRhjJ  {Rh j)  (44) 

j  j  j 

where  the  S  superscript  is  written  explicitly  to  denote  a  measured  sheet  (cm'  )  parameter. 
Eor  the  case  of  a  two  layer  model,  we  label  the  layers  as  layer  1  (the  bulk  or  epilayer)  and 
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layer  2  (the  degenerate  interfacial  layer).  Writing  the  Hall  mobility  and  concentration  in 
terms  of  Rh  and  aand  recasting  the  resulting  equations  in  the  form  of  the  summations 
above  yields  expressions  that  can  be  solved  for  the  mobility  and  concentration  in  the 
layer  of  interest  (layer  1): 


_  (Mh^H  ~ 


(45) 


(46) 


The  values  of  the  mobility  and  concentration  in  the  degenerate  interfacial  layer 
(jU2,ng2ld)  are  determined  by  taking  the  values  of  the  measured  mobility  and 
concentration  in  the  low  temperature  limit,  where  the  degenerate  layer  is  dominant.  As 
carriers  in  the  non-degenerate  bulk  layer  are  frozen  out  with  lowering  temperatures,  the 
carrier  concentration  measurements  approach  the  value  of  the  degenerate  carrier  density, 
allowing  the  determination  of  the  degenerate  interfacial  carrier  density.  The  application 
of  the  two  layer  model  presented  above  produces  mobility  and  concentration  curves 
which  are  amenable  to  conventional  donor/acceptor  model  fitting.  The  effects  of  the 
degenerate  conduction  layer  and  the  resulting  correction  on  a  representative  data  set  are 
shown  in  Figure  1 1  and  Figure  12. 
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Figure  11.  Hall  measured  mobility  data  from  sample  A342,  pre-  and  post-correction. 
The  effect  of  the  degenerate  layer  in  skewing  the  mobility  measurement  downward  is 

apparent. 


As  expected,  the  effect  of  removing  the  degenerate,  high  defect  density  region’s 
contribution  to  the  carrier  mobility  is  to  increase  the  estimate  of  the  mobility  in  the  bulk 
region.  This  result  is  compatible  with  the  existence  of  a  high  defect  density  interfacial 
region  and  a  higher  quality  bulk  region. 
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Figure  12  .  Hall  carrier  concentration  for  sample  A342,  pre-  and  post-correction.  The 
effect  of  the  degenerate  interfacial  layer  is  apparent  at  low  temperatures  (high  1/T).  The 
corrected  data  are  recognized  as  a  standard  double  shallow  donor  case. 


The  corrected  Hall  concentration  data  in  the  bulk  region,  as  shown  in  Figure  12, 
reflect  a  typical  n  versus  T  curve  for  a  non-degenerate,  two  donor  dominant 
semiconductor.  The  corrected  carrier  concentration  in  the  bulk  GaN  is  seen  to  be 
significantly  lower  at  all  temperature  values. 

After  the  Hall  data  has  been  corrected  for  the  interfacial  layer  effects,  the  donor 
and  acceptor  concentrations  are  determined  by  fitting  a  simple  donor  activation  model  to 
the  data.  The  development  of  this  model  relies  heavily  on  the  single  donor  analysis 
presented  by  Look  [68],  although  the  specific  case  of  two  or  three  single  shallow  donors 
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is  not  explicitly  covered  in  this  reference.  The  two  donor  case  is  of  interest  because 


If  the  range  of  measurement  temperatures  is  such  that  the  Fermi  level  may  be 
assumed  to  vary  by  only  a  small  amount  in  the  upper  regions  of  the  band-gap,  then  we 
may  assume  that  donors  or  acceptors  more  than  a  few  kT  from  this  range  are  temperature 
independent,  either  being  completely  ionized  or  un-ionized  depending  on  their  band-gap 
location.  Since  the  Hall  measurement  temperatures  used  in  this  study  vary  from 
approximately  20K  to  320K,  it  is  reasonable  to  treat  the  acceptors  as  being  completely 
ionized  throughout.  We  may  thus  dismiss  with  the  sum  over  k,  being  unable  to 
differentiate  between  the  k  acceptor  levels  in  any  case,  and  define  an  effective  acceptor 
concentration 
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(50) 


k  k 

The  effective  acceptor  concentration  is  thus  just  a  constant  to  be  subtracted  from  the 
carrier  concentration  over  the  temperature  ranges  considered. 

The  total  carrier  concentration  may  be  related  to  the  Fermi  energy  by 

n  =  (51) 

where  Nc  is  the  conduction  band  density  of  states,  N(.  =  h~^ ,  and  Sc  is  the 

conduction  band-edge  energy.  The  value  of  Nc  is  computed  to  be  4.98x10'"^  K^^'^cnT^ 
for  GaN  (m*  =  0.22 mg).  Defining  the  activation  energy  of  a  site  with  respect  to  the  value 
of  the  band-gap  energy  (So)  so  that  =  S^  -Sj ,  and  combining  equations  7,8  and  10, 
the  net  carrier  concentration  as  a  function  of  temperature  may  be  determined: 


^  Dm _ 

n  EaJkT  A 


(52) 


m  ^  •  2  Nr 


where  the  degeneracy  ratio  (gg  /  gi )  is  taken  as  Vi  for  the  shallow  donor  states.  The 
carrier  concentration  for  two  shallow  donors  and  an  undetermined  number  of  deep 
acceptors,  in  the  low  to  moderate  temperature  regime,  is  thus 
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where  the  donor  concentrations  {Ndi  and  Ndi)  and  activation  energies  {Edi  and  Edi)  may 
differ  for  donors  1  and  2.  Equation  12  is  parameterized  by  these  donor  concentrations 
and  activation  energies,  as  well  as  by  the  effective  acceptor  concentration  {Nf  ).  Given 
these  parameters,  n  may  be  calculated  numerically  as  a  function  of  the  temperature.  In 
order  to  develop  a  sense  of  the  dependencies,  examples  of  this  fitting  procedure  are 
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shown  in  Figure  13  for  one,  two  and  three  donor  models  with  various  acceptor 
concentrations. 


2  3.8  5.6  7.4  9.2  1  1  12.8  14.6  16.4  18.2  20 

1000/T(K-i) 


Figure  13.  Examples  of  the  donor  fitting  model,  curves  labeled  A  are  one  donor  models, 
B  are  two  donor  models,  and  C  curves  are  three  donor  models.  Curves  labeled  1  are  for 
acceptor  densities  of  l.SxlO'*’  cm~^  and  the  curves  labeled  with  2  are  for  acceptor 

concentrations  of  5  x  lO'^  cm~^ . 


Since  the  EPR  measurements  in  this  study  are  made  at  measurement  temperatures 
below  lOK,  virtually  all  of  the  donor  sites  are  in  their  neutral  charge  states,  the 
conduction  band  electrons  having  been  “frozen  out”.  This  should  hold  for  even  the 
shallowest  donors.  Since  EPR  measures  a  single  site  for  each  spectral  component,  which 
is  consistent  with  shallow  hydrogenic  donor  sites  [33],  the  shallow  donor  measured  Hall 
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state  may  tentatively  be  identified  with  the  EPR  measured  donor  state.  This  correlation 
can  be  experimentally  verified  by  comparison  of  the  Hall  shallow  donor  concentrations 
and  the  EPR  measured  concentration.  These  values  should  scale  linearly  if  the  EPR 
shallow  donor  site  is  truly  associated  with  the  shallowest  Hall  measured  state.  Under  this 
assumption,  changes  in  deeper  state  populations  should  not  affect  the  EPR  measurement 
of  the  shallow  donor  state. 
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IV.  Experimental  Procedures 


A  description  of  the  experimental  setup  and  procedures  must  include  discussion 
of  the  irradiation  experiments  as  well  as  the  spectroscopy  experiments.  Additionally, 

Hall  Effect  experiments  were  conducted  in  parallel  with  the  spectroscopy  experiments, 
and  should  be  considered  part  of  the  overall  experimental  program.  Supporting 
experimental  methodologies  include  both  sample  preparation  and  properties,  and 
dosimetry  considerations. 

Sample  Preparation  and  Handling 

The  samples  used  in  this  dissertation  were  provided  by  a  variety  of  sources, 
primary  among  which  are  Air  Force  Research  Laboratory  (AFRL)  Materials  and 
Manufacturing  Directorate  (ML).  AFRL/ML  samples  were  provided  by  Dr.  Joe  Van 
Nostrand  (A342,  A346,  A350,  A351,  and  A363)  and  Dr.  David  Look  (SB0009B). 
Samples  provided  by  Dr.  Van  Nostrand  (A###)  were  grown  in  the  late  1990’s  through 
early  2000’ s  at  AFRL/ML  via  MBE,  and  are  silicon  doped  n-type  GaN  on  sapphire 
substrates.  The  sample  provided  by  Dr.  Look  (SB0009B)  was  acquired  by  AFRL/ML 
from  Samsung  for  use  in  Hall  effect  measurements  and  is  a  220  |im  free-standing  GaN 
layer  grown  via  HVPE  [69].  This  HVPE  sample  is  smaller  than  any  of  the  MBE  samples 
and  is  the  only  sample  of  this  type  readily  available;  therefore,  experimental  procedure 
for  this  sample  is  different  from  other  samples  as  there  is  no  “control  sample”  to  compare 
against.  Several  additional  samples  were  initially  provided  by  Dr.  Van  Nostrand  with 
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18  3 

carrier  concentrations  greater  than  10  cm  .  These  samples  were  later  determined  to  be 
inappropriate  for  EPR  measurements  and  were  not  used  thereafter. 

Upon  receipt  of  samples  they  were  sent  to  the  packaging  group  of  Mr.  Larry 
Callaghan  (AFRL/ML)  to  be  diced  into  segments  of  appropriate  size  to  fit  into  EPR 
sample  tubes,  nominally  3mm  x  10mm.  Samples  were  separated  and  labeled  by 
subsection  (A342-1,  A342-2,  etc).  Different  subsections  were  used  to  provide  nearly 
identical  experimental  and  control  samples  -  typically,  one  subsection  of  a  sample  was 
irradiated  while  an  accompanying  subsection  was  preserved  in  a  virgin  state  to  compare 
against  the  irradiated  sample  under  near  identical  experimental  conditions.  An  additional 
subsection  (nominally  4mm  x  4mm)  was  sectioned  from  each  sample  for  use  in  Hall 
effect  measurements.  Sample  subsection  physical  characteristics  are  listed  in  Table  5  for 
samples  of  interest.  These  carrier  concentration  measurements  and  mobility 
measurements  were  taken  from  documents  supplied  by  Dr.  Van  Nostrand  with  the  Si 
doped  samples,  and  were  repeated  on  the  subsections  designated  as  Hall  measurement 
samples  at  a  later  date.  The  concentration  and  mobility  data  shown  in  Table  5  are  taken 
from  the  measurements  provided  with  the  samples  from  AFRL/ML. 
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Table  5.  Sample  Physical  Characteristics 


Sample 

Mass 

(gm) 

Nominal  Epitaxial 
Thickness  (/jm) 

RT  Carrier 
Concentration  (cm-3) 

RT  Mobility 
(cm /Vs) 

A342-1 

0.0182 

2.0 

2.44E-I-17 

68.5 

A342-2 

0.0393 

2.0 

2.44E-I-17 

68.5 

A346-1 

0.0522 

2.0 

3.81E-I-17 

276 

A346-2 

0.0514 

2.0 

3.81E-I-17 

276 

A350-1 

0.0379 

2.0 

2.72E-I-17 

141.9 

A350-2 

0.0178 

2.0 

2.72E-I-17 

141.9 

A351-1 

0.0344 

2.0 

3.99E-I-17 

A351-2 

3.99E-I-17 

207.5 

A363-1 

■IIlEIIgM 

6.78E-I-16 

104.9 

A363-2 

0.0383 

2.0 

6.78E-I-16 

104.9 

SB0009B-1 

Not  available 

Not  available 

SB0009B-2 

0.0114 

Not  available 

Not  available 

Samples  were  stored  out  of  direct  UV  exposure  (predominantly  sunlight)  at  room 
temperature  in  clean  sample  holders  until  irradiated.  Following  irradiation,  samples  were 
loaded  into  EPR  sample  tubes  immediately  and  thereafter  were  kept  immersed  in  liquid 
nitrogen  (LN2)  until  EPR  or  Hall  measurements  could  be  performed.  This  cryogenic 
storage  technique  was  used  in  order  to  ameliorate  annealing  effects.  Eor  both  EPR  and 
Hall  measurements,  low  temperatures  are  required  and  the  samples  were  exposed  to  the 
temperature  profiles  required  of  the  measurement  process  -  typically  very  fast  reduction 
from  room  temperature  to  approximately  4K  in  the  case  of  EPR  measurements  or 
reduction  from  room  temperature  to  approximately  20K  followed  by  a  slow  warm-up  to 
room  temperature  in  the  case  of  Hall  measurements.  In  both  cases,  the  time  at  room 
temperature  prior  to  performing  the  experimental  measurements  is  kept  at  the  lowest 
practicable  level.  Eor  EPR  measurements,  this  time  is  on  the  order  of  eight  to  15  minutes 
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total  room  temperature  exposure  prior  to  measurement.  For  Hall  measurements,  the  time 
is  longer  and  there  is  some  localized  sample  heating  due  to  the  requirement  to  apply 
contacts  to  the  samples  after  the  irradiations  are  performed. 

Following  measurement,  EPR  samples  are  replaced  in  LN2  for  cold  storage  prior 
to  additional  EPR  measurements  or  are  left  at  room  temperature  in  the  case  of  isochronal 
annealing  studies.  Samples  that  are  to  be  annealed  at  room  temperature  are  typically 
annealed  in  a  quartz  sample  tube  under  ambient  atmosphere. 

Sample  cleaning  is  accomplished  by  swabbing  with  methanol  or  acetone  (where 
required)  followed  immediately  by  rinsing  with  de-ionized  water.  Samples  were 
typically  cleaned  in  this  fashion  previous  to  first  irradiation,  when  contaminated  by  the 
vacuum  grease  or  rubber  cement  typically  used  for  sample  mounting,  or  when 
contamination  was  visible  on  the  samples. 

Irradiation  Experiments 

Irradiations  were  all  performed  at  the  Wright  State  University  (WSU)  Van  de 
Graff  (VDG)  facility  and  were  performed  by  Maj  Greene  under  the  guidance  of  Dr.  Gary 
Earlow,  of  the  WSU  faculty.  All  operators  are  trained  on  VDG  operations  and  safety,  as 
well  as  completing  the  WSU  radiation  safety  training  course. 

The  VDG  at  WSU  is  a  low-  to  mid-energy  accelerator,  routinely  operating  in  the 
500keV  to  1.8  MeV  range,  with  beam  currents  less  than  30  |lA.  Beam  uniformity  is 
typically  good,  with  a  temporal  current  deviation  estimated  at  approximately  ±3%  and  a 
temporal  energy  deviation  estimated  at  approximately  ±5%  (P-V).  These  values  are 
derived  from  operator  experience  and  observation.  Spatial  beam  uniformity  over  a  2  cm 
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by  2  cm  square  is  estimated  by  the  WSU  VDG  facility  staff  at  ±  2-3%  from  optical 
measurements  of  irradiated  plastics. 

The  irradiation  chamber  consists  of  an  evacuated  (210'^  torr)  aluminum  cylinder, 
equipped  with  a  magnetic  beam  steering  system  and  an  aperture  for  control  of  secondary 
emissions.  The  chamber  is  capped  with  a  cryogenically  cooled  vacuum  cap  which 
incorporates  the  sample  mounting  stage  as  well  as  electrical  pass-through  lines  for 
electrical  and  temperature  measurements.  In-situ  electrical  measurements  of  GaN 
samples  were  not  performed  during  any  of  the  experiments;  however,  temperatures  were 
cryogenically  controlled  and  remotely  monitored  during  irradiations  when  possible. 
Irradiations  in  this  study  were  performed  with  the  beam  parameters  listed  in  Table  6. 
Irradiation  flux  levels  were  limited  to  approximately  20  pA  (1.2x10^"^  electrons/sec)  as  an 
attempt  to  control  dose  rate  dependent  effects  [70]. 


Table  6.  VDG  Irradiation  Parameters 


Irradiation  Date  / 
Samples  Irradiated 

Beam  Current/Flux 

Particle  Energy 

Total  Fluence 

12  May  03 

A342-1,  A363-1  A350-1 
A351-1,SB0009B-1 

20  pA/  1.2x10^^  s'^ 

1.0  MeV 

1x10^^  e'  cm'^ 

17  Jun  03 

A350H1,  A351H1 

20  pA/  1.2x10^^  s'^ 

1.0  MeV 

1x10^^  e'  cm'^ 

11  Apr  03 

A342H1,  A342-1,  A351-1 

10  pA/ 6x10^^  s'^ 

1.0  MeV 

1x10^^  e'  cm'^ 

11  Mar  03 

A342-1 

10pA/6xl0‘^  s'^ 

1.0  MeV 

1x10^^  e’  cm^ 

30  Sep  02 

GaN  /  Sapphire 

20  pA/  1.2x10^^  s'^ 

1.5  MeV 

2x10^^  e'  cm'^ 

30  Jul  02 
GaN/Sapphire 

10pA/6xl0‘^  s'^ 

0.5  MeV 

1x10^^  e'  cm'^ 
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The  cryogenically  cooled  sample  stage  and  vacuum  cap  discussed  above  was 
machined  for  this  experiment  by  the  AFIT  model  shop  in  order  to  cool  the  samples 
during  irradiation  and  permit  the  temperature  and  electrical  measurements  to  be  made 
while  maintaining  the  vacuum  integrity  of  the  system.  This  sample  stage  is  shown  in 
Figure  14. 


Figure  14.  Cold  head  and  sample  mount  assembly. 


The  original  WSU  VDG  setup  used  a  sample  mount  which  was  water  cooled,  the 
mount  being  cooled  primarily  to  prevent  heat  build-up  in  the  vacuum  chamber  itself. 

Due  to  reports  of  defect  annealing  at  room  temperatures  in  GaN  [37]  and  estimates  of 
sample  temperatures  as  high  as  200°  C,  the  new  system  was  designed  to  perform  cooling 
of  the  sample  by  liquid  nitrogen  (LN2).  Further  modifications  of  the  system  were  made 
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to  insure  coolant  flow  (prevent  vapor  lock)  and  to  accurately  determine  the  sample  stage 
temperature  during  irradiation  (inclusion  of  an  embedded  resistance  thermal  device  or 
RTD).  In  the  final  configuration,  the  cold  head  achieved  minimum  sample  stage 
temperatures  of  84K,  and  was  able  to  maintain  sample  stage  temperatures  below  95K 
during  20  pA  irradiations.  Temperature  control  was  demonstrated  over  periods  of  3-6 
hours,  after  which  LN2  dewars  must  be  replaced.  Temperature  profiles  using  this  cold 
head  design,  as  measured  by  an  embedded  Resistive  Temperature  Device  (RTD)  are 
shown  in  Figure  15. 


Figure  15.  Typical  temperature  profile  for  long  irradiations.  Data  is  from  a 

18  2 

10  eVcm  irradiation  at  approximately  20  pA  beam  current  performed  on 


12  May  2003. 
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Beam  Current  (pA) 


In  a  typical  experimental  run,  the  VDG  system  would  be  “exercised”  on  the  day 
prior  to  the  irradiation  by  setting  up  and  running  the  system  at  the  desired  energy  and 
current  for  1-2  hours.  This  procedure  allows  the  tuning  of  the  energy  and  current  without 
inclusion  of  the  sample,  as  this  process  may  take  considerable  time  before  stability  is 
achieved.  After  a  stable  setup  is  achieved,  the  system  is  allowed  to  run  for  1-2  hours  in 
order  to  check  stability  and  because  the  system  becomes  more  stable  as  it  is  allowed  to 
run  at  the  desired  set  point.  This  procedure  allows  the  sample  irradiation  to  begin  at  the 
desired  setpoint  on  the  following  day,  minimizing  dose  and  dose  rate  excursions. 

After  the  exercise  of  the  system,  the  vacuum  chamber  is  opened  and  the  samples 
mounted  to  the  sample  stage.  Typically,  up  to  three  samples  (approximately  10mm  x 
3mm)  are  mounted  on  the  stage  at  one  time  to  provide  comparison  of  samples  under 
identical  irradiation  conditions.  Samples  were  often  accompanied  by  3.175mm  x 
3.175mm  TLD-400  chips  and  larger  10mm  x  10mm  (approximate)  material  samples  for 
Hall  Effect  measurements.  The  mounting  was  initially  accomplished  by  means  of 
applying  a  small  amount  of  Apiezon  Type  N  vacuum  grease.  This  vacuum  grease  is  not 
EPR  active  in  its  unirradiated  state;  however,  upon  irradiation  the  vacuum  grease  takes  on 
a  crystalline  form  and  exhibits  an  EPR  absorption  peak  (see  Experimental  Results  section 
for  details).  Due  to  EPR  interference  from  the  irradiated  vacuum  grease  signal,  the  cold 
head  assembly  was  modified  to  mechanically  mount  samples.  The  majority  of  sample 
irradiations  were  performed  without  grease  mounting  using  this  physical  clamping 
system.  Only  very  small  portions  of  the  samples  were  shielded  by  the  aluminum 
mounting  bars;  no  appreciable  effect  on  the  EPR  measurements  is  expected  due  to  the 
small  percentage  of  the  total  volume  being  shielded. 
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After  sample  mounting,  the  stage  is  affixed  to  the  vacuum  chamber  and  the 


chamber  is  pumped  down  to  the  desired  pressure  for  approximately  12-24  hours.  Before 
irradiation  on  the  following  day,  the  LN2  connections  are  fitted  and  the  instrumentation 
cabling  run  from  the  control  room  for  experiment  monitoring.  During  irradiation,  all 
personnel  are  evacuated  from  the  room  and  the  VDG  room  is  alarmed  to  prevent 
personnel  exposure  to  the  high  secondary  (x-ray)  radiation  levels. 

Beam  current  and  total  fluence  are  measured  by  current  collection  and  integration  via  the 
conducting  end-cap  and  sample  mount.  This  is  performed  by  an  electrical  connection 
which  leads  to  an  ammeter  and  signal  integrator  in  the  control  room.  This  charge 
collection  method  is  the  basis  of  the  measurement  of  electron  fluence  values. 

Irradiations  were  performed  to  the  specifications  of  irradiation  energy  and  dose 
shown  in  Table  7. 


Table  7.  Irradiation  Specifications 


Irradiation 

Batch 

Energy 

Dose 

Samples 

Purpose 

02/03  -  01 

1.0  MeV 

0.44  Mrad(GaN) 
1-10^®  eVcm^ 

A342 

EPR 

NRL/DC 

03/03  -01 

1.0  MeV 

0.44  Mrad(GaN) 
1-10^^  eVcm^ 

A342,  A350 

EPR 

NRE/DC 

04/03  -01 

1.0  MeV 

4.41  Mrad(GaN) 
1-10^’  eVcm^ 

A342,  A350, 

EPR 

NRE/DC 

05/03  -  01 

1.0  MeV 

4.41  Mrad(GaN) 
1-10^’  eVcm^ 

A342H1, 

A363H1, 

A350H1 

Hall 

Measurement 

05/03  -  02 

1.0  MeV 

4.41  Mrad(GaN) 
add  200  um  dose 
1-10^^  eVcm^ 

A342,  A346, 
A350,  A351, 
SB0009B 

EPR 

NRE/DC 
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After  the  conclusion  of  irradiation,  the  sample  is  removed  from  the  vacuum 
chamber  (after  appropriate  safety  checks  have  been  performed)  and  removed  from  the 
sample  stage.  The  samples  are  immediately  placed  into  quartz  holders  and  submerged  in 
LN2  for  transport.  Samples  are  kept  at  cryogenic  temperatures  throughout  the  transport 
and  measurement  processes  in  order  to  reduce  the  impact  of  thermal  annealing  until 
annealing  measurements  are  performed. 

Dosimetry 

All  dosimetry  during  the  irradiations  was  performed  by  means  of  current 
integration.  The  sample  stage  and  supporting  structure  are  electrically  grounded  from  the 
Van  de  Graff  generator  structure  and  beam  tube,  allowing  a  measurement  of  the  charge 
deposition  at  the  sample.  The  current  generated  from  the  electron  beam  interaction  with 
the  stage  is  monitored,  measured,  and  integrated  for  a  running  measurement  of  the  total 
beam  fluence  on  the  target  region.  TLD-400  thermoluminescent  dosimeters  were 
investigated  as  a  possible  means  of  conducting  independent  dosimetry;  however,  at  the 
high  doses  (>  10  Mrad)  used  in  this  experiment  the  TLD  response  curve  was  insufficient 
to  provide  good  dosimetry.  A  suitable  TLD  material  for  the  energy  and  dose  values  used 
was  not  identified,  leading  to  the  decision  to  calculate  deposited  doses. 

Dosimetry  calculations  were  performed  to  translate  the  measured  charge 
deposition  (which  results  in  a  measure  of  the  total  fluence,  (|))  into  an  estimate  of  the 
deposited  dose  in  the  sample  of  interest.  These  calculations  were  performed  using  the 
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TIGER  Monte  Carlo  electron  transport  codes  [71]  to  determine  deposited  energy  in  the 
material  as  a  function  of  the  total  fluence. 

Calculation  of  the  electron  stopping  power  of  the  various  materials  under  study 
was  performed  via  the  XGEN  code  (part  of  the  TIGER  package)  from  the  material 
parameters  shown  in  Table  8.  These  material  parameters  were  used  to  compute  the 
stopping  powers  shown  in  Eigure  16.  Silicon  is  computed  as  a  reference  point  for  dose 
calculations. 


Table  8.  Material  Parameters  used  for  Dose  Calculations 


Material 

Weight  Fractions 

Density  (g/cm  ) 

Si 

Si  -  1.0 

N/A 

2.33 

GaN 

Ga- 0.8327 

N- 0.1673 

6.15 

AI2O3 

A1  -  0.5293 

0  -  0.4707 

3.98 

Sample  geometries,  source  characteristics,  and  monte  carlo  parameters  (histories, 
batches,  etc)  were  input  following  the  calculation  and  tabulation  of  the  material  stopping 
powers  and  the  one  dimensional  simulation,  tiger.exe,  was  used  to  estimate  the  energy 
deposited  in  the  GaN  layer.  This  calculated  dose  value  is  dependent  upon  both  the  total 
fluence  and  the  energy  of  the  electron  beam.  Beams  were  assumed  to  be  oriented 
perpendicular  to  the  GaN  face  of  the  samples,  and  were  assumed  to  be  comprised  of 
monoenergetic  electrons.  Doses  were  computed  for  a  variety  of  beam  energies,  and  since 
the  code  output  is  in  energy  deposition  units  of  {MeV  ■  cm^  %  ■  source  particle)  the  total 
dose  is  easily  scaled  with  the  fluence. 
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Figure  16.  Electron  stopping  powers  as  calculated  by  XGEN  for  materials  of  interest. 

The  dose  deposition  profile  in  a  typical  GaN  epilayer  of  2.0  pm  on  a  sapphire 
substrate  is  shown  in  Eigure  17.  The  majority  of  the  absorbed  dose  is  clearly  in  the 
sapphire  substrate  (note  relative  thicknesses).  All  quoted  dose  values  from  this  point 
forward  are  the  GaN  absorbed  dose,  calculated  by  integration  of  the  deposited  dose  over 
the  GaN  sample  thickness  unless  otherwise  noted. 


76 


Figure  17.  Dose  deposition  profiles  for  GaN  (2  pm)  on  AI2O3  (1  mm)  for  particle 
energies  of  0.5,  1.0,  and  1.5  MeV.  Calculation  was  performed  using  200,000  histories  in 

20  batches  in  TIGER. 


After  scaling  by  the  total  fluence,  and  applying  the  proper  unit  conversions,  a  set 
of  dose  curves  (dose  in  MRad(GaN)  vs  (|))  can  be  calculated.  The  calculated  curves  for 
fluences  of  10^^  to  10^*  eVcm^  at  energies  of  0.5,  1.0,  and  1.5  MeV  for  the  2  pm  films 
discussed  above  and  the  200  pm  free  standing  GaN  sample  are  shown  in  Figure  18. 
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Figure  18.  Dose  curves  (MRad(GaN)  vs  (|))  produced  via  TIGER.  Results  for  2|im  GaN 
films  on  AI2O3  are  shown  as  solid  lines  and  filled  symbols.  Results  for  200  |im  free¬ 
standing  layers  are  shown  with  dashed  lines  and  open  symbols. 


The  values  displayed  in  Figure  1 8  are  used  as  the  dose  estimates  throughout  the 
dissertation.  All  dose  values  listed  have  been  calculated  using  the  appropriate  beam 
energies,  sample  thicknesses  and  geometries,  and  fluence  levels. 

Spectroscopy 

Spectroscopy  was  performed  using  an  X-band  (9.5  GHz)  EPR  spectrometer 
manufactured  by  Bruker  Inc.  and  provided  by  AFRE/ME.  Due  to  equipment  difficulties 
with  the  liquid  helium  (EHe)  cooling  system,  the  majority  of  the  EPR  spectroscopy  was 
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performed  on  a  similar  Bruker  EPR-300  system  located  at  Naval  Research  Laboratory 
(NRL),  DC  in  Washington  DC.  This  system  was  provided  by  Dr.  William  Carlos. 

Prior  to  use,  the  various  EPR  systems  were  calibrated  using  standard  weak  pitch 
sources  to  determine  the  stability  and  accuracy  of  the  g-value  determination.  The  g-value 
measured  was  determined  to  be  accurate  within  ±0.002  for  the  samples  examined. 

To  provide  physical  mounting  of  the  samples  within  the  microwave  cavity  and 
maintain  good  thermal  contact  between  the  sample  and  LHe  coolant  flow  standard 
Wilmad  EPR  Sample  tubes  (quartz,  4mm  ID)  were  modified  by  grinding  a  coolant  hole 
in  the  bottom  of  the  sample  tube  and  annealing.  Tubes  were  abraded  until  the  form  of 
Eigure  19  was  achieved,  and  then  annealed  at  1000°  C  for  approximately  1  hour.  This 
arrangement  allows  a  rectangular  sample  of  approximately  3mm-4mm  width  to  rest  in  the 
bottom  of  the  tube  as  shown,  while  LHe  is  free  to  flow  through  the  tube  opening  and  cool 
the  sample  directly.  Samples  mounted  in  this  manner  are  not  only  efficiently  cooled,  but 
the  samples  are  not  required  to  be  mounted  using  contact  cement,  vacuum  grease,  or 
other  substances  which  could  produce  an  EPR  signature  of  varying  magnitude. 
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Figure  19.  Sample  mount  and  cooling  flow  arrangement.  Open-ended  tube  arrangement 
shown  is  required  to  cool  samples  to  near  4K. 

Relative  intensity  measurements  require  a  spin  standard  to  provide  a  reference 
signal  that  does  not  interfere  with  the  signal  to  be  measured  and  that  can  be  measured 
simultaneously  with  the  sample  of  interest.  These  requirements  were  met  by  the  quartz 
sample  tubes  used  with  the  EPR  measurement  system.  In  particular,  each  sample  tube 
was  found  to  exhibit  a  characteristic  signal  at  g=2. 007+0. 002  which  could  be  used  to 
compare  the  relative  signal  intensities  of  samples  measured  in  the  same  tube.  The  signals 
from  these  tubes  are  independent  of  angular  orientation  and  are  not  easily  saturated  at  low 
temperatures,  thus  forming  an  ideal  spin  reference  for  relative  magnitude  comparisons. 
Sample  tubes  were  labeled  and  measurements  involving  spin  comparisons  were 
performed  in  identical  sample  tubes  at  identical  locations  within  the  microwave  cavity. 
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Signal  peak  scaling  was  performed  in  the  post-measurement  analysis  phase  using  this 
spin  standard  as  a  reference. 

The  spectrometer  is  equipped  with  a  goniometer  to  measure  the  angular 
displacement  of  the  sample  with  respect  to  a  fixed  (arbitrary)  reference  point.  Alignment 
of  the  sample  with  respect  to  the  magnetic  field  axis  is  manually  accomplished  (which 
produces  an  alignment  to  within  ±15°  typically).  Determination  of  the  sample  angle 
with  respect  to  the  magnetic  field  is  accomplished  in  the  data  analysis  phase  by 
investigating  the  symmetry  points  of  the  signal  anisotropy.  A  goniometer  was  not 
available  on  the  NRL/DC  Bruker  ESR-300;  thus,  angular  measurements  were  limited  to 
relatively  large  values  of  A0  and  were  of  lower  precision,  typically  on  the  order  of  ±  15°. 

The  sample  cavity  is  cooled  to  a  level  between  50K  and  lOOK  prior  to  introducing 
the  sample  tube  and  sample.  This  is  required  to  minimize  the  temperature  cycling  of  the 
sample  during  the  handling  and  cooldown  process.  The  sample  is  cooled  to  the  operating 
temperature  of  approximately  4K  within  30-40  minutes  during  typical  operation.  As 
sample  cooling  is  provided  in  both  the  AFRL/ML  and  NRL/DC  spectrometers  by  a  LHe 
system,  the  cryostats  associated  with  the  spectrometer  (as  well  as  accompanying  transfer 
lines)  must  be  maintained  at  vacuums  on  the  order  of  1x10’^  torr  or  below. 

Cryostat  and  transfer  line  arrangements  were  similar  for  both  EPR  systems  used, 
featuring  evacuated  cryostat  chambers  specially  configured  to  mate  with  the 
spectrometer’s  microwave  cavity.  The  cryostat  chambers  typically  house  a  thermal 
sensor  (thermocouple  or  RTD  sensor),  heater  units,  wiring  for  heaters  and  sensors,  and 
transfer  tubes  for  the  LHe  coolant  as  well  as  exhaust  gas.  Such  cavities  must  be  pumped 
to  vaccum  prior  to  each  use,  due  to  outgassing  and  moderate  leak-up  rates.  This  was 
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accomplished  by  means  of  a  turbo  pump  with  cold  trap  in  the  case  of  the  NRL/DC  setup, 


and  a  dry  turbo  pump  for  the  AFRL/ML  spectrometer. 

Typical  spectrometer  parameters  are  shown  in  Table  9  and  were  optimized  for 
detection  of  signals  due  to  impurities  or  defects  in  solid  state  materials.  Spectrometer 
parameters  were  initially  based  on  suggested  settings  from  Bruker  and  were  optimized  by 
experimentation  with  the  samples  under  consideration. 


Table  9.  Spectrometer  Parameters 


Brukei 

(AFRl 

rEMX 

L/ML) 

Dr.  Carlos  Spectrometer 

NRL/DC 

Parameter 

Value 

Parameter 

Value 

Microwave 

Frequency 

9.48  GHz 

Microwave 

Frequency 

9.51  GHz 

Microwave 

Power 

1.0-  10.0  mW 

Microwave 

Power 

1.0-  10.0  mW 

Modulation 

Frequency 

100  kHz 

Modulation 

Frequency 

100  kHz 

Modulation 

Amplitude 

3G 

Modulation 

Amplitude 

2.85  GHz 

Conversion 

Time 

40-163  msec 

Conversion 

Time 

40-163  msec 

Time 

Constant 

81  -  326  msec 

Time 

Constant 

81  -  326  msec 

Receiver 

Gain 

10^ 

Receiver 

Gain 

10^ 

Operating 

Temperature 

4.2  K 

Operating 

Temperature 

4.2  K 

The  spectrometer  parameters  of  Table  9  were  chosen  (or  recorded)  based  upon 
experimental  constraints  and  experience.  Parameters  were  chosen  based  upon  the 
rationale  and  constraints  below  in  the  following  paragraphs. 
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Microwave  frequency  is  a  characteristic  of  the  microwave  cavity,  waveguides, 
and  microwave  generation  system  and  is  thus  not  amenable  to  modification;  therefore,  it 
is  simply  measured  and  recorded. 

The  microwave  power  is  user  selectable  and  is  determined  based  upon  the 
competing  requirements  to  maximize  signal  strength  (higher  power)  and  to  avoid 
saturation  of  the  absorption  signal  at  high  powers.  Typically,  the  signals  under 
observation  in  this  research  project  were  undetectable  with  power  levels  of  less  than  1.0 
mW;  saturation  began  appear  at  powers  greater  than  10  mW. 

The  modulation  amplitude  is  set  at  the  maximum  level  that  does  not  affect  the 
lineshape  of  the  measured  signal  by  producing  distortion  or  excessive  non-symmetric 
linewidths.  Higher  amplitudes  are  found  to  substantially  increase  the  SNR  of  the 
measurement  up  to  the  point  of  the  onset  of  signal  distortion.  In  the  samples  under 
investigation,  3  Gauss  was  found  to  provide  good  SNR  without  distorting  the  lineshape. 

Conversion  time  and  the  time  constant  are  related,  in  that  the  conversion  time 
essentially  represents  the  time  the  system  spends  integrating  at  each  measurement  step, 
while  the  time  constant  is  basically  a  noise  filter  on  the  input  signal.  The  value  of  the 
conversion  time  is  based  in  part  on  the  inherent  signal  strength  of  the  sample  and  upon 
experimental  constraints.  It  was  found  that  a  ratio  of  approximately  2/1  between  the  time 
constant  and  conversion  time  produces  reasonably  smooth  absorption  signatures  without 
noticeable  smearing  or  distortion  of  the  signal. 

Receiver  gain  is  typically  set  as  high  as  is  practicable  to  avoid  signal  channel 
overflows.  The  advent  of  a  signal  channel  overflow  is  signaled  by  a  message  to  the  user 
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on  both  spectrometers  used  in  this  research  -  allowing  the  run  to  be  re -performed  with  a 
lower  receiver  gain  selected. 

Hall  Effect  Measurements 

Hall  measurements  were  performed  by  Mr.  Tim  Cooper  of  AFRL/ML  on  samples 
provided  by  Maj  Greene.  Temperature  dependant  Hall  (T-Hall)  measurements  were 
performed  on  a  Lakeshore  Model  7507  system  equipped  with  a  closed  cycle  He  cooling 
system  for  temperature  stabilization  between  15K  and  320K.  The  samples  were  kept  at 
LN2  temperature  up  until  the  process  of  applying  (or  reapplying)  contacts  was  begun. 
Indium  solder  contacts  were  applied  by  hand  by  Mr.  Cooper,  and  were  not  subjected  to 
any  additional  treatments.  Samples  used  for  Hall  measurements  typically  measured 
approximately  4mm  x  4mm  square.  Measurements  were  performed  at  5K  or  lOK 
increments  from  approximately  20K  to  320K. 
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V.  Experimental  Results 


Sapphire  (AI2O3)  Results 

Sapphire  is  a  nearly  constant  companion  to  the  materials  investigated  in  this 
experiment  due  to  its  presence  as  a  ubiquitous  substrate  material.  The  presence  of 
sapphire,  or  more  properly  of  corundum  (impurity  free  AI2O3)  [72]  is  of  interest  to  the 
EPR  spectroscopist  due  to  the  large  number  of  strong  paramagnetic  absorption  peaks 
characteristic  of  this  material. 

Attempts  to  find  reference  materials  or  published  data  on  the  properties  of 
characteristic  EPR  signals  of  corundum  proved  futile,  and  the  lack  of  reference  materials 
was  verified  by  discussions  with  experienced  researchers  in  this  field  [73].  Accordingly, 
pre-irradiation  and  post-irradiation  baseline  measurements  of  simple  corundum  substrates 
were  conducted  to  determine  their  EPR  spectra  and  behaviors. 

When  samples  were  appropriately  cooled  (requiring  modification  of  quartz 
sample  tubes  to  prevent  insulation  of  the  samples)  it  was  found  that  the  sapphire  signals 
were  typically  strongly  saturated  at  the  1.0-10.0  mW  microwave  powers  typically  used  in 
this  experimental  effort.  Under  these  conditions,  sapphire  lines  were  not  evident  in  the 
sample  spectra  at  temperatures  less  than  30K.  As  sample  temperatures  approach  25-30K 
and  above,  the  presence  of  spectral  lines  attributed  to  sapphire  is  again  noticed,  and 
sapphire  signals  may  be  present  in  a  highly  saturated  form  (low  magnitude)  at  even  lower 
temperatures  when  long  integration  times  are  used. 
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Sapphire  substrate  spectra  measured  over  a  range  of  magnetic  field  angles  (-80°- 


80°)  are  shown  in  Figure  20.  These  spectra  were  obtained  at  operating  temperatures  of 
20K-30K,  near  the  level  at  which  saturation  would  begin  to  occur.  The  most  noticeable 
feature  of  these  spectral  lines  are  the  pronounced  anisotropies  with  field  angle  rotation. 
This  strong  anisotropy  provides  a  means  of  differentiating  these  signals  from  the  GaN 
signals  of  interest.  The  signals  observed  in  the  sapphire  substrates  are  also  very  narrow, 
allowing  differentiation  from  GaN  spectral  lines  which  are  typically  wider. 


Figure  20.  Wide  field  scans  of  sapphire  substrate  materials  at  approximately  30K. 
Marked  angular  anisotropy  is  evident  in  this  series  of  scans. 
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Post-irradiation  sapphire  samples  showed  no  additional  signals  at  electron 

17  2 

fluences  up  to  10  eVcm  .  Sapphire  spectra  ceased  to  be  an  experimental  interferent  after 


initial  cooling  problems  were  addressed  by  the  tube  modifications  discussed  earlier. 
Sapphire  spectra  are  recorded  here  as  a  reference  for  future  investigations. 

MBE  GaN  on  Sapphire  Substrate  Results 
Pre-Irradiation  EPR  Spectroscopy 

Pre-irradiation  characterization  of  the  initial  high  carrier  density  samples 
conducted  at  AFRL/ML  yielded  signals  that  were  attributable  only  to  the  sapphire 
substrate  materials.  Carrier  concentrations  in  the  samples  initially  procured  and 
evaluated  from  May  2002  through  October  2002  were  greater  than  1x10  cm'  ,  which 
was  determined  later  [74]  to  be  too  high  for  measurement  of  the  shallow  donor  signal.  At 
carrier  concentrations  above  1x10  cm'  ,  the  shallow  donor  signal  is  degraded, 
presumably  due  to  screening  of  the  donor  sites  by  the  onset  of  degeneracy  [75]  and  the 
resultant  temperature  independent  carrier  concentrations.  Samples  A342,  A363,  A346, 
A350  and  A35 1  were  obtained  to  address  this  issue,  all  of  which  have  room  temperature 
carrier  concentrations  at  or  below  1x10  cm'  . 

Pre-irradiation  EPR  characterization  of  samples  A342-1,  SB0009B-1,  KD40296, 
and  A3 50-1  was  performed  in  February  2003  at  the  NRL/DC  EPR  spectrometer 
laboratory.  Prior  to  irradiation,  EPR  signals  were  determined  to  occur  at  approximately 
g=1.95  (anisotropic)  and  g=2.007  (isotropic)  in  both  GaN/sapphire  samples  and  in  free¬ 
standing  GaN  in  quartz  sample  tubes.  Signal  magnitude  and  width  were  variable  for  the 
g=1.95  signal  (hereafter  referred  to  as  SDi),  varying  between  different  samples  and 
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experimental  runs.  Additionally,  this  signal  is  slightly  anisotropic  so  that  the  signal 
location  and  width  shifted  very  slightly  depending  upon  the  sample’s  angular  orientation 
with  respect  to  the  magnetic  field  vector.  As  discussed  above,  the  shallow  donor  signal 
in  sample  A342  was  very  weak  due  to  its  high  carrier  density  (~  1.2x10'*  cm”*).  Signal 
magnitude  for  the  g=2.007  (hereafter  referred  to  as  tube  signal  or  TOi  signal)  varied  with 
experimental  runs  and  was  generally  more  uniform  in  width  than  the  SDi  signal,  in 
addition  to  being  narrower.  No  anisotropy  was  noted  for  the  TOi  signal  in  any  of  the 
cases  examined.  Representative  baseline  measurements  for  all  samples  examined  are 
presented  in  figures  Figure  30  through  Figure  35,  along  with  the  corresponding  post¬ 
irradiation  measurements. 

The  constant  g-value  of  the  TOi  signal  (within  the  limitations  of  microwave 
frequency  uncertainties  and  errors)  and  the  relatively  consistent  signal  peak  width  were, 
when  considered  along  with  the  almost  universal  presence  of  this  signal,  strong  indicators 
that  this  signal  was  due  to  a  background  source  rather  than  a  component  of  the  samples 
under  examination.  The  assumption  that  this  signal  originated  in  the  sample  tubes  used 
to  mount  the  samples  was  easily  verified  by  measurement  of  the  signal  in  a  “tube-only” 
spectrum.  The  consistency  of  the  TOi  signal  and  its  parameters  make  it  ideal  for  use  in 
normalizing  various  spectral  signals  for  comparative  analysis.  A  representative  spectrum 
of  this  tube  signal  is  shown  in  Figure  21. 

The  signal  occurring  at  g=1.95  is  identified  in  the  literature  as  a  signal  arising 
from  shallow  (effective  mass-like)  donors  in  GaN  [33]  due  to  band-structure 
considerations.  The  SDi  signal  was  found  to  exist  in  Si  doped  samples  of  GaN  on 
sapphire  substrates  as  well  as  in  nominally  undoped  freestanding  GaN  layers.  The  signal 
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resonance  locations  were  not  appreciably  different  in  the  intentionally  Si  doped  and 
nominally  undoped  samples. 
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Figure  21.  Representative  EPR  scan  of  the  signal  originating  in  quartz  sample  tube  (tube 
#2).  The  resultant  signal  is  isotropic  and  not  easily  saturated.  This  spectrum  was 

obtained  at  approximately  4K. 


Baseline  (virgin)  EPR  measurements  were  performed  on  the  remainder  of  the 
samples  (A346,  A363,  A350,  A351,  and  A342)  during  March  and  April  2003  at 
NRE/DC.  These  measurements  were  consistent  with  the  previous  set,  showing  the  SDi 
signal  and  TOi  signal  in  all  samples  examined.  The  magnitude  of  the  SDi  signal  was 
variable  across  samples  in  all  of  these  measurements  after  scaling  to  the  TOi  signal 
magnitude.  The  majority  of  these  measurements  were  performed  using  the  same  sample 
tube  (referred  to  as  tube  #2),  which  allows  consistent  scaling.  Samples  measured  using 
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different  sample  tubes  are  found  to  not  scale  consistently,  but  scale  factors  were  not 
recorded  for  other  sample  tubes.  Control  subsections  of  most  samples  were  EPR 
baselined  in  this  timeframe  for  comparison  with  the  irradiated  sample  subsections. 

Post-Irradiation  EPR  Spectroscopy 

17  2 

After  irradiation  (initially  at  1.6  MeV,  10  electrons/cm  ),  EPR  spectrometry  was 
again  performed  on  the  initial  samples  to  determine  if  any  new  EPR  spectral  lines  had 
been  induced  in  the  material  by  the  irradiation.  Irradiations  of  the  intial  batch  of  GaN 
samples  (n  >  10'* I were  conducted  during  the  last  week  of  September  2002.  The 
appearance  of  several  large  sapphire  signals  is  indicative  of  thermal  insulation  of  the 
samples  in  these  early  measurements.  Results  from  these  irradiations  showed  a  strong 
signal  imparted  at  approximately  g=2.0025  that  was  not  present  in  any  of  the  pre¬ 
irradiation  spectra.  These  signals  are  shown  in  Eigure  22  as  a  function  of  the  angle  from 
the  instrument’s  magnetic  field  vector  to  the  sample  substrate  plane.  The  stability  of  the 
signal  under  rotation  with  respect  to  (WRT)  the  magnetic  field  immediately  distinguishes 
it  from  the  native  sapphire  signals,  which  display  an  extreme  anisotropy  under  rotation. 
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Figure  22.  Spectra  from  irradiated  GaN  on  sapphire  (9-30-02).  Spectra  are  shown  for 
various  angles  (0)  with  respect  to  the  instrument  magnetic  field. 

Initially,  these  signals  were  thought  to  be  defects  induced  in  the  GaN  layer  from 
irradiation,  an  assignation  bolstered  by  their  resemblance  to  the  native  SDi  signal  and  the 
absence  of  any  corresponding  signals  in  both  the  baseline  (pre-irradiation)  measurements 
and  measurement  of  sapphire  substrates  irradiated  simultaneously  with  the  GaN  samples. 
The  signal  (shown  in  Figure  23)  persisted  after  a  thorough  cleansing  of  the  sample  with 
acetone  to  remove  external  traces  of  the  vacuum  grease  used  to  mount  the  sample  during 
irradiation. 

Measurement  of  the  temperature  dependence  of  the  signal  showed  that  the  signal 
persisted  even  up  to  near  room  temperature  with  only  minor  changes  in  the  signal 
magnitude.  This  finding  is  incompatible  with  a  signal  resulting  from  radiation  induced 
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defect  sites  in  the  GaN  layer,  due  to  the  absence  of  phonon  mediated  lattice  interactions 
which  reduce  the  state  lifetimes  and  lead  to  signal  saturation  as  temperature  is  increased 
[76].  The  signal  was  found  to  be  isotropic  with  respect  to  rotation  in  the  magnetic  field  to 
within  the  experimental  uncertainty.  This  was  also  inconsistent  with  the  expected 
spectrum  of  a  defect  associated  with  the  GaN  lattice,  where  some  anisotropy  should  be 
noted  due  to  interactions  with  the  Wurtzite  structure  lattice  [77]. 


Figure  23.  EPR  absorption  signal  imparted  to  GaN/sapphire  signals  after  irradiation  (9- 
30-02).  This  figure  shows  the  non-symmetric  form  of  the  signal,  with  a  slight  hump  on 
the  high-g  side  of  the  signal,  and  a  deeper  low-g  side,  characteristic  of  multiple 

overlapping  absorption  signals. 


In  light  of  these  discrepancies,  further  investigations  were  conducted,  including 
the  irradiation  of  sapphire  substrates  and  the  vacuum  grease  used  to  mount  the  samples 
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during  irradiation.  These  irradiations  were  performed  during  the  first  week  of  October, 

2002. 

Investigation  of  the  irradiated  vacuum  grease  and  sapphire  substrates  indicated 
that  vacuum  grease  irradiated  without  the  benefit  of  shielding  by  the  substrate  material 
tends  to  crystallize.  This  crystallization  was  found  to  be  present  on  all  of  the  GaN  and 
sapphire  samples  irradiated  during  the  study;  however,  the  crystallization  was  slightly 
more  prevalent  on  the  GaN  sample  examined  than  on  the  sapphire  samples.  After 
cleansing  with  acetone  to  remove  the  crystallization,  a  discoloration  of  the  GaN  sample 
underneath  the  location  of  the  crystallized  grease  was  noted;  this  was  not  found  on  the 
sapphire  samples.  Spectra  from  samples  of  vacuum  grease  that  had  been  irradiated  to  the 
same  levels  as  the  targets  of  interest  were  obtained  and  compared  to  the  data  of  Figure 
23. 
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Figure  24.  EPR  spectrum  of  crystallized  vacuum  grease  after  irradiation.  The  signal 
shown  here  is  identical  to  that  identified  in  earlier  GaN  samples  (see  Figure  23) 


It  is  clear  from  the  data  of  Figure  24  that  the  signal  of  interest  in  Figure  23  is  due 
to  the  irradiated  vacuum  grease.  Apiezon  Type  N  vacuum  grease  was  chosen  for  its  lack 
of  paramagnetic  signals,  and  is  well  known  as  an  EPR-suitable  mounting  and  sealing 
choice  in  the  EPR  spectroscopy  community.  There  is  no  previous  record  or  publication 
of  the  signals  that  arise  from  this  material  subsequent  to  electron  irradiation. 

Upon  recognition  of  the  nature  of  the  irradiated  vacuum  grease  signal, 
experimental  methods  were  adjusted  to  obviate  the  requirement  for  grease-based 
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mounting  and  mechanical  methods  were  developed,  as  discussed  previously  in  the 
Experimental  Procedures  section. 

After  subtraction  of  the  vacuum  grease  absorption  spectra,  no  signals  were  noted 
in  the  initially  irradiated  samples.  In  these  initial  measurements,  neither  the  shallow 
donor  signal  nor  any  signals  originating  from  the  irradiation  were  noted.  Changes  in 
samples  and  experimental  procedures  were  prompted  due  to  the  lack  of  spectral  features, 
including: 

•  low  temperature  (77K  -  lOOK)  irradiation  and  storage, 

•  acquisition  of  samples  with  lower  carrier  concentration  [78],  and 

•  modification  of  sample  tubes  to  allow  better  coolant  flow. 

None  of  the  new,  lower  dopant  density,  samples  had  measured  carrier 

concentrations  above  5x10  eVcm  ,  with  the  exception  of  A342.  Measurements  of  the 
new  samples  were  conducted  starting  in  February  2003  at  NRL/DC  due  to  the  failure  of 
the  LHe  transfer  line  on  the  AFRL/ML  spectrometer’s  cryostat  and  resultant  loss  of 
cooling  capability.  Measurement  of  the  EPR  spectra  of  these  new  samples  yielded  the 
expected  results  discussed  in  the  previous  section,  Pre-Irradiation  EPR  Spectroscopy. 

Following  preparation  of  the  lower  carrier  concentration  samples,  irradiations 
were  performed  to  approximately  1x10^^  eVcm^  at  1.0  MeV  on  the  previously  EPR 
characterized  sample  A342-1  on  10-11  March  2003  at  the  WSU  VDG  Facility.  The 
sample  was  kept  cooled  during  irradiation  (however,  no  sample  temperature 
measurements  are  available)  and  were  stored  in  an  EN2  bath  until  EPR  measurements 
were  performed  on  13-14  March  2003.  The  coolant  flow  rate  for  this  irradiation  is 
questionable,  and  it  is  possible  that  the  sample  was  not  maintained  at  cryogenic 


95 


temperatures  for  the  duration  of  the  irradiation  run.  No  temperature  profiles  for  these 
irradiations  are  available.  These  March  2003  post-irradiation  EPR  measurements  were 
performed  at  NRL/DC. 

EPR  measurements  on  the  irradiated  material,  when  compared  to  previous 
measurements  on  the  identical  sample  subsections,  were  inconclusive  regarding  a  change 
in  the  magnitude  of  the  shallow  donor  signal.  Comparison  to  un-irradiated  control 
sample  subsections  of  the  same  sample  showed  no  discernable  change  in  the  magnitude 
of  the  shallow  donor  signal,  nor  was  any  change  in  the  pre-  and  post-irradiation  EPR 
spectra  of  single  sample  subsections  noted. 

Comparisons  between  irradiated  sample  subsections  and  the  corresponding  virgin 
control  sample  subsections  are  shown  in  Eigure  25.  The  data  shown  in  Eigure  25  display 
several  unusual  characteristics:  spurious  out  of  phase  signals,  large  baseline  shifts,  and 
widely  varying  linewidths  and  lineshapes  in  the  SDi  and  TOi  spectral  lines.  Problems 
with  phase-locking  in  the  lock-in  amplifier  are  suspected  for  these  features.  After  scaling 
by  the  appropriate  TOi  signal  magnitudes  the  SDi  signal  magnitudes  of  the  scans  in  the 
irradiated  sample  are  found  to  vary  by  up  to  86%.  A  similar  analysis  in  the  virgin  sample 
exhibits  a  93%  change  in  the  SDi  measurements.  These  results  indicate  that  the  spectra 
recorded  are  inconsistent  and  should  not  be  used  for  comparative  analyses. 
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Figure  25.  Comparison  of  virgin  and  irradiated  GaN  on  sapphire 

(Sample  Aa342-l/2).  Measurements  were  conducted  on  the  same  day.  Note  that  the 
data  of  the  lower  plot  display  spurious  signals  (at  g=1.93  and  g=1.985)  as  well  as 
differences  in  TOi  linewidths  and  lineshapes.  Magnetic  field  angles  of  0°  to  90°  are 

shown  in  each  plot. 


Irradiations  of  samples  A342,  A350  and  A35 1  were  performed  on  April  1 1*, 

2003  at  a  fluence  level  of  approximately  1x10  eVcm  and  an  electron  energy  of  1.0 
MeV.  The  irradiations  were  conducted  with  an  improved  stage  design  to  allow  better 
coolant  flow;  additionally,  the  sample  stage  temperature  was  monitored.  Temperature 
excursions  from  the  desired  cryogenic  levels  were  limited  to  disassembly  and  handling 
times.  Temperature  and  beam  current  profiles  for  this  irradiation  are  shown  in  Figure  26. 
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Samples  were  maintained  at  LN2  temperatures  during  storage  and  transit,  while  EPR 


measurements  were  conducted  at  NRL/DC  on  14  and  15  April  2003. 


Figure  26.  Temperature  and  beam  current  profile  for  1 1  April  2003  Irradiation 
(^  ~  lx  10'^  e~ ! crn^  ).  Temperature  scale  is  on  the  left  axis,  beam  current  scale  is  on  the 
right  axis.  Temperature  spike  at  0.75  hours  is  due  to  loss  of  chamber  vacuum. 


The  pre-  and  post-irradiation  measurements,  as  well  as  the  irradiated  sample  comparisons 
to  virgin  control  samples  were  inconclusive  in  that  no  change  in  the  signal  magnitude 
could  be  supported  by  the  data  provided.  The  single  sample  exhibiting  variation  in  the 
SDi  magnitude  was  A35 1-1,  in  which  case  the  SDi  may  have  decreased  by 
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approximately  5-10%;  however,  the  uncertainty  in  these  measurements  is  such  that  no 
conclusive  results  could  be  determined.  Comparisons  of  the  pre-  and  post-irradiation 
EPR  measurements  for  these  samples  are  shown  in  Figure  27  and  Figure  28.  The  data 
from  sample  A342  are  once  again  plagued  with  large  baseline  fluctuations  and  noise;  the 
percentage  change  from  the  lowest  to  highest  magnitude  SDi  peaks  in  the  irradiated 
sample  is  91%.  When  the  SDi  signal  in  the  irradiated  material  is  averaged  over  available 
scans  the  magnitude  is  1.97±1.41,  demonstrating  extreme  measurement  fluctuation. 


Figure  27.  Comparison  of  virgin  and  irradiated  {(j)  ~  10'^  e  /  cm^ )  GaN  on  Sapphire 
(Sample  A342).  Measurements  conducted  under  identical  conditions,  same  day. 
Magnetic  field  angles  of  0°  to  90°  are  shown  in  each  plot. 
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The  A342  virgin  data  average  to  1 . 10+0.20  over  the  displayed  scans.  These  data  are 
considered  too  unreliable  to  use  in  comparative  analyses.  Sample  A351,  shown  in  Figure 
28,  is  considered  a  better  estimate  of  the  relative  changes  due  to  irradiation  at  these 
fluences. 
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Figure  28.  Comparison  of  virgin  and  irradiated  (^-lO'^e  !  crn^)  GaN  on  Sapphire 
(Sample  A351).  Measurements  conducted  under  identical  conditions,  same  day. 
Variation  in  the  g=1.95  signal  location  is  due  to  rotational  anisotropy.  Experimental 
difficulties  resulted  in  collection  of  only  two  good  spectra  in  the  irradiated  sample. 
Magnetic  field  angles  of  0°  to  90°  are  shown  in  each  plot.  Angular  measurements  in  this 

sample  were  only  available  at  0°  and  90°. 
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The  data  of  Figure  28  are  less  comprehensive  than  in  later  experiments  due  to 
experimental  difficulties  in  acquiring  good  sample  data  with  a  well  resolved  TOi  signal 
for  comparison.  As  a  result  only  two  high-quality  spectra  were  obtained  for  the  irradiated 
A351-1  sample,  both  in  the  B  Lc  (0°)  orientation.  Similar  difficulties  were  encountered 
in  obtaining  spectra  for  the  virgin  material,  where  two  good  spectra  were  measured  at  90° 
and  one  good  spectrum  at  0°. 

The  irradiation  fluence  that  the  samples  were  exposed  to  (1x10  eVcm'  )  is 
expected  to  be  at  the  lower  threshold  of  detectability  from  previous  work  and  theoretical 
determinations,  given  an  assumed  defect  introduction  rate  of  0. 1-0.3  cm'^  [79].  The  null 
result  prompted  the  decision  to  perform  another  irradiation  at  a  higher  fluence  level 
(approximately  1x10^^  eVcm^  at  1.0  MeV)  which  occurred  on  12  May  2003.  This 
irradiation  run  was  monitored  for  sample  stage  temperature  throughout  the  course  of  the 
run,  and  the  samples  were  held  in  an  LN2  bath  during  storage  and  transit,  as  well  as 
holding  all  irradiated  and  control  samples  in  dark  containers.  Temperature  and  beam 
current  profiles  for  the  irradiation  of  12  May  2003  are  shown  in  Figure  29. 
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Figure  29.  Temperature  and  beam  current  profile  for  12  May  2003  Irradiation 
( ^  ~  1 X 10'*  / cm^  ).  Temperature  scale  is  the  left  axis,  beam  current  scale  is  on  the 

right  axis. 


Comparison  of  irradiated  samples  with  virgin  control  samples  for  A342-1,  A350- 
1,  A351-1,  and  A363-1  was  performed  on  15  and  16  May  2003  at  NRL/DC. 
Comparisons  of  the  shallow  donor  signal  magnitude  (post-correction)  in  the  irradiated 
and  virgin  samples  following  this  irradiation  showed  remarkable  differences.  The 
magnitude  of  the  shallow  donor  signal  in  every  irradiated  sample  (for  which  a  consistent 
dose  had  been  administered)  decreased  with  respect  to  the  pre-irradiation  runs  or  the 
virgin  control  sample  magnitude.  The  decrease  in  the  shallow  donor  signal  magnitudes 
are  demonstrated  in  Figure  30  through  Figure  35  for  GaN  epilayers  on  sapphire 
substrates.  Post- irradiation  data  for  sample  A342-1  is  compared  to  its  pre-irradiation 
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state  (A342-1  virgin)  as  well  as  to  a  control  subsection  (A342-2)  due  to  difficulties  in 
obtaining  high-quality  spectra  from  subsection  A342-2.  Other  samples  are  compared  to 
the  corresponding  control  subsection. 

Most  measurements  made  on  irradiated  or  virgin  samples  showed  reasonable 
uniformity  after  system  response  corrections  were  made.  Sample  measurements  made  at 
various  angles  to  the  magnetic  field  showed  variations  of  25%  to  30%  in  the  integrated 
signal  intensity  (see  following  figures).  This  variation  was  consistent  across  samples, 
with  the  larger  magnitude  peaks  always  occurring  in  the  B  Lc  orientation. 
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Figure  30.  Comparison  of  virgin  and  irradiated  (^~  lO'^  e  I crn^)  GaN  on  Sapphire 
(Sample  A363).  Measurements  conducted  under  identical  conditions,  same  day. 
Variation  in  the  g=1.95  signal  location  is  due  to  rotational  anisotropy.  Magnetic  field 
angles  of  0°  to  90°  are  shown  in  each  plot. 
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Figure  31.  Comparison  of  virgin  and  irradiated  (^~  lO'*  e  I crn^)  GaN  on  Sapphire 
(Sample  A350).  Measurements  conducted  under  identical  conditions,  same  day. 
Variation  in  the  g=1.95  signal  location  is  due  to  rotational  anisotropy.  Magnetic  field 
angles  of  0°  to  90°  are  shown  in  each  plot. 


After  the  initial  post-irradiation  measurements  were  made  on  sample  A350, 
showing  nearly  complete  removal  of  the  shallow  donor  signal,  a  simple  annealing 
experiment  was  conducted.  The  sample  was  allowed  to  sit  unperturbed  at  room 
temperature  (293K)  for  17  hours  (overnight)  and  was  then  re-measured  by  the 
spectrometer  under  the  same  conditions  as  for  the  previous  day’s  runs.  The  sample 
showed  recovery  of  the  shallow  donor  signal,  with  a  recovery  of  approximately  25%  of 
the  original  signal  magnitude,  from  about  12%  following  irradiation.  This  annealing  is 
shown  in  Figure  32. 
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Figure  32.  Comparison  of  A350  after  17  hours  room  temperature  annealing  with  virgin 
and  irradiated  (^~  10'*  e~ Icm^)  GaN  on  Sapphire.  Variation  in  the  g=1.95  signal 
location  is  due  to  rotational  anisotropy.  The  lower  curve  is  after  annealing.  Magnetic 
field  angles  of  0°  to  90°  are  shown  in  each  plot. 
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Figure  33.  Comparison  of  virgin  and  irradiated  (^~  lO'^  e  I crn^)  GaN  on  Sapphire 
(Sample  A351).  Measurements  conducted  under  identical  conditions,  same  day. 
Variation  in  the  g=1.95  signal  location  is  due  to  rotational  anisotropy.  Magnetic  field 
angles  of  0°  to  90°  are  shown  in  each  plot. 


The  number  of  sample  subsections  on  the  sample  stage  for  the  12  May  2003 
irradiation  run  was  larger  than  for  previous  irradiations,  due  to  the  lengthy  duration  of  the 
irradiation.  While  samples  were  placed  in  positions  initially  expected  to  be  within  the 
diameter  of  the  electron  beam,  sample  A342  was  placed  at  the  extreme  outside  edge  of 
the  2  cm  beam  aperture.  Upon  examination  of  the  EPR  spectra  for  this  sample,  it  was 
noted  that  sample  A342  represented  a  deviation  from  the  EPR  signal  reductions  noted  in 
the  other  exposed  samples.  This  result  can  be  explained  by  the  partial  shielding  of  the 
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sample  by  the  beam  aperture  during  the  course  of  irradiation,  significantly  reducing  the 
fluence  on  this  sample  from  that  of  the  more  centrally  located  samples. 


Figure  34.  Comparison  of  virgin  and  irradiated  (^~  lO'^  e  I cm^)  GaN  on  Sapphire 
(Sample  A342-1/2).  Measurements  conducted  under  identical  conditions,  same  day. 
Variation  in  the  g=1.95  signal  location  is  due  to  rotational  anisotropy.  Magnetic  field 
angles  of  0°  to  90°  are  shown  in  each  plot. 
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Figure  35.  Comparison  of  virgin  and  irradiated  (^~  lO'*  e  I crn^)  GaN  on  Sapphire 
(Sample  A342-1).  Measurements  conducted  under  similar  experimental  conditions,  on 
identical  sample  subsections.  Variation  in  the  g=1.95  signal  location  is  due  to  rotational 
anisotropy.  Magnetic  field  angles  of  0°  to  90°  are  shown  in  each  plot. 


The  shallow  donor  signal  magnitude  in  the  post-irradiation  samples  was  found  to 
be  a  function  of  the  signal  magnitude  (and  thus  donor/dopant  density)  in  the  virgin 
material.  The  maximum  signal  decrease  was  noted  in  samples  A363  and  SB0009B  where 
the  signal  was  completely  removed,  a  decrease  of  approximately  88%  was  observed  for 
sample  A350  after  scaling  for  differing  sample  volumes,  and  the  minimum  signal 
decrease  was  noted  in  sample  A351,  where  the  signal  decreased  by  approximately  a 
factor  of  2.  The  relationship  between  carrier  concentration  (or  donor  density)  and  EPR 
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signal  magnitude  for  virgin  and  irradiated  samples  is  shown  in  Figure  36.  The  data  of 
Figure  36  have  not  been  scaled  by  sample  volume. 


EPR  Signal  Magnitude  vs  Sample  Carrier  Concentration 
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Figure  36.  Comparison  of  scaled,  integrated  EPR  signals  (pre-  and  post-irradiation) 
versus  post-growth  room  temperature  carrier  concentration  as  determined  by  Hall 
measurements.  Note  that  sample  A342  was  probably  under-dosed  in  this  exposure. 
Sample  SB0009B  has  been  adjusted  for  the  increased  volume  of  the  sample. 


One  wide  magnetic  field  scan  (approx  1450G-5500G)  was  performed  on  sample 
A363-1  after  irradiation  to  detect  any  new  signals  that  may  have  resulted  from  the 
irradiation.  This  scan  is  shown  in  Figure  37.  The  new  signal  imparted  just  above  4000G 
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was  scanned  repeatedly,  and  although  it  showed  up  in  several  experimental  runs,  was 
found  to  be  relatively  unstable,  making  measurement  and  analysis  difficult. 


Figure  37.  Wide  magnetic  field  scan  (s  ±  c)of  sample  A363-1  after  irradiation 
{(j)  =  10'*^  ! crn^).  The  sharp  signal  appearing  at  approximately  3400  G  is  the 

ubiquitous  tube  signal.  Magnetic  field  units  are  used  as  the  dependent  variable  due  to  the 

width  of  the  scan. 

The  data  of  Figure  37  at  first  glance  seems  to  exhibit  a  periodic  structure,  which 
was  determined  in  post-measurement  analysis  to  correspond  to  a  hyperfine  interaction 
with  a  single  Ga  atom,  with  atomic  abundances  of  approximately  39%  for  ^^Ga  and  61% 


ill 


7 1 

for  Ga.  The  broad  spectrum  of  the  Ga  hyperfine  interaction  seems  to  be  superimposed 
upon  a  second  signal  near  3100G  that  shows  no  resolved  hyperfine  structure.  A  sharp 
peak  at  approximately  3380G  is  still  in  apparent,  and  is  once  again  attributed  to  the 
sample  tube.  A  more  complete  analysis  of  the  hyperfine  structure  can  be  found  in  the 
Analysis  and  Discussion  section  of  this  dissertation.  This  same  hyperfine  structure  was 
also  determined  to  be  present  in  sample  SB0009B  (free-standing  HVPE  GaN)  after  an 
identical  irradiation  (see  following  section  for  details). 

Free-Standing  GaN  Sample  Results 

Pre-Irradiation  EPR  Spectroscopy 

Sample  SB0009B  is  a  single  sample  (no  subsectioning  was  possible)  of  relatively 
small  size.  This  fact  has  driven  much  of  the  experimental  procedure  as  it  pertains  to  this 
particular  sample.  More  virgin  baseline  spectra  were  taken  on  SB0009B  than  on  the 
other  samples  and  multiple  irradiations  were  not  performed. 

Sample  SB0009B-1  was  characterized  by  EPR  spectroscopy  on  6  Eebruary  2003 
and  14  March  2003  in  its  un-irradiated  state.  These  “baseline”  EPR  measurements  were 
performed  in  duplicate  on  this  sample  in  particular  due  to  the  lack  of  sufficient  material 
for  a  control  sample.  Therefore  the  data  presented  are  pre-irradiation  and  post-irradiation 
comparisons  on  the  same  sample,  SB0009B-1. 

This  sample  varied  in  several  aspects  from  the  GaN  epilayers  discussed 
previously,  it  is  HVPE  grown  rather  than  MBE  grown  and  the  substrate  and  interfacial 
region  have  been  removed.  The  EPR  spectra  from  this  sample  differed  in  several 
respects  from  the  MBE  grown  epilayer  samples.  The  SDi  signal  in  this  material,  while 
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retaining  the  same  g-value  and  anisotropic  behavior  as  in  the  MBE  samples,  is  markedly 
wider  in  the  HVPE  samples.  The  SDi  signal  here  is  also  present  at  slightly  higher 
temperatures,  with  the  signal  actually  being  easier  to  obtain  at  near  20K  than  at  the  4K 
temperature  where  the  MBE  sample  data  was  taken.  The  overall  integrated  intensity  of 
the  SDi  signal  is  found  to  be  much  greater  for  the  HVPE  samples  than  for  the  MBE; 
however,  the  greater  width  of  the  absorption  signal  in  the  HVPE  material  reduces  the 
apparent  magnitude  of  a  first-derivative  plot.  Increased  signal  integrated  intensity  in  this 
sample  is  attributed  to  the  greater  thickness  of  the  sample,  approximately  100  times  that 
of  the  MBE  grown  epilayers. 

Post-Irradiation  EPR  Spectroscopy 

The  results  of  post-irradiation  EPR  measurements  on  SB0009B  mirror  those  of 
the  previously  discussed  MBE  grown  samples.  The  complete  disappearance  of  the  SDi 
signal  is  noted  following  irradiation  of  10  eVcm  at  1.0  MeV,  with  measurements 
performed  at  both  4K  and  22K.  This  is  the  sole  irradiation  for  this  sample;  a  comparison 
of  the  pre-  and  post-irradiation  EPR  spectra  appear  in  Eigure  38  and  Eigure  39  for 
measurement  temperatures  of  4K  and  22K  respectively. 
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Figure  38.  Comparison  of  virgin  and  irradiated  (^~  lO'^  e  I crn^)  free-standing  GaN 
(Sample  SB0009B-1).  Measurements  conducted  under  similar  experimental  conditions, 
same  sample  subsection.  Variation  in  the  g=1.95  signal  location  is  due  to  rotational 
anisotropy.  Measurements  were  conducted  at  4K  sample  temperature. 
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Figure  39.  Comparison  of  virgin  and  irradiated  {(/)-  10'*  e  I crn^)  free-standing  GaN 
(Sample  SB0009B-1).  Measurements  conducted  under  similar  experimental  conditions, 
same  sample  subsection.  Variation  in  the  g=1.95  signal  location  is  due  to  rotational 
anisotropy.  Measurements  were  conducted  at  22K  sample  temperature. 


A  wide  magnetic  field  scan  was  performed  on  sample  SB0009B-1  after  irradiation 
to  search  for  the  appearance  of  signals  due  to  the  irradiation;  however,  baseline 
measurements  were  not  performed  over  the  entire  magnetic  field  range  scanned.  The 
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data  from  this  wide  field  scan  is  shown  in  Figure  40.  A  strong  resemblance  to  the  wide- 
field  scan  of  sample  A363-1  after  irradiation  (Figure  37)  is  immediately  obvious  in  the 
data  of  SB009B-1.  This  wide-scan  data  and  the  structure  displayed  therein  will  be 
analyzed  more  fully  in  the  Analysis  and  Discussion  section. 


Figure  40.  Wide  magnetic  field  scan  (s  ±  c)  of  sample  SB0009B-11  after  irradiation 
{(j)  =  10'*^  Icm^).  The  sharp  signal  appearing  at  approximately  3400  G  is  the 
ubiquitous  tube  signal.  Magnetic  field  units  are  used  as  the  dependent  variable  due  to  the 

width  of  the  scan. 


116 


Supporting  Experimental  Results 


Pre-Irradiation  Haii  Effect  Measurements 

Samples  were  subjected  to  Hall  measurements  previous  to  irradiation  to  determine  the 
baseline  values  of  carrier  concentration  and  carrier  mobility  in  their  as-grown  state.  Most 
sample  data  showed  some  impact  from  the  presence  of  a  conductive  interface  channel 
which  was  corrected  during  analysis  following  the  method  of  Look,  et  al  [67] .  Details  of 
the  conductive  channel  correction  are  outlined  in  the  Theoretical  Considerations  section, 
and  Hall  measurement  results  are  shown  in  Figure  41  and  Figure  42  for  virgin  materials. 
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Figure  41.  Corrected  Hall  mobilities  of  virgin  materials  used  in  the  current  study.  Good 
Hall  data  was  not  available  on  SB0009B.  Measurements  are  limited  to  the  30K-320K 
regime  due  to  the  use  of  20K  measurements  to  perform  multiple-layer  corrections. 
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The  Hall  mobilities  of  Figure  41  indicate  that  the  MBE  samples  are  all  of 
intermediate  quality  as  judged  by  the  peak  mobilities  of  less  than  350  cm  V  sec'  for  all 
samples.  All  samples  shown  have  a  typical  |Ih(T)  profile,  peaking  near  160K. 

Corrected  carrier  concentrations  of  the  MBE  grown  samples  show  that  the  range 
of  room  temperature  concentrations  is  approximately  8x10^^  cm'^  up  to  7.2x10^^  cm'^,  as 
is  shown  in  Eigure  42.  Two  donor  models  fitted  to  the  measured  Hall  concentrations 
were  used  to  determine  the  shallow  donor  densities  of  interest  in  all  samples  [20] .  These 
fitted  data  are  shown  in  Eigure  43  through  Eigure  47  for  virgin  samples;  shallow  donor 
concentrations  and  activation  energies  determined  from  fitting  of  this  data  are 
summarized  in  Table  10,  appearing  on  page  124. 
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Figure  42.  Mult-layer  corrected  Hall  carrier  concentrations  of  virgin  materials  used  in  the 
current  study.  Good  Hall  data  was  not  available  on  SB0009B.  Hall  data  is  limited  to  the 
30K-320K  regime  due  to  the  need  to  use  the  lowest  T  points  in  the  multi-layer  analysis. 
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Figure  43.  Pre-irradiation  temperature  dependent  Hall  data  for  sample  A363.  Model 
(line)  is  for  a  two  donor  fit  with  (Ndi=6.3E16,  Edi=22.9  meV,  Nd2=5.5E16, 

ED2=79.2meV,  Na=2.3E15). 


120 


Figure  44.  Pre-irradiation  temperature  dependent  Hall  data  for  sample  A350.  Model 
(line)  is  for  a  two  donor  fit  with  (Ndi=2.18E17,  Edi=8.2  meV,  Nd2=4E17,  ED2=50.2meV, 

Na=1E15). 
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Figure  45.  Pre-irradiation  temperature  dependent  Hall  data  for  sample  A35 1 .  Model 
(line)  is  for  a  two  donor  fit  with  (Ndi=3.58E17,  Edi=5.7  meV,  Nd2=6.6E17, 

ED2=47.2meV,  Na=1.0E15). 
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Figure  46.  Pre-irradiation  temperature  dependent  Hall  data  for  sample  A342.  Model 
(line)  is  for  a  two  donor  fit  with  (Ndi=1.28E18,  Edi=14.2  meV,  Nd2=5.0E16, 

ED2=79.2meV,  Na=1.0E16). 
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Figure  47.  Pre-irradiation  temperature  dependent  Hall  data  for  sample  SB0009B.  Model 
(line)  is  for  a  two  donor  fit  with  (Ndi=4.8E15,  Edi=22.0  meV,  Nd2=3.6E15, 

ED2=70.2meV,  Na=2.1E15). 


Calculated  shallow  donor  concentrations  and  activation  energies  derived  from  the 
data  of  Eigure  43  through  Eigure  47  are  summarized  in  Table  10  below. 


Table  10.  Calculated  donor  densities  and  energies  derived  from  two-donor  fit 


Sample 

Calculated 

SD  Concentration  (cm'^) 

Calculated 

SD  Activation  Energy  (meV) 

SB0009B 

4.8E15 

22.0 

A363 

6.3E16 

22.9 

2.18E17 

8.2 

A351 

3.58E17 

5.7 

A342 

1.28E18 

4.5 
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Post-Irradiation  Haii  Effect  Measurements 


Samples  A342,  A350  and  A35 1  were  subjected  to  Hall  effect  measurements 

17  2  18  2 

following  irradiation  at  10  eVcm  and  10  eVcm  .  These  data  were  fit  using  the 
procedure  described  above.  In  general,  Hall  carrier  concentrations  fell  dramatically  after 
irradiation,  as  in  Figure  48  and  Figure  49. 

In  the  higher  dose  =  lO'^  e~  I cm^ )  measurements,  the  Hall  carrier 

concentration  data  in  each  sample  show  three  different  slope  regimes,  indicating  the 
introduction  of  an  intermediate  energy  donor.  Donor  fits  for  both  samples  were 
performed  using  three  donors,  providing  a  good  match  to  the  experimental  data.  We  will 
be  primarily  concerned  with  the  shallowest  donors  for  comparison  with  EPR  data,  which 
would  not  include  this  intermediary  energy  donor  given  the  unshielded  activation  energy 
of  Eij2  ~  30-50meV  . 
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Figure  48.  Post-irradiation  temperature  dependent  Hall  data  for  sample  A350.  For  the 
(j)  =  lO'^  e~  ! crn^  case  the  donor  fit  yields  (Ndi=1.9E17,  Edi=6.0  meV,  Nd2=3.60E17, 
ED2=50meV,  Na=1.0E16).  Eor  the  ^  =  lO'^  e~  I crn^  case,  fit  parameters  are 
(Ndi=2.8E16,  Edi=10.8  meV,  Nd2=2.10E17,  ED2=27.5meV,  Na=1.0E16). 


126 


Figure  49.  Post-irradiation  temperature  dependent  Hall  data  for  sample  A350.  For  the 
(f)  =  10'^  e~  I cni^  case  the  donor  fit  yields  (NDi=3.38el7,  Edi=5.3  meV,  Nd2=7.60E17, 
ED2=50.2meV,  Na=1.0E15).  Eor  the  (j)  =  10'*^  e~  I  crn^  case,  fit  parameters  are 
(Ndi=1.65E17,  Edi=8.5  meV,  Nd2=1.5E17,  ED2=31.7meV,  Na=1.5E16).  Eow 

-5 

temperature  data  (10  /T  >25)  in  this  case  exhibited  instabilities  and  may  not  be  reliable. 


Post  irradiation  peak  carrier  mobilities  are  decreased  in  all  of  the  samples 
examined.  Mobility  reduction  is  obvious  in  the  high  dose  sample  runs,  with  the 
mobilities  reduced  by  a  factor  of  approximately  three  in  sample  A350  and  a  factor  of  two 
in  sample  A351.  As  shown  in  figures  50  and  51,  the  low  temperature  mobilities  are 
greatly  reduced,  indicating  that  defect  scattering  is  the  dominant  feature  as  expected. 
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Figure  50.  Measured  Hall  mobility  for  sample  A350,  at  three  different  radiation  doses. 


Figure  51.  Measured  Hall  mobility  for  sample  A351,  at  three  different  radiation  doses. 
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Unfortunately,  with  such  low  mobility  samples  (|l  <  400  cm  /V  sec)  the  primary 
scattering  function  is  presumably  the  crystal  structure  quality  rather  than  ionized  impurity 
scattering.  This  fact  limits  the  use  of  mobility  modeling  via  numerical  solutions  of  the 
Boltzmann  transport  equation  [79].  However,  it  is  apparent  from  the  data  of  Figure  50 
and  Figure  5 1  that  the  greatest  impact  on  the  mobility  is  in  the  low-temperature  regime 
(T<150K)  which  is  the  region  in  which  ionized  impurity  scattering  dominates,  indicating 
that  the  introduction  of  charged  scattering  sites  has  resulted  from  the  irradiation. 
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VI.  Analysis  and  Discussion 


Discussion  of  the  experimental  data  and  analyses  leading  to  the  application  of  the 
data  to  understanding  radiation  interactions  with  GaN  is  presented  in  the  following 
section.  This  discussion  and  supporting  analyses  are  presented  in  three  phases:  Analysis 
of  Hall  effect  measurements,  discussion  and  analysis  of  the  EPR  data  with  regards  to  the 
shallow  donor  signal,  and  analysis  of  the  resolved  hyperfine  structure  noted  previously. 

Hall  Effect  Measurements 

Analysis  of  the  donor  and  carrier  concentration  data  provided  via  Hall 
measurements  is  critical  to  the  proper  analysis  of  the  EPR  measurements  and  so  will  be 
presented  first.  While  the  Hall  measurement  data  is  important  to  the  EPR  analysis,  it  is 
not  the  primary  focus  of  this  research  program.  In  view  of  the  supporting  nature  of  the 
Hall  measurements,  precedence  is  given  to  analyses  which  directly  support  the 
interpretation  of  the  EPR  data.  Hall  data  analyses  which  bear  most  directly  upon  the 
interpretation  of  EPR  results  are  determination  of  the  donor  activation  energies  and  of  the 
shallow  donor  concentration.  These  data  will  be  used  to  properly  quantify  the  EPR 
observations  of  the  shallow  donor  sites  and  derive  the  underlying  relationships  governing 
radiation  interactions  with  these  sites. 

The  donor  concentrations  determined  by  the  donor  fitting  model,  outlined  in  the 
Theoretical  Considerations  section,  varied  by  roughly  two  orders  of  magnitude  across  the 
five  primary  samples  used  in  this  experiment.  Over  this  range,  activation  energies  of  the 
shallow  donors  were  determined  to  change  with  the  calculated  donor  concentration  in  a 
roughly  linear  fashion  (from  about  5meV  up  to  22  meV)  due  to  ionized  donor  screening. 
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The  magnitude  of  this  screening  effect  is  parameterized  by  the  value  of  the  empirical 
screening  parameter  (o;^  ),  which  is  constant  for  a  particular  semiconductor  material 
[66] .  The  variation  in  the  activation  energy  with  donor  concentration  is  given  by 

E^=E^,-a^N^'^  (54) 

where  Edo  is  the  activation  energy  in  the  dilute  concentration  limit.  The  value  of 
reported  in  the  literature  via  calculation  is  2.1x10“^  meV  -cm  [80];  however,  measured 
values  of  the  screening  parameter  taken  from  multiple  samples  reported  in  the  literature 
yield  values  ranging  from  approximately  2.8x10^^  meV -cm  to  3.1x10“^  meV -cm 
[81].  A  linear  regression  fit  of  the  activation  energies  versus  donor  concentrations 
measured  in  the  current  study  gives  a  measured  value  of  the  screening  parameter  of 
3.07 1 X 10“^  meV  -  cm  and  a  dilute  activation  energy  of  29.6  meV. 
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Figure  52.  Fit  of  the  empirical  screening  equation  with  data  from  samples  SB0009B, 
A342,  A363,  A350,  and  A351.  The  fit  yields  values  of  =  3.071x10“^  meV  ■  cm  and 

=  29.6  meV  . 

This  calculated  shallow  donor  activation  energy  agrees  well  with  the  values  of  the 
unscreened  activation  energy  in  several  studies  [67], [34], [80], [81]  where  reported  values 
vary  from  25meV  to  30.7meV.  The  demonstrated  agreement  with  previously  measured 
values  of  GaN  shallow  donor  activation  energy  and  screening  parameter  validate  the 
results  of  the  donor  fit  employed  in  this  study,  in  which  the  results  are  clearly  dependent 
upon  the  quality  of  the  multilayer  correction  applied  to  the  Hall  data  prior  to  the  use  of 
the  donor  fitting  model.  Both  of  these  analysis  methods  are  therefore  validated  by  the 
agreement  of  these  calculated  parameters  with  established  values.  With  the  shallow 
donor  fitting  and  multilayer  correction  technique  validated  by  measurement  of  the 
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accepted  Edo  and  Cfo,  the  results  of  these  analyses  can  be  used  with  some  confidence  that 
the  underlying  analysis  is  sound. 

The  pre-irradiation  and  post-irradiation  Hall  mobility  data  discussed  previously 
(samples  A350  and  A351,  figures  50  and  51)  were  suggestive  of  the  implantation  of 
scattering  sites  in  the  material  as  a  result  of  electron  irradiation,  although  the  initial  data 
indicate  the  actual  mobilities  are  limited  by  crystalline  faults  produced  during  the  growth 
process.  Peak  low  temperature  mobilities  decrease  by  approximately  10%  and  70%  of 
the  pre-irradiation  values  for  fluences  of  10  eVcm  and  10  eVcm  respectively.  Since 
the  bulk  of  the  decrease  occurs  in  the  lower  temperature  regime  (T<150K)  where  ionized 
impurity  scattering  is  the  dominant  scattering  mechanism  in  virgin  materials  [82],  it  is 
reasonable  to  conclude  that  this  mobility  reduction  is  due  to  an  increase  in  the  population 
of  charged  defect  centers.  Without  high  mobility  samples,  quantitative  data  regarding 
this  charged  defect  concentration  is  unavailable;  even  with  better  samples,  modeling  of 
Hall  effect  measurement  data  is  insufficient  to  determine  the  identities  of  these  centers. 
EPR  measurements  may  be  able  to  identify  some  of  the  acceptor  states  so  that  a 
quantitative  study  may  be  undertaken. 

Post-irradiation  Hall  effect  measurements  made  on  samples  A350  and  A35 1  show 
the  same  activation  energies  for  the  shallow  donors  (<20  meV)  and  generally  show  donor 
concentrations  that  drop  as  a  function  of  the  irradiation  dose.  Available  sample  points 
(A350  and  A351)  show  that  the  shallow  donor  concentration  drops  by  approximately 
5.5%  of  the  original  concentration  at  10  eVcm  incident  fluence  and  51.8%  of  the  pre- 
irradiation  concentration  at  10  eVcm  in  sample  A351;  in  sample  A350,  the 
corresponding  shallow  donor  concentration  reductions  are  12.8%  and  85.5%, 
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respectively.  In  sample  A35 1  the  reduction  corresponds  to  a  linear  function, 


approximately  1:1  with  the  incident  fluence,  and  in  sample  A351  the  dose  dependence 
appears  to  be  approximately  a  ratio  of  1:1.5.  These  linear  functions  are  shown  in  Figure 
53. 


Figure  53.  Shallow  donor  concentrations  as  a  function  of  the  radiation  fluence. 


EPR  Spectroscopy 

EPR  spectra  of  the  samples  investigated  show  a  linear  response  in  the  magnitude 
of  the  SDi  signal  with  the  calculated  shallow  donor  density.  This  relationship  is  noted 
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only  after  the  SDi  spectra  have  been  scaled  by  the  TOi  signal  and  corrected  for  the 
volume  differences  between  samples.  After  these  corrections  are  applied  and  the  data  are 
compared  with  the  calculated  shallow  donor  concentrations  (not  the  overall  carrier 
concentrations)  in  the  virgin  material,  the  relationship  between  the  measured  values  of 
SDi  integrated  area  and  donor  concentration,  as  determined  by  modeling  of  the  Hall 
effect  data,  are  given  by 

=  (2.262  X 10  crn^  )n (55) 
where  value  of  the  proportionality  constant  varies  with  the  scaling  factors  used,  being 
different  if  a  different  mass  scaling  factor  were  chosen.  The  correlation  coefficient  (R  ) 
of  this  fit  is  0.9924,  indicating  an  extremely  linear  relationship.  Supporting  data  for  this 
relationship  are  shown  in  Figure  54.  The  quality  of  the  linear  relationship,  expected  from 
a  theoretical  standpoint,  is  an  indication  that  the  corrections  and  scalings  applied  to  the 
data  have  not  resulted  in  large  errors.  Similar  linear  dependences  on  carrier  concentration 
in  the  EPR  spectra  of  Si  doped  GaN  has  been  noted  in  [83]. 
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Figure  54.  SDi  scaled  EPR  magnitude  plotted  against  the  calculated  shallow  donor 
concentration.  Least  squares  fit  is  shown  as  dotted  line. 


Measurement  uncertainty  in  the  EPR  magnitudes  was  found  to  increase  with 
increasing  SDi  signal  magnitude.  This  is  not  unexpected,  since  the  scaling  process 
involves  division  by  the  integrated  area  of  the  TOi  signal,  which  is  relatively  constant  in 
magnitude;  therefore,  as  the  SDi  signal  grows  in  magnitude  the  spectrometer’s  signal 
channel  response  scales  all  of  the  data  down  to  avoid  signal  channel  overflows  and  the 
magnitude  of  the  TOi  signal  is  reduced  proportionally  As  the  SDi  signal  magnitude 
increases  and  the  relative  magnitude  of  the  TOi  signal  is  decreased,  the  quality  of  the  fit 
used  to  estimate  the  peak’s  integrated  intensity  suffers,  increasing  the  uncertainty  in  the 
magnitude  of  the  resulting  scaled  data. 
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The  demonstrated  linear  relationship  between  the  EPR  signal  magnitude  and  the 
shallow  (29  meV)  effective  mass  donor  concentration,  as  determined  by  Hall 
measurement,  at  29  meV  supports  the  contention  that  the  EPR  SDi  signal  is  associated 
with  this  Hall  donor  site.  If  the  SDi  site  were  associated  with  some  conglomeration  of 
relatively  shallow  donor  sites,  the  observed  linear  relationship  would  be  expected  to 
reflect  the  changing  ratios  of  the  various  donor  concentrations.  The  ratios  of  the  two 
shallowest  donors  varied  from  0.75  to  1.84  in  the  samples  examined;  more  importantly, 
the  sum  of  both  donor  concentrations  in  sample  A35 1  results  in  total  donor 
concentrations  over  10  cm'  ,  which  is  incompatible  with  the  good  quality  EPR  spectra 
taken  on  this  sample. 

The  single  undoped  HVPE  sample  is  presented  in  the  previous  data  set;  it  falls 
neatly  on  the  linear  fit,  although  at  the  lower  donor  densities  (~  10^^  cm'^)  that  this 
sample  represents  the  fit  quality  is  relatively  insensitive  to  large  relative  changes  in  the 
donor  concentration.  Still,  the  inclusion  of  the  HVPE  material  in  this  linear  series  is 
suggestive  of  an  underlying  relationship.  While  the  donor  type  in  the  HVPE  sample  is 
certainly  different  than  the  Si  observed  in  the  MBE  samples,  the  unscreened  activation 
energy  in  this  sample  is  near  30  meV,  consistent  with  either  the  nitrogen  vacancy  defect 
or  the  oxygen  impurity  at  30-35  meV  [84].  In  either  case,  the  EPR  spectra  of  the  shallow 
donors  vary  only  in  the  observed  linewidth  and  the  assignment  of  the  SDi  signal  to  a 
single  shallow  donor  site  may  be  assumed. 

Having  demonstrated  the  linear  relationship  of  the  EPR  signal  magnitude  and  the 
corresponding  donor  concentration,  analysis  of  the  post-irradiation  EPR  data  is 
simplified.  The  observed  reduction  of  the  SDi  signal  can  now  be  confidently  linked 
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directly  to  a  corresponding  drop  in  the  concentration  of  the  29.6  meV  donor  population. 
Reduction  in  the  SDi  signal  intensities  in  all  samples  following  irradiation  at  an  electron 
fluence  of  10  eVcm  was  demonstrated  in  the  Experimental  Results  section,  resulting  in 
signal  decreases  ranging  from  approximately  50%  in  the  higher  donor  concentration 
samples  to  complete  disappearance  of  the  SDi  signal  in  the  lower  donor  density  samples. 
This  would  indicate  that  the  shallow  donor  concentration  has  fallen  below  the  10  -10 
cm'  EPR  detection  threshold. 

These  EPR  results,  particularly  those  in  the  native  GaN  samples,  immediately 
indicate  that  the  native  shallow  donor  at  g~1.95  cannot  be  due  to  a  simple  defect  such  as 
the  nitrogen  vacancy  (Vn)  as  proposed  by  earlier  authors  [18,85  ,86],  since  the  primary 
effect  of  the  irradiation  is  to  induce  vacancy  -  interstitial  pairs  on  both  sublattices 
simultaneously  [87].  The  production  of  a  defect  primarily  responsible  for  the  native 
shallow  donor  signal  (Vn)  would  increase  the  magnitude  of  the  signal,  regardless  of  the 
intentional  doping,  whereas  the  signal  magnitude  is  shown  to  decrease  in  every  instance. 
This  finding  is  in  harmony  with  more  recent  theoretical  studies,  which  point  toward  the 
inclusion  of  some  sort  of  impurity  doping  rather  than  a  native  defect  [88]  and  is  one  of 
two  recent  studies  which  offer  experimental  evidence  for  this  theory  [89].  The  HVPE 
sample  results  srongly  support  this  assertion,  since  in  the  MBE  Si  doped  material  the 
dopant  donor  states  could  conceivably  mask  the  effect  of  the  native  donor  states;  no  such 
process  is  possible  in  the  undoped  n-type  HVPE  sample.  This  result  adds  to  the  growing 
consensus  that  n-type  auto-doping  is  due  to  impurities  rather  than  defects. 

The  decrease  of  the  shallow  donor  signal  points  to  the  interaction  of  the  shallow 
donor  in  both  doped  and  intrinsic  GaN  (whether  intentionally  doped  Si,  or  unintentionally 
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introduced  impurity  donors  such  as  O)  with  a  direct  radiation  product  such  as  vacancies 
or  interstitials.  As  discussed  previously,  this  interaction  can  take  the  form  of  either 
compensation  of  shallow  donors  due  to  the  formation  of  stable  acceptor  traps  due  to  the 
irradiation  or  direct  interaction  of  the  donors  and  radiation  induced  defects  via  the 
formation  of  donor-defect  complexes  by  the  oppositely  charged  donor  and  mobile  defect 
sites.  Obvious  candidates  for  this  complexation  interaction  are  interstitials  (Gai  or  Ni) 
due  to  the  higher  mobility  of  these  defects  in  comparison  to  the  corresponding  vacancies 
[20].  The  results  of  this  study,  when  combined  with  the  results  of  lower  energy 
irradiations  conducted  roughly  in  parallel  by  Look  [20]  point  to  some  combination  of 
donor-defect  complexation  and  acceptor  trap  compensation.  Look’s  low  energy 
irradiations  (~  300  keV)  were  designed  to  produce  nitrogen  sublattice  damage  only;  Hall 
effect  measurements  in  both  this  study  and  Look’s  experiment  show  similar  decreases  in 
shallow  donor  populations.  Single  sample  Hall  data  alone  are  unable  to  unambiguously 
determine  whether  the  shallow  donor  populations  are  being  complexed  or  compensated; 
while  a  series  of  irradiations  on  samples  with  different  donor  concentrations  should  be 
able  to  identify  if  the  complex  formation  process  is  dominant,  there  is  no  obvious  means 
to  determine  the  relative  contributions  of  these  processes  to  the  donor  concentration 
reduction  in  the  case  where  both  processes  are  present. 

Annihilation  of  the  shallow  donor  EPR  signature  is  indicative  of  the  direct 
passivation  (transition  to  a  diamagnetic  deep  state)  or  compensation  of  the  shallow  donor 
sites,  ruling  out  the  possibility  of  a  secondary  impurity  being  forced  out  of  the  lattice,  as 
was  suggested  as  a  less  likely  possibility  by  Look  [20].  If  a  secondary  impurity  donor 
such  as  H  were  present,  and  the  donor  concentration  were  dropping  due  to  the  expulsion 
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of  the  donor  from  the  lattice  during  the  irradiation  process,  then  the  strong  linear 
dependencies  of  the  EPR  data  with  the  fluence  and  donor  concentration  would  not  be 
expected  and,  more  importantly,  the  observed  room  temperature  annealing  behavior 
would  not  occur.  The  observed  post  irradiation  data  of  Figure  55  on  page  143  indicate 
that  passivation  by  a  radiation  induced  acceptor  state  must  occur,  since  the  post¬ 
irradiation  measurements  indicate  a  constant  level  of  donor  reduction  independent  of  the 
initial  donor  concentration.  It  is  important  to  understand  that  this  constant  decrease  is  a 
signature  of  acceptor  compensation  but  does  not  rule  out  the  presence  of  defect-donor 
complex  formation. 

A  plausible  model  for  this  donor  passivation  posits  the  complexation  of  the  single 
donor  state  (O  or  Si)  with  the  nitrogen  vacancy,  which  is  expected  to  be  a  single  acceptor 
[90],  accompanied  by  compensation  of  the  remaining  donor  sites  by  acceptors  not 
involved  in  the  donor-defect  complexation  process.  While  this  complexation  presumably 
occurs  at  the  irradiation  temperatures  of  80-100K  where  many  of  the  shallow  donors  are 
ionized,  if  the  complexation  process  reduces  the  carrier  concentration  by  transforming 
these  donor  sites  into  deep  states  this  change  should  be  preserved  at  the  EPR 
measurement  temperatures  of  4K-10K.  The  rate  constant  associated  with  the  formation 
of  defect  complexes  is  expected  to  be  a  function  of  the  temperature  in  two  regards:  higher 
temperatures  should  lead  to  higher  defect  mobilities,  and  at  temperatures  sufficient  to 
ionize  the  donor  states  a  coulombic  attraction  exists  between  the  donor  and  acceptor  sites. 
This  introduces  the  possibility  that  the  contribution  of  the  compensating  acceptor  sites 
may  be  determined  by  performing  a  series  of  irradiations  at  different  irradiation 
temperatures  to  vary  the  value  of  the  complexation  rate  constant.  This  approach  would 
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presumably  require  lower  temperatures  than  are  currently  available  during  irradiation  - 
down  to  4.2K  or  lower  to  prevent  any  complexation  from  occurring. 

Regardless  of  the  relative  proportions  of  the  two  competing  processes  involved  in 
annihilating  the  shallow  donor  sites,  when  both  compensating  centers  and  complexes  are 
present  the  decrease  in  the  donor  concentration  scales  in  a  one  to  one  fashion  with  the 
total  defect  concentration  up  to  the  saturation  point.  The  difference  between  the  pre¬ 
irradiation  and  post- irradiation  data  is  thus  equal  to  the  total  defect  concentration  as 
discussed  earlier.  In  the  case  of  the  10  eVcm  irradiation,  the  decrease  is  determined  to 

17  3  1 

be  2.09x10  cm'  ,  which  implies  a  defect  introduction  rate  of  approximately  0.2  cm'  . 
This  compares  favorably  with  previously  measured  values  using  only  Hall  measurements 
and  with  theoretical  values  [6].  When  Hall  measurements  alone  are  considered, 
performed  on  samples  irradiated  at  fluences  of  lO'^  e~  /  cm^ ,  3x10'^  /  cm^ ,  and 

10'*  e~  I  crn^  a  damage  constant  of  approximately  0.145  cm''  is  determined,  in  relatively 
good  agreement  with  the  EPR  results. 

While  the  data  at  the  time  of  this  writing  are  not  overwhelming,  the  observed 
linear  post-irradiation  trend,  coupled  with  the  observation  that  the  slopes  of  the  pre-  and 
post-irradiation  lines  are  nearly  identical,  disputes  the  attribution  of  the  donor  decrease  to 
complex  formation.  For  reasonable  values  of  a  (the  complex  formation  rate  constant) 
there  should  be  some  differentiable  departure  from  the  slope  of  the  pre-irradiation  curve. 
Instead  of  this  behavior,  the  slopes  of  the  high-dose  exposure  show  an  almost  exact 
correspondence  to  the  pre-irradiation  values.  This  is  in  contrast  to  Look’s  explanation  of 
the  donor  concentration  drop  as  being  due  solely  to  complex  formation  [20] ;  however,  the 
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irradiations  for  these  two  studies  were  conducted  at  different  temperatures  (80K  in  the 
present  studies,  over  300K  in  Look’s  study)  which  could  drastically  impact  the  value  of 
the  a  parameter.  It  is  quite  possible  that  our  measurements  exhibit  incomplete  complex 
formation  and  the  resulting  acceptor  compensation,  while  the  complex  formation  process 
at  higher  temperatures  proceeds  much  more  quickly,  involving  the  great  majority  of 
radiation  induced  defects  in  the  complex  formation  process.  Since  Look’s  irradiations 
were  performed  at  lower  energies  (-350  keV)  it  is  also  possible  that  we  are  seeing  a 
greater  contribution  from  the  action  any  acceptors  produced  by  Ga  sublattice  damage. 

Our  data  also  indicate  some  level  of  room  temperature  annealing,  leading  to  the 
belief  that  the  nitrogen  interstitial,  which  should  be  much  more  mobile  than  the 
corresponding  vacancy  and  has  been  shown  to  readily  form  complexes  in  GaN  [91],  is  a 
primary  source  of  any  complex  formation.  Additionally,  this  defect  is  expected  to  act  as 
a  single  acceptor,  and  the  recombination  of  vacancy-interstitial  pairs  at  relatively  low 
temperatures  (300K)  would  reduce  the  compensating  acceptor  population  as  well. 
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Figure  55.  Pre-  and  post- irradiation  values  of  the  scaled  EPR  signal  intensity  as  a 
function  of  the  measured  Hall  carrier  concentration.  Black  stars  represent  the  post¬ 
irradiation  value  predicted  for  o  =  2.34x10  'Vm^.  The  point  labeled  A351-lI(Apr)  is 
the  only  irradiated  data  point  at  ^=10'^  e~  /cm^ ,  all  other  irradiated  values  are  at 

^=10‘V“ /cm^ 

Data  from  sample  A35 1  also  validate  the  assumption  that  the  ESR  signal 
dependence  on  fluence  is  governed  by  equation  (20)  which  predicts  a  linear  relationship 
between  fluence  and  shallow  donor  annihilation  for  low  recombination  rates.  Reduction 
of  the  EPR  SDi  signal  with  increasing  doses  is  demonstrated  by  the  data  shown  in  Eigure 
56,  where  the  SDi  peaks  are  isolated  and  compared  for  sample  A351  for  three  values  of 
absorbed  dose. 
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Figure  56.  SDi  spectra  in  sample  A351  at  fluences  of  Oe  /cm^,  10'^  e  I  crn^ ,  and 
lO'^  e~  I cm^ .  Data  marked  with  filled  dots  was  collected  in  April  2003  and  the  data 
marked  with  open  squares  was  collected  in  May  2003.  Peaks  have  been  arbitrarily  placed 
along  the  x  axis  by  varying  amounts  to  ease  comparison. 


The  data  of  Figure  56  indicate  that  the  SDi  signal  magnitudes  decrease  in  a 
regular  fashion  with  increasing  electron  dose;  however,  the  linear  relationship  is  not 
immediately  obvious  since  the  data  shown  are  the  measured  absorption  derivative 
functions.  When  the  integrated  peak  magnitudes  are  examined,  the  data  for  electron 
fluences  of  0,  iO  ,andi0  eVcm  indicate  that  a  linear  relationship  between  fluence  and 
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EPR  magnitude  exists  with  slope  defined  by  the  radiation  damage  constant  and  the 
spectrometer  response  function.  This  dependence  is  illustrated  by  the  linear  fit  of  Figure 
57,  which  has  an  value  of  0.9898. 


Figure  57.  Finear  fit  to  the  dose  dependence  of  the  shallow  donor  passivation  in  sample 
A351.  The  R^  value  of  the  linear  fit  is  0.9898. 


The  passivation  process  is  also  observed  in  the  Hall  effect  measurements  as  well, 
where  the  resulting  drop  in  the  carrier  concentration  is  seen  to  decrease  after  irradiation. 
In  samples  A350  and  A351  the  percentage  decrease  in  the  carrier  concentration  measured 
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by  Hall  effect  is  closely  related  to  the  signal  magnitude  reduction  in  the  EPR 


measurements,  as  is  shown  in  Figure  58. 


Figure  58.  Comparison  of  EPR  SDi  magnitudes  and  Hall  concentration  in  A351  and 

18  2 

A350,  before  and  after  irradiation  to  10  eVcm  .  The  ratio  of  post-irradiation  to  pre¬ 
irradiation  values  is  shown  over  the  irradiated  column  in  each  set. 


For  sample  A350,  the  Hall  carrier  concentration  due  to  the  29.6  meV  shallow 
donor  decreases  by  approximately  88.4%  following  irradiation  at  a  dose  of  10^^  eVcm^ 
whereas  the  EPR  SDi  signal  magnitude  decreases  by  91.8%.  In  sample  A351,  the 
correspondence  is  even  closer,  with  the  Hall  carrier  concentration  due  to  the  shallow 
donor  dropping  by  53.9%  and  the  SDi  signal  decreasing  by  55.1%.  This  correlation  is 
apparent  only  after  the  proper  scaling  is  applied  to  the  EPR  data  and  the  shallow  donor 
concentration  is  extracted  from  the  Hall  measurements.  This  is  one  more  indication  of 
the  identification  of  the  SDi  signal  as  the  30  meV  donor.  The  close  degree  of  correlation 
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in  these  measurements  is  not  surprising  in  a  theoretical  sense,  since  these  values  must 
scale  together;  however,  in  an  experimental  sense  it  is  gratifying,  but  somewhat 
surprising,  to  observe  such  a  close  correlation  following  the  corrections  and  scaling 
performed  independently  on  each  data  set. 

The  conclusion  to  be  drawn  from  the  above  analysis  is  that  the  primary  effect  on 
the  carrier  concentration  of  the  sample  due  to  irradiation  at  these  levels  is  to  simply 
passivate  by  complex  formation  and/or  compensation  of  the  previously  active  shallow 
donors,  whether  they  are  a  shallow  implanted  donor  or  an  unintentionally  introduced 
impurity  donor.  Introduction  of  additional  donors  are  viewed  as  second  order  effects  in 
relation  to  the  magnitude  of  this  passivation  process. 

Evidence  of  a  third  donor  in  the  Hall  data  occurs  at  the  highest  irradiation  fluence 
in  samples  A350  and  A351.  Analysis  of  these  data  points  yields  an  estimate  of  the 
unscreened  activation  energy  of  approximately  107-117  meV,  which  could  be  associated 
with  the  nitrogen  vacancy  at  around  100  meV  (see  Figure  7,  page  23).  This  assessment  is 
only  a  suggestion;  more  data  will  be  required  to  demonstrate  that  these  are  the  nitrogen 
vacancy  donor  sites. 

As  alluded  to  previously,  a  single  annealing  data  point  was  collected,  representing 
the  degree  of  SDi  signal  recovery  over  17  hours  of  room  temperature  exposure,  shown  in 
Figure  59.  This  data  point  is  clearly  not  intended  to  be  an  exhaustive  investigation  of  the 
annealing  behavior,  but  is  indicative  of  the  involvement  of  some  relatively  mobile 
radiation  induced  defect  or  defects  in  the  passivation  process.  The  SDi  signal  intensity  is 
demonstrated  to  recover  to  approximately  24.1%  of  the  virgin  signal  strength,  an 
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Figure  59.  Measured  annealing  of  the  radiation  effect  in  sample  A350.  Annealing  was  a 

room  temperature  over  a  17  hour  period. 

increase  of  1.97  times  the  pre-annealing  signal  strength  of  12.2%. 

Hyperfine  Structure  Analysis 

After  noting  marked  similarity  between  the  data  of  Figure  37  and  Figure  40  and 
defect  hyperfine  structure  previously  measured  by  EPR  in  GaAs  [92]  [93],  a  hyperfine 
coupling  analysis  was  conducted.  After  initially  fitting  the  spectra  presented  in  these 
reference  papers  with  the  Breit-Rabi  model  to  confirm  its  accuracy  as  implemented,  the 
spectra  observed  in  this  study  were  fitted  using  the  same  general  technique. 

Analysis  of  resolved  hyperfine  structure  begins  with  a  determination  of  the 
appropriate  spin  Hamiltonian  and  proceeds  by  solving  for  the  resulting  eigenvalues. 
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allowable  transitions,  and  accompanying  resonant  frequencies.  In  many  cases  the  Breit- 
Rabi  formula  [62]  can  be  used  to  determine  the  solutions  of  the  eigenvalue  problem. 

The  development  and  discussion  of  this  model  were  presented  earlier,  and  at  this  point 
only  the  application  of  this  model  to  the  observed  experimental  data  will  be  discussed. 

Energy  levels  of  the  spin  Hamiltonian  are  computed  via  the  Breit-Rabi  formula 
using  the  nuclear  g-values  and  nuclear  spins  for  the  two  stable  isotopes  of  Ga,  ^^Ga  and 

7 1 

Ga,  which  occur  with  isotopic  abundances  of  60.1%  and  39.9%  respectively.  The 
nuclear  parameters  of  these  isotopes  were  presented  in  Table  4  on  page  33.  As  discussed 
previously,  the  published  g-values  of  the  isotopes  impose  an  additional  constraint  on  the 
data  analysis,  in  that  for  the  ^^Ga  /  ^^Ga  case 


S  N 


0.787 


(56) 


so  that  the  hyperfine  coupling  constants  associated  with  these  isotopes  must  occur  in  this 
constant  ratio,  no  matter  their  values.  Using  the  nuclear  properties  of  these  Ga  isotopes, 
the  energy  eigenvalues  calculated  by  equation  32  (the  Breit-Rabi  formula)  are  computed 
as  a  function  of  the  external  magnetic  field  strength.  These  eigenvalues  are  plotted  in 
Figure  60  for  gallium  hyperfine  coupling  with  a  single  electron.  The  transition  energy  of 
the  microwave  frequency  is  calculated  (e  =  h  (9.5  GHz)=6.295  x  10“^"^  J  )  and  resonant 
values  of  the  magnetic  field  are  found  numerically,  as  denoted  by  the  vertical  arrows  in 
Figure  60. 
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Figure  60.  Breit-Rabi  levels  modeled  for  and  ^^Ga.  Vertical  arrows  represent  the 
resonant  magnetic  field  strengths  corresponding  to  the  resonant  microwave  energy 

(9.50  GHz) 


Resonant  field  locations  having  been  identified  numerically,  the  model  is 
completed  by  summing  of  Gaussian  derivative  functions  centered  at  each  of  the  resonant 
field  values.  These  Gaussian  derivative  functions  are  summed  for  each  of  the  isotopic 
constituents  to  form  a  spectrum  associated  with  that  isotope,  and  the  resulting  isotopic 
spectra  are  weighted  by  the  appropriate  isotopic  abundances  and  summed  as  shown  in 
Figure  61  to  produce  the  total  absorption  spectrum. 
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B  (Gauss) 


Figure  61.  Model  output  showing  placement  of  equal  width  Gaussian  derivative 
functions  at  the  calculated  resonant  field  values.  Curves  are  shown  for  the  transitions  due 
to  ^^Ga  transitions  (lower  curve),  ^^Ga  transitions  (middle  curve),  and  the  weighted  sum 
of  these  transitions  (upper  curve)  plotted  with  the  experimental  data  for  sample  A363. 


The  final  results  of  the  Breit-Rabi  analysis  are  shown  in  Figure  62  along  with  the 
experimental  data  from  sample  A363  for  comparison.  Data  for  sample  SB0009B  shows  a 
precisely  identical  structure,  with  the  high  field  peaks  again  being  obscured  due  to  greater 
high  field  cavity  absorption. 

Model  parameters  for  both  samples  (A363  and  SB0009B)  were  found  to  be 
identical,  with  g  =  2.116+0.005,  =  O.OSdOTO.OOlcm'* ,  ™A  =  0.0674±0.001cm"' 

and  a  Gaussian  line  width  of  approximately  130  G.  The  ratio  of  the  isotopic  hyperfine 
splitting  constants  determined  above  =0.78?)  matches  the  ratio  of  Ga  isotopic 
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nuclear  g- values  g j g  =0.787)  and  indicates  that  the  splitting  observed  is  due  to 
hyperfine  coupling  with  a  Ga  atom.  Coupling  with  a  Ga  atom  could  indicate  any  of  the 
following  states: 

•  an  electron  bound  to  a  gallium  interstitial, 

•  a  bound  state  of  a  nitrogen  vacancy  coupling  with  gallium  nearest 
neighbors,  or 

•  an  electron  in  a  bonded  state  of  a  complex  involving  gallium. 

Fortunately,  literature  references  exist  that  have  measured  similar  paramagnetic  states  in 
related  materials  (GaP  and  AlGaAs)  by  EPR  and  in  GaN  by  related  techniques. 

The  similarity  of  the  hyperfine  coupling  constant  determined  in  this  study 
(®A  =  1588.9±30M//z)  with  that  reported  by  Kennedy  and  Spencer  in  AlGaAs  via  EPR 
measurements  ( =  1498±30M//z )  and  identified  as  a  Ga  interstitial  suggest  that  the 
signal  under  investigation  is  also  Ga  interstitial.  Eurthermore,  the  observed  hyperfine 
splitting  constant  is  in  excellent  agreement  with  values  found  via  ODMR  by  Einde,  et  al 
[94]  (®A  =  1580±50M//z)  and  Bozdog,  et  al  [95]  (®A  =  \600±50MHz)  in  electron 
irradiated  GaN  and  attributed  by  both  authors  to  an  interstitial  atom,  possibly 
complexed  with  another  defect  or  impurity,  due  to  similarities  with  previous  studies  in 
related  materials  (GaP  [96]  and  AlGaAs  [97]). 

An  additional  absorption  signal  overlaid  with  the  hyperfine  spectrum  is  noted  at 
g=2.2;  however,  no  hyperfine  structure  is  attributable  to  this  spectral  feature.  It  is 
possible  that  it  may  be  a  signal  arising  from  free  nitrogen  interstitials  or  vacancies,  which 
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would  be  expected  to  exhibit  a  much  more  closely  spaced  hyperfine  structure  than  for  Ga, 
possibly  producing  the  relatively  broad,  structureless  spectral  feature  observed. 


Figure  62.  Plot  A:  Sample  A363  widescan  data  (baseline  subtracted).  PlotB:  Breit-Rabi 
formula  model  with  A,®  =  0.0530±0.001cm  '  andg,,  =2.116+0.005  ,  summed  for  ^^Ga 
and  ^^Ga.  Vertical  arrows  are  added  to  aid  in  comparison  of  the  plots. 

Given  the  good  agreement  in  the  measured  hyperfine  coupling  constants  in  this 
study  and  in  previous  ODMR  work  in  GaN  [94], [95]  it  is  reasonable  to  conclude  that  the 
post-irradiation  defect  observed  here  is  of  the  same  nature  as  those  observed  via  ODMR. 
The  observation  of  the  Ga,  spectrum  following  irradiation  along  with  the  simultaneous 
passivation  of  the  SDi  spectral  line  is  a  further  indication  that  the  passivation  of  the 
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dominant  shallow  donor  is  due  primarily  to  the  action  of  mobile  nitrogen  interstitials, 
although  a  lack  of  quantitative  data  regarding  the  introduction  of  the  Ga  interstitial  makes 
this  assessment  speculative. 

When  viewed  in  light  of  the  observed  donor  complexation  process  discussed 
previously,  this  observation  is  particularly  exciting.  For  the  first  time,  a  means  exists  to 
directly  observe  the  differing  rates  of  Ga  and  N  sublattice  damage  simultaneously  in  a 
single  sample. 
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VII.  Conclusions 


Experimental  Conclusions 

The  controversial  origin  of  the  native  n-type  doping  in  GaN  is  not  due  to  defect 
doping  by  simple  point  defects,  such  as  the  nitrogen  vacancy.  The  primary  GaN  shallow 
donor  sites  with  activation  energy  of  29.6  meV  are  observed  to  decrease  in  density 
following  irradiation  in  HVPE  and  MBE  grown  GaN  samples  with  shallow  donor 
concentrations  from  4.8x10'^  cm'^  to  1.38x10'^  cm'^,  which  is  incompatible  with  a 
point  defect  source  for  the  native  shallow  donor.  Native  n-type  doping  in  GaN  thus  must 
be  impurity-driven,  rather  than  defect-driven. 

It  is  clear  from  these  irradiation  studies  that  the  shallow  donor  sites  in  GaN  are 
effectively  “passivated”  as  a  result  of  low  energy  electron  irradiation  in  moderate  doses. 
This  passivation  is  theorized  to  be  the  result  of  compensation  by  radiation-induced 
acceptor  states  (presumably  N  interstitials),  possibly  in  combination  with  complex 
formation  between  the  mobile  negatively  charged  nitrogen  interstitial  and  positively 
charged  shallow  donor  sites.  This  theory  is  supported  by: 

•  the  observed  dependence  of  the  passivation  on  the  initial  shallow  donor  density, 

•  the  observation  of  room  temperature  annealing  indicative  of  a  mobile  defect, 

•  the  observed  reduction  in  both  the  Hall  measured  shallow  donor  density  and 
accompanying  decrease  in  the  EPR  absorption  peak  identified  as  the  shallow 
donor  in  GaN, 
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•  and  the  simultaneous  observation  of  positively  identified  Ga  interstitials  in  large 
concentrations  in  these  samples. 

Compensation  and  complexation  are  presented  as  competing  processes,  both  of  which 
decrease  the  shallow  donor  concentration.  Donor-defect  complex  formation  without 
accompanying  acceptor  compensation  is  tentatively  rejected  as  an  explanation  for  the 
observed  donor  density  reduction  based  upon  the  dependence  of  the  reduction  on  initial 
donor  concentration.  Damage  constants  for  this  process  are  determined  to  be  between 
0.145  cm'^  and  0.209  cm'^  as  measured  by  Hall  and  EPR  respectively. 

First  observation  of  the  gallium  interstitial  hyperfine  spectrum  in  GaN  is  reported 
and  is  supported  by  modeling  of  the  observed  spectrum.  This  defect  is  observed 
following  high-dose  irradiation.  The  identification  of  the  gallium  interstitial  spectrum 
allows  direct  computation  of  the  defect  introduction  rate  in  electron  irradiated  GaN.  With 
this  observation,  radiation  induced  sublattice  damage  in  GaN  can  be  clearly  identified 
and  measured  in  single  samples  for  both  sublattices.  The  first  simultaneous  identification 
of  gallium  and  nitrogen  sublattice  damage  effects  verifies  the  prediction  that  gallium  and 
nitrogen  sublattice  damage  are  expected  at  1.0  MeV  irradiation  energies. 

Proposed  Experimental  Direction 

Several  avenues  of  research  immediately  suggest  themselves  as  a  consequence  of  the  data 
garnered  in  this  experimental  effort.  A  few  of  the  possibilities  of  highest  interest  are: 

•  Measurement  of  the  ratios  of  passivation  (Ni  induced)  vs  Gai  concentration  as 
functions  of  particle  energy,  dose,  and  angle  to  parametrize  sublattice  damage 
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constants.  Identification  of  defect  introduction  rates  on  each  sublattice  as  a 


function  of  the  radiation  characteristics  is  now  possible. 

•  An  experimental  investigation  of  radiation  effects  in  different  donor  density 
samples  with  an  emphasis  on  variations  in  the  irradiation  temperatures  and  the  use 
of  extremely  low  temperatures.  This  could  perhaps  reveal  the  relative 
contributions  of  acceptor  concentration  and  donor-defect  complex  formation  to 
the  donor  density  decrease. 

•  A  more  thorough  annealing  study  of  both  defects,  to  determine  sublattice 
annealing  rates  and  monitor  formation  of  intermediary  states 

•  A  more  thorough  analysis  of  the  Ga  interstitial  hyperfine  spectra  including 
performing  angular  measurements  to  identify  elements  of  the  hyperfine  tensor  and 
background  levels  in  virgin  materials. 

The  measurement  and  assessment  of  the  gallium  and  nitrogen  interstitial  effects  in 
irradiated  GaN  set  the  stage  for  effective,  detailed  analyses  of  defect  formation  and 
dynamics  in  GaN  through  the  coupled  use  of  the  EPR  and  Hall  experimental 
methodologies.  This  effort  should  be  pursued  in  order  to  better  understand  the 
underlying  physical  processes  governing  radiation  interactions  in  GaN  materials. 


157 


Bibliography 


1  Lester  Eastman  and  Umesh  Mishra,  IEEE  Spectrum  Online  (May  2002),  edited  by 
Glenn  Zorpette. 

2  E.  Shen,  et  al,  IEEE  Elec.  Dev.  Eett.  22(10),  457  (2001). 

3  Yi-Geng  Wu,  et  al,  IEEE  Elec.  Dev.  Eett.  48(3),  586  (2001). 

4  E.  Binet,  J.Y.  Duboz,  E.  Rosencher,  E.  Scholz,  and  V.  Harle,  Appl.  Phys.  Eett. 

69(2),  1202(1996). 

5  Thomas  Gehrke,  Manufacturing  GaN,  Nitronex  Corporation  publication. 

6  D.C.  Eook,  D.C.  Reynolds,  J.W.  Hemsky,  J.R.  Sizelove,  R.E.  Jones,  and  R.J. 
Molnar,  Phys.  Rev.  Eett.  79(12),  2273  (1997). 

7  Committee  on  the  Assessment  of  Combat  Hybrid  Power  Systems,  National 
Research  Council,  Combat  Hybrid  Power  System  Component  Technologies: 
Technical  Challenges  and  Research  Priorities  (2003). 

8  Triangle  TechJournal,  Nitronex  Closes  on  $1 1.3  Million  in  bunding,  August 
2003. 

9  D.C.  Eook,  Z.  Eang,  and  E.  Polenta,  MRS  Internet  J.  Nitride  Semicond.  Res.  5S1, 
WI0.5  (2000). 

10  Morko9,  Nitride  Semiconductors  and  Devices,  Springer  Verlag,  1999. 

1 1  Chin  et  al.  Journal  of  Applied  Physics,  75,  7365  (1994). 

12  Weitzel,  Wide  Bandgap  Semiconductor  RE  Power  Devices,  1997  IEEE 
Symposium  on  Compound  Semiconductors,  pg  479. 

13  Morkoc  et  al.  Prog.  Quant.  Electron.  20,  361  (1996). 

14  Bougrov  et  al.  Properties  of  Advanced  Semiconductor  Materials,  Wiley  2001. 

15  Suzuki  et  al,  Phys.  Rev.  B,  52  (11),  pp  8132  -  8139. 

16  Eanciulli  et  al,  Phy.  Rev.  B,  48  (20),  pp  15144-15147. 


158 


17  Suzuki  et  al,  Electron  and  Hole  Effective  masses  of  GaN  and  AIN,  in  Properties, 
Processing  and  Applications  of  Gallium  Nitride  and  Related  Semiconductors, 
INSPEC  1999.  pp  175-178. 

18  H.P.  Maruska  and  J.J.  Tietjen,  Appl.  Phys.  Eett.  15,  327  (1969). 

19  J.E.  Van  Nostrand,  J.  Solomon,  A.  Saxler,  Q.H.  Xie,  D.C.  Reynolds,  and  D.C. 
Eook,  J.  Appl.  Phys,  87(12),  8766  (2000). 

20  D.C.  Eook,  G.C.  Earlow,  P.J.  Drevinsky,  D.E.  Bliss  and  J.R.  Sizelove,  pending 
publication. 

21  Y.  Nakano  and  T.  Kachi,  Appl.  Phys.  Eett.  79(11),  1631  (2001). 

22  Orton  et  al.  Acceptors  in  GaN  and  related  compounds,  in  Properties,  Processing 
and  Applications  of  Gallium  Nitride  and  Related  Semiconductors,  INSPEC  1999. 
pp  98-100. 

23  Keller,  et  al,  IEEE  Transaction  on  Electron  Devices  48(3),  552. 

24  Ambacher,  et  al,  J.  Appl.  Phys.  85(6),  3222  (1999). 

25  Showalter,  BMDO  Spring  Update  2001  Internet  Edition, 
http://www.acq.osd.mil/bmdo/bmdolink/html/update/spr01/updhor3.htm. 

26  S.S.  Park,  I-W.  Park,  and  S.H.  Choh,  Jpn.  J.  Appl.  Phys.,  Part  2  39,  El  141  (2000). 

27  Jin-Cheng  Zheng,  Hui-Qiong  Wang,  A.T.S.  Wee  and  C.H.A.  Huan,  Appl.  Phys. 
Eett.  79(11),  1643  (2001). 

28  DenBaars,  SP.  ITRI,  Inc  Report:  “High  Temperature  Electronics  in  Japan”, 
sponsored  by  Office  of  Naval  Research  and  National  Science  Eoundation  (1998). 

29  S.C.  Jain,  M.  Willander,  J.  Nayaran,  and  R.  Van  Overstraeten,  J.  Appl.  Phys. 
87(3),  965  (2000). 

30  Nasser,  et  al.  Journal  of  Microwaves  and  Optoelectronics,  2(3),  22  (2001). 

31  Debasis  Sengupta,  J.  Phys.  Chem.  B  107(1),  297  (2000). 

32  G.B.  Stringfellow,  Organometallic  Vapor  Phase  Epitaxy:  Theory  and  Practice 
(Academic  Press,  San  Diego,  2002). 

33  Carlos,  et  al.  Phys.  Rev.  B,  48  (24),  17878  (1993). 


159 


34  D.C.  Look,  G.C.  Farlow,  PJ.  Drevinsky,  D.F.  Bliss,  and  J.R.  Sizelove,  pending 
publication  (2003). 

35  Nord,  et  al.  Nucl.  Inst,  and  Meth.  In  Phys.  Res.  B,  202,  93  (2003). 

36  Look,  DC.  Private  conversation,  12  April  2002. 

37  K.H.  Chow,  C.D.  Watkins,  V.  Harle,  and  F.  Scholz,  Phys.  Rev.  B  55  R10177 
(1997). 

38  D.C.  Look,  Phys.  Stat.  Sol.  B  228(1),  293  (2002). 

39  Atherton,  NM  Principles  of  Electron  Spin  Resonance,  edited  by  Ellis  Horwood 
and  T.J.  Kemp  (Prentice  Hall,  New  York,  1993),  pg  1. 

40  A.R.  Boate,  J.R.  Morton,  and  K.E.  Preston,  J.  Mag.  Res.  24,  259  (1976). 

41  B.  Bleaney,  in  Hyperfine  Interactions,  edited  by  Arthur  J.  Ereeman  and  Richard 
B.  Erankel  (Academic  Press,  New  York,  1967)  chap.  I,  pg  41. 

42  John  A.  Weil,  James  R.  Bolton,  and  John  E.  Wertz,  Electron  Paramagnetic 
Resonance:  Elementary  Theory  and  Applications  (Wiley-Interscience,  New  York, 
1994)  pg  386. 

43  B.  Eanghanki,  Electron  Paramagnetic  Resonance  of  Process  Induced  Defects  in 
Silicon,  dissertation,  Paderborn  University  (2001). 

44  M.W.  Bayerl,  M.S.  Brandt,  and  M.  Stutzmann,  Phys.  Stat.  Sol.  (a)  159(2),  R5 
(1997). 

45  John  A.  Weil,  James  R.  Bolton,  and  John  E.  Wertz,  Electron  Paramagnetic 
Resonance:  Elementary  Theory  and  Applications  (Wiley-Interscience,  New 
York,  1994). 

46  Charles  Poole,  Electron  Spin  Resonance:  A  Comprehensive  Treatise  on 
Experimental  Techniques  (Wiley,  New  York,  1983). 

47  Atherton,  NM  Principles  of  Electron  Spin  Resonance,  edited  by  Ellis  Horwood 
and  T.J.  Kemp  (Prentice  Hall,  New  York,  1993),  chap  4,  pg  130. 

48  John  A.  Weil,  James  R.  Bolton,  and  John  E.  Wertz,  Electron  Paramagnetic 
Resonance:  Elementary  Theory  and  Applications  (Wiley-Interscience,  New 
York,  1994)  App.  B,  pg  447. 

49  Atherton,  NM  Principles  of  Electron  Spin  Resonance,  edited  by  Ellis  Horwood 
and  T.J.  Kemp  (Prentice  Hall,  New  York,  1993),  chap  4,  pg  23. 


160 


50  Author  conversation  with  Dr.  Mary  Ellen  Zvanut  (August  2001). 

51  Glaser,  et  al.  Phys.  Rev.  B,  51  (19),  13326  (1995). 

52  M.  Linde,  S.J.  Uftring,  G.D.  Watkins,  V.  Harle,  and  F.  Scholz,  Phys.  Rev.  B,  55 
(16),R10177  (1997). 

53  Bruker  Spin  Resonance  Almanac  2003,  published  by  Bruker  Inc. 

54  John  A.  Weil,  James  R.  Bolton,  and  John  E.  Wertz,  Electron  Paramagnetic 
Resonance:  Elementary  Theory  and  Applications  (Wiley-Interscience,  New 
York,  1994)  Ch  3,  pg  73. 

55  Look  DC  Electrical  Characterization  of  GaAs  Materials  and  Devices  (Wiley, 
Chichester,  1989)  pp  1-2. 

56  M  Ahoujja,  Y.K.  Yeo,  R.L.  Hengehold,  and  J.E.  Van  Nostrand,  Appl.  Phys.  Lett. 
80(8),  1382  (2002). 

57  Charles  Kittel,  Introduction  to  Solid  State  Physics  (John  Wiley  and  Sons,  New 
York,  1966)  pg313. 

58  H.  Boudinov,  S.O.  Kucheyev,  J.S.  Williams,  C.  Jagadish,  and  G.  Li,  Appl.  Phys. 
Lett.  78(7),  943  (2001). 

59  A. I.  Titov  and  S.O.  Kucheyev,  J.  Appl.  Phys.  92(10),  5740  (2002). 

60  Atherton,  NM  Principles  of  Electron  Spin  Resonance,  edited  by  Ellis  Horwood 
and  T.J.  Kemp  (Prentice  Hall,  New  York,  1993),  chap  2,  pp  49-50. 

61  Cohen-Tannoudji,  Diu,  Laloe,  Quantum  Mechanics,  (Wiley-Interscience,  New 
York,  1977)  pp  643-653. 

62  G.  Breit  and  LI.  Rabi,  Phys.  Rev.  38,  2082  (1931). 

63  D.A.  Steck,  Rubidium  87  D  Line  Data,  Theoretical  Division  Los  Alamos  National 
Laboratory  (2001). 

64  Cohen-Tannoudji,  Diu,  Laloe,  Quantum  Mechanics,  (Wiley-Interscience,  New 
York,  1977)  pp  1035-1039. 

65  Atherton,  NM  Principles  of  Electron  Spin  Resonance,  edited  by  Ellis  Horwood 
and  T.J.  Kemp  (Prentice  Hall,  New  York,  1993),  chap  2,  pg  321. 


161 


66  Morko^  H.  Materials  Science  and  Engineering  R  33,  135-207  (2001). 

67  Look,  DC  and  Molnar,  RJ,  Appl.  Phys.  Lett.  70(25),  3377  (1997). 

68  Look,  DC  Electrical  Characterization  of  GaAs  Materials  and  Devices,  edited  by 
D.V.  Morgan  and  H.L.  Grubin  (Wiley,  Chichester,  1989),  pp  107-131. 

69  S.S.  Park,  I-W.  Park,  and  S.H.  Choh,  Jpn.  J.  Appl.  Phys.,  Part  2  39,  LI  141  (2000). 

70  Author  conversations  with  Dr.  Gary  Larlow,  Wright  State  University  (January 
2003). 

7 1  Integrated  TIGER  Series  of  Coupled  Electron/Photon  Monte  Carlo  Transport 
Codes  System,  Radiation  Safety  Information  Computational  Center,  Oak  Ridge 
National  Laboratory. 

72  C.A.  Bates,  J.P.  Bentley,  B.L.  Jones,  and  W.S.  Moore,  J.  Phys.  C:  Solid  State 
Phys.  3,  570(1970). 

73  Author  discussions  with  Dr.  Mary  Ellen  Zvanut  (August  2001). 

74  Author  discussions  with  Dr.  William  Carlos,  Naval  Research  Laboratory, 
Washington  DC  (October  2002). 

75  M.  Ahoujja,  Y.K.  Yeo,  R.L.  Hengehold,  and  J.E.  Van  Nostrand,  Appl.  Phys.  Lett. 
50(8),  1382  (2002). 

76  C.  Poole  Jr  and  H.  Larwell,  Theory  of  Magnetic  Resonance,  (Wiley,  New  York, 
1987). 

77  L.K.  Koschnick,  K.  Michael,  J-M.  Spaeth,  B.  Beaumont  and  P.  Gibart,  MRS 
Internet  J.  Nitride  Semicond.  Res.  1,  9  (1996). 

78  Author  discussions  with  Dr.  Bill  Carlos,  Naval  Research  Laboratory  Washington 
DC  (September  2002). 

79  Author  discussions  with  Dr.  David  Look,  Wright  State  University  Semiconductor 
Research  Center  (August  2003). 

80  C.  Mavroidis,  J.  Harris,  M.J.  Kappers,  C.J.  Humphries,  and  Z.  Bougrioua,  J.  Appl. 
Phys.  93(11),  9095  (2003). 

81  D.  Huang,  et  al.  Solid  State  Elec.  45,  711  (2001). 


162 


82  Look,  DC  Electrical  Characterization  of  GaAs  Materials  and  Devices,  edited  by 
D.V.  Morgan  and  H.L.  Grubin  (Wiley,  Chichester,  1989),  pp  77-78. 

83  Palczewska,  et  al.  MRS  Internet  J.  Nitride  Semicond.  Res.  3,  45  (1998). 

84  A.J.  Ptak,  et  al,  Appl.  Phys.  Lett.  79(17),  2740  (2001). 

85  M.  Ilegems  and  H.C.  Montgomery,  J.  Phys.  And  Chem.  Of  Sol.  34,  885  (1973). 

86  E.R.  Glaser,  T.A.  Kennedy,  H.C.  Crookham,  J.A.  Ereitas,  Jr,  M.  Asif  Khan,  D.  T. 
Olson,  and  J.N.  Kuznia,  Appl.  Phys.  Lett.  63(19),  2673  (1993). 

87  Bozdog  et  al,  Phys  Rev  B,  59(19)  12479  (1999). 

88  J.  Neugebauer  and  C.G.  Van  de  Walle,  Appl.  Phys.  Lett.  69,  503  (1996). 

89  K.D.  Greene,  J.R.  Petrosky,  and  D.C.  Look,  pending  publication. 

90  P.  Boguslawski,  E.L.  Briggs,  and  J.  Bernholc,  Phys.  Rev.  B  51,  17255  (1995). 

91  K.H.  Chow,  G.D.  Watkins,  Akira  Usui,  and  M.  Mizuta,  Phys.  Rev.  Lett.  85(13), 
2761  (2000). 

92  H.J.  von  Bardeleben  and  J.C.  Bourgoin,  Phys.  Rev.  B  33(4),  2890  (1986). 

93  H.J.  von  Bardeleben,  J.C.  Bourgoin,  and  A.  Miret,  Phys.  Rev.  B  34(2),  1360 
(1986). 

94  Matthias  Linde,  S  J  Uftring,  G  D  Watkins,  V  Harle,  and  L  Scholz,  Phys.  Rev.  B 
55(16),  R10177  (1997). 

95  C  Bozdog,  et  al.  Phys.  Rev.  B,  59(19),  12479  (1999). 

96  KM  Lee,  Defects  in  Electronic  Materials,  MRS  Symposia  Proceedings  No.  104, 
edited  by  M  Stavola,  SJ  Pearton,  and  G  Davies  (Materials  Research  Society, 
Pittsburgh,  1988). 

97  A  Kennedy  and  M  G  Spencer,  Phys.  Rev.  Lett.  57(21),  2690  (1986). 


163 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
0MB  No.  0704-0188 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering 
and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of 
information,  including  suggestions  for  reducing  this  burden  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite 
1204,  Arlington,  VA  22202-4302,  and  to  the  Office  of  Management  and  Budget,  Paperwork  Reduction  Project  (0704-0188),  Washington,  DC  20503. 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


2.  REPORT  TYPE  3.  DATES  COVERED 


1 .  REPORT  DATE  (DD-MM-YYYY) 

08-09-2003 


4.  TITLE  AND  SUBTITLE 


Doctoral  Dissertation 


Aug  2001  -  Sep  2003 


5a.  CONTRACT  NUMBER 


ELECTRON  PARAMAGNETIC  RESONANCE  SPECTROSCOPY  AND 
HALL  EFFECT  STUDIES  OF  THE  EFFECTS  OF  LOW  ENERGY 
ELECTRON  ENERGY  ON  GALLIUM  NITRIDE  (GAN) 


5b.  GRANT  NUMBER 


5c.  PROGRAM  ELEMENT  NUMBER 


6.  AUTHORISI 


5d.  PROJECT  NUMBER 


Greene,  Kevin,  D.,  Major,  USAF 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 


7.  PERFORMING  ORGANIZATION  NAMEISI  AND  ADDRESSIES) 

Air  Force  Institute  of  Technology 

Graduate  School  of  Engineering  and  Management  (AFIT/EN) 
2950  P  Street,  Building  640 

WPAFB  OH  45433-7765  DSN:  785-4539 


9.  SPONSORING/MONITORING  AGENCY  NAMEIS)  AND  ADDRESSIESI 

Defense  Threat  Reduction  Agency/CSNP 

Attn:  Mr.  Bill  Daitch 

John  J.  Kingman  Road 

Fort  Belvoir,  VA  22060-6201 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


AFIT/DSP/ENP/03-02 


10.  SPONSOR/MONITOR'S  ACRONYM(S) 


1 1 .  SPONSOR/MONITOR'S  REPORT 
NUMBERISI 


12.  DISTRIBUTION/AILABILITY  STATEMENT 

APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED 


14.  ABSTRACT 

The  nature  of  native  donors  in  GaN,  types  and  interactions  of  radiation-induced  defects,  and  damage  creation  coefficients  for  1.0 
MeV  electron  irradiation  have  been  ascertained  by  the  concerted  application  of  electron  paramagnetic  resonance  spectroscopy  and 
Ilall  effect  measurements  to  virgin  and  electron-irradiated  GaN  epilayers.  Samples  produced  via  molecular  beam  epitaxy  and  hydridi 
vapor  phase  epitaxy,  both  silicon  doped  and  nominally  undoped,  were  subjected  to  Van  de  Graff  generator  produced  monoenergtic 
electron  beams  with  total  fluences  of  10''16  -  10^18  electrons/cm^2.  Nitrogen  vacancies  are  rejected  as  a  possible  cause  of  n-type 
conductivity  in  nominally  undoped  GaN  due  to  the  decrease  of  shallow  donor  populations  following  irradiation.  Damage  constants 
for  a  donor-defect  complex  formation  or  defect  compensation  are  found  to  be  approximately  0.15  cm^-1.  Identification  by  resolved 
hyperfine  splitting  is  accomplished  for  radiation-induced  Ga  interstitial  complexes. 


15.  SUBJECT  TERMS 

Electron  Paramagnetic  Resonance,  Gallium  Nitrides,  Electron  Irradiation,  Semiconductors,  Hall  Effect,  Radiation  Damage,  III-V 
Compounds 


16.  SECURITY  CLASSIFICATION  OF:  | 

a.  REPORT 

b.  ABSTRACT 

c.  THIS  PAGE 

u 

u 

u 

17.  LIMITATION  OF 


19a.  NAME  OF  RESPONSIBLE  PERSON 

James  C.  Petrosky,  ETC,  USA  (ENP) 


19b.  TELEPHONE  NUMBER  linclude  area  code! 

(937)255-3636  x  4600 


Standard  Form  298  (Rev.  8/98) 

Prescribed  by  ANSI  Std.  239.18 


